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In technologically important materials, spin is more important than orbital motion 
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Chapter 1 


MAGNETIC FORCE 


The most fundamental aspect of magnetic 
fields is that forces and torques act on mag- 
nets and currents. The phenomenon of force 
and torque interactions between magnetic poles 
was investigated by Coulomb and Cavendish at 
the end of the 18th century. Oersted’s discov- 
ery of the magnetic effect of electric currents 
led Ampere and Faraday to invent the first mo- 
tors in the 1820’s. The experimental relations 
of forces between currents were determined and 
reported by Ampere. The results of those ex- 
periments are the basis of the notion of a mag- 
netic force field and its effect on magnetic poles 
and on electric currents. Magnetic force field, 
magnetic pole strength, and magnetic moment 
are quantities which are defined in this chapter. 


1.1 CURRENT FORCES 


Two parallel wires have a force of attraction 
which is proportional to the product of the two 
electric currents and inversely proportional to 
the separation of the wires. Figure 1.1 depicts 
two such currents flowing out of the paper. The 
force on the right-hand current, exerted over a 
length Lz is: 


we KB (1.1) 


where K is a constant depending on the unit 
system. For currents flowing in opposite direc- 
tions the force is repulsion, and when currents 
flow perpendicular to one another they exert 
no force on each other. Ampere demonstrated 
these effects and Maxwell formulated them by 
developing the notions of field theory. To deal 
with this action at a distance, Maxwell con- 
ceived of a force field which would act on a cur- 
rent dependent upon the current direction as 


vg (ovt) 


Figure 1.1: current forces 


well as the length of a conductor which carries 
the current. 

Taking Ampere’s experimental results, con- 
cisely stated in equation 1.1, and predicating 
a force field which also deals with antiparallel 
currents and perpendicular currents was not a 
simple task. To do this Maxwell invented vec- 
tors. Briefly, to account for Ampere’s observa- 
tions, the field from the first current had to be 
perpendicular to the second, and the force per 
unit length on the second current had to be a 
cross-product of the first current’s field and the 
second current. 

Br is the symbol which is used here to des- 
ignate the force field, a vector. The force field 
arising from current flowing in a long straight 
wire, as J, in figure 1.1, is perpendicular to 
both the current and to the radius vector from 
the wire to the point where the field is acting. 
This field is constant at a particular radius, di- 
rected in the right-hand sense around current 
I,: Pointing the right-hand thumb in the direc- 
tion of the source current (J), the curled fin- 
gers indicate the direction of the resulting field 
(counter-clock wise). 

In the MKS unit system, the force field acting 
on current 2 because of current 1 is Br,,: 


a Bol 
Br,, = Org tt x 1,,,) (1.2) 
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r 
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1, (vt 
Figure 1.2: magnetic force field 


where ys, 1a the permeability of free space, yu, = 
4x x 10-7, and the boldfaced 1 indicates a unit 
vector. The unit vector 1;, is in the direction 
of the current J;, and the unit vector 1,,, is 
the radial direction from current J; toward the 
field point. The direction of the force field is 
then the vector cross-product of the direction 
of the source current and the direction from the 
source current to the field point. 

The units of the force field is newtons per 
amp-meter, so the the constant K of equation 
1.1 can be seen to be p,/27 in the MKS unit 
system. generalizing from this example, the el- 
ement of force on an element of length dl car- 
rying a current J in a magnetic force field Br 
[Lorentz’s equation] is: 


d¥ = Id1x Br (1.3) 


1.2 MAGNETS 


Magnets are attracted to iron products such as 
refrigerator doors. Using a pair of magnets, 
you can easily demonstrate that they can at- 
tract and repel one another. This is explained 
by magnetic poles: a magnet has two poles of 
opposite polarity, and with a pair of magnets 
the poles that are alike repel one another while 
unlike poles attract each other. 
- Directions on the earth’s surface are found 
with a compass, which points toward the geo- 
graphical north pole, or more accurately toward 
the magnetic pole in the north. The compass is 
a needle magnet on a pivot. The end of the nee- 
dle which points to the north is defined as the 
magnet’s north pole. This leads us to conclude 
that the earth’s magnetic pole in the north is 
actually a south pole, as it attracts the com- 
pass’s north pole. | 

Attraction and repulsion are forces which act 
on pairs of magnets: Two north poles repel each 
other equally, and two south poles repel each 
other equally. The north pole of one magnet 


R.B. Yarbrough 


© 15 May 1990 2 


and the south pole of another attract each other 
equally. The north and south poles of the same 
magnet exert forces, but they can only cause de- 
formation of the magnet, not cause motion. For 
8 two magnet system, the result is two forces 
acting on each pole of each magnet. See fig- 
ure 1.3 where the pole forces are shown acting 
equally and oppositely on each pair of poles as 
single vectors, with summed forces indicated on 
each pole. The individual forces are inversely 
proportional to the square of the separation. 


Fs, 


Ss Magnet 1 WN 


bead 


tN 
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Figure 1.3: magnet forces 


Suppose two identical magnets have pole 
strengths P and lengths |. Placed on a line 
with north poles adjacent and and separated by 
a distance d as in figure 1.4. Coulomb’s exper- 
iments showed that the forces between two like 
magetic poles is a repulsion proportional to the 
product of the pole strengths and inversely pro- 
portional to the separation between the poles. 


exgrce Fas 
Figure 1.4 


Then, using Coulomb’s Law, the forces indi- 
cated are: 


Fyn = kP?/d? 
Fsn = Fns — —kP?/(d+ 1)? 
Fss = kP?/(d+2l)? 
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The total force (to the right) on the right-hand 
magnet can be calculated from this expression: 
F = Fnn+Fsn+Fns+Fss 


P2 P2 P2 
k—— — 2k-———— + k-—___ 
( + yap 


d2 d+l) 


and the force on the left-hand magnet is equal 
but to the left. The constant k depends on 
the unit system and the definition of the pole 
strength P. The notion of force acting at a 
distance from the source of the force is called 
a field, and in particular a force field. Here 
with two magnets we find such a phenomenon. 
Each magnet experiences force because of the 
other magnet. Then each magnet experiences 
the other as a source of a force field. Each pole 
can be visualized as having the force acting on 
it as the product of the force field at that point 
and its own pole strength P. Because this is a 
magnetic force, it is of the same nature as the 
field caused by an electric current, Br. From 
figure 1.4, the north pole of the right hand mag- 
net would experience a force field of to the right 
of magnitude: 


Bp =kP/d? —kP/(d+1)? 


where the first term on the right-hand side is 
due to the left- hand magnet’s north pole, and 
the second is due to the left-hand magnet’s 
south pole. 

Thus the field due to an isolated magnetic 
pole is directed from that pole toward the point 
where the field is being calculated. The field is 
proportional to the pole strength and inversely 
proportional to the square of the distance from 
the pole to the point. We define the vector from 
the source (pole) to the field point as r,y and 
its magnitude as r,y. The vector field at the 
field point is Br: 


Bri=k 


3 73 (1.4) 


~ 


and the unit vector is then seen to be: 


Loy = rsy/|rsy| 


For a group of poles, the field is calculated 
by vectorally adding the effects of each. 

The force field due to a magnet acts on other 
magnets as well as on electric currents. This 
requires pole strength be defined to result in 
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correctly calculated forces on current carrying 
conductors. Dimensionally the force field Bp is 
in newtons per amp- meter in eq. 1.2, and must 
be in units of newtons per unit pole in equation 
1.4. Thus it is necessary that the P be in units 
of amp-meter. From eq. 1.4 it is necessary that 
the units of kP be henrys x amps to result in 
a force field in Teslas or Webers/m?. Then k 
must be henries per meter, and so must contain © 
to Which has the correct units. 


1.3 MAGNETS AND CURRENT 
LOOPS 


The field generated by a long straight wire car- 
rying a current is quite different from that of 
a magnet. In order to obtain the correct con- 
stants in the expression for the force field due 
to a magnetic pole, we must find an equiva- 
lency with current. The equivalency is between 
the far field of a magnet and that of a loop of 
current. The constant k is determined by com- 
paring the fields generated by each. We begin 
with the field due to a magnet. 


1.3.1 FIELD OF A MAGNET 


Figure 1.5: magnet field 


A bar magnet consists of a north pole +P 
(a positive magnetic charge) and a south pole 
—P (negative magnetic charge) separated by 
a length L. To determine the magnetic field 
caused by a bar magnet or “dipole”, we use a 
form of Coulomb’s law. The coordinate sys- 
tem most convenient for this situation is the 
spherical systern which uses the radial vector r, 
measured from the origin - which will be the 
middle of the bar magnet; an azimuthal an- 
gle, ¢, which is measured from an arbitrary line 
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through the center of the magnet, but perpen- 
dicular to the axis joining the two poles; and a 
polar angle 8, measured from the half of the 
polar axis passing through the north pole of 
the magnet. The coordinates for our analysis 
is shown in fig. 1.5. the north pole is on the 
polar axis at r = L/2, 6 = 0 and the south pole 
is at r = L/2, 8 = 180°. The application of 
Coulomb’s Law, obtains the following expres- 
sion for the field at (r,6,¢): 


kP1,(r- 5 cos 6) 


(r2 + © — rL cos 6)3/? 
kPl.(r+ s cos 6) 
(r24+ £ + rL cos 6)3/? 
kP14(4 sin 6) 
(r2 + £ — rL cos 6)3/? 
kP1,(4 sin 0) 


(r24+ © + rL cos 6)3/2 


Br = 


The far field approximation is that r >> L. 
Then, making the approximation that 


(l1+2)-"/?x1-—nz/2 
we have the far field expression: 


2 cos 6 sin @ 
me 1, + 16} 


Br =kPL{ (1.5) 


1.3.2 FIELD OF A CURRENT LOOP 


The field of a loop of current and the field of a 
bar magnet become identical when the distance 
from the source becomes large with respect to 
the dimensions of the source, i.e. L for the bar 
magnet (r >> L) and the radius for the cur- 
rent loop. The analysis of the current loop, 


Figure 1.6: current loop field 


which is centered at the origin of the spheri- 
cal coordinate system, with the current on a 
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circle around the polar axis circulating in the 
right-hand sense, as indicated in figure 1.6, is 
most readily solved using the vector magnetic 
potential formulation of the Biot-Savart Law, 
and then taking the curl to obtain the force 
field. This is a common example in many text 
books on electromagnetic fields and is beyond 
the scope of this course. The results, when the 
distance from the source ris much greater than 
a, the radius of the current loop, the force field 
is found to be: 

sin 0) 


2 cos 6 
4nr3 ar 4nr3 16} 


which again is a far-field expression. 


Br = pona’l{ (1.6) 


1.3.3 MAGNETIC MOMENT 


A comparison of equations 1.5 and 1.6 shows 
that equivalency is between (1) the product of 
the bar magnet’s pole strength and pole sepa- 
ration, PL, and (2) the current loops product 
of current and area of the loop, mr?J. These 
terms are both called the magnetsc moment: 


MAGNETIC MOMENTS: %a’I equivalent to PL 


The magnetic moment is a vector, which for 
a bar-magnet points in the direction from the 
south pole (—P) toward the north pole (+P), 
while for the current loop the direction is found 
from the right hand rule applied to the current: 
cur] the fingers of the right hand in the direction 
of the current, and the extended thumb points 
in the direction of the magnetic moment. This 
is equivalent to crossing the radius direction 
with the direction of current. The magnetic 
moment is symbolized as m, and 1s defined as: 


Pligs. 00) 
IAlyzy (1.8) 
Again comparing equations 1.5 and 1.6, it is 
seen that the constant k for magnetic pole 


forces and fields in coulomb’s law equations 
must be : 


BAR MAGNET: mM = 
CURRENT LOOP: m = 


k= - IN THE MKS SYSTEM 


Then the magnetic force field due to an isolated 
magnetic pole, for the MKS unit system (from 
equation 1.4) is: 


(1.9) 


1. MAGNETIC FORCE 


1.4 TORQUE ON A MAGNETIC 
MOMENT 


Calculation of torque on a magnetic dipole can 
be done using either a magnet with Coulomb’s 
Law, or a current loop with equation 1.3. In 
either case, the torque is found to be: 

TORQUE T=m~x Br (1.10) 
Exercise 1.1: derive equation 1.10 for a bar 
magnet. 


Exercise 1.2: Derive equation 1.2 for a square 
current loop with an area a? and a current I. 


1.5 ELECTROMAGNETS 


The term magnet refers to permanent magnets 
that continue the same behavior over a long pe- 
riod of time. On the other hand, it is possible 
to magnetize materials such as iron and steel so 
that they become temporary magnets. An iron 
rod can be magnetized by applying an electric 
current to a coil of wire wound around it. It 
will be strongly magnetic as long as the current 
continues, but will lose most of its magnetism 
when the current is stopped. Such magnets are 
called electromagnets and their strength is pro- 
portional to the current. After the current has 
stopped they may have a small remnant mag- 
netism. These effects can be measured through 
forces acting on an electromagnet in the vicinity 
of another magnet. 


1.6 SUMMARY 


The definitions given here for magnetic pole 
strength and for magnetic moment have been 
made so the force on a pole and the torque on 
the magnetic moment are both calculated from 
the force field Br. The units for magnetic pole 
strength are amp-meters and the units for mag- 
netic moment are amp-meters2. The only unit 
system employed has been the MKS unit sys- 
tem. The units of the force field are directly 
determined as from the force equations 

The magnetic force field field, Br, has been 
defined for the following cases: 
1. Where it arises from a current in a long 
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straight wire (Ampere’s Law, from eq. 1.2): 


_ Mol; 
Br,, = Dar,g pir 


2. from an isolated magnetic pole (Coulomb’s 
Law, from eq. 1.7: 


3. Where it is the far-field arising from a mag- 
netic dipole (from eq. 1.5 and 1.6): 


Br= 


mons {21, cos 6 + le sin 4) } (1.11) 


4 
The force on an element of electric current is 
part of the Lorentz force equation as given in 
eq. 1.3: 
dF = Idl x Br 


The torque on a magnetic dipole moment m 
is given in eq. 1.10: 


T=mx Br 


Chapter 2 


MAGNETIC INDUCTION 


2.1 EXPERIMENTAL BASIS 


The statement that lke poles of magnets re- 
pel each other and unlike poles attract can be 
proven from the fact that it is possible to find 
three poles which mutually repel each other, 
but it is not possible to find three poles which 
mutually attract each other. 

From the earliest times magnets were known 
for their ability to attract (but not repel) cer- 
tain metals, most notably iron. This shows that 
a magnetic pole induces an unlike pole in the 
closest part of nearby iron. As magnetic poles 
occur only in pairs (i.e. magnetic monopoles 
have not been found), then a magnetic pole will 
also induce a like pole in the most remote part 
of the nearby iron. As these induced poles are 
equal in strength, the result will be an attrac- 
tion of the iron to the closer pole of the magnet. 

Oersted’s discovery of the magnetic effect of 
electric currents eventually led Faraday to in- 
vestigate whether magnetic fields could cause 
or induce electric current in a conductor. He 
found this could be done either by changing 
the current that created the field or by mov- 
ing a conductor relative to a magnet or other 
field source. Henry had independently discov- 
ered the effect of the changing current. Fara- 
day carefully described his experiments, includ- 
ing all the failures. This record keeping proved 
invaluable to Maxwell when he set out to un- 
derstand electromagnetism. 


It was in 1831 that Faraday discovered that 
electric current was induced in a wire which 
formed a closed circuit when the current in a 
nearby wire was changed. Faraday saw the in- 
duced current was the response to an electromo- 
tive force (abbreviated emf). Ohm’s discovery 
of the law of resistance in 1827 provided the key 


to understanding the relation between current 
and emf. 

Faraday attributed the induced emf to a flux 
arising from the magnetic field of the current 
in the nearby wire. He envisioned the induced 
emf being caused by changing this flux. The 
induced current was proportional to the emf, 
by Ohm’s Law. | 

The emf is equal to the negative time rate of 
change of the magnetic ” flux”: 


emf= -4 (2.1) 


The symbol ¢ (phi) indicates the quantity mag- 
netic flux. The emf has units of electric po- 
tential, volts in the MKS unit system. Thus 
the units of the magnetic flux ¢ is volt-seconds, 
which is also called webers. 

Equation 2.1 is the most fundamental form 
of Faraday’s Law of electromagnetic induction 
and serves as the definition for magnetic flux. 


2.2 MAGNETIC FLUX DENSITY 


Magnetic flux is a flow concept which involves 
no tangible material. However it is useful to 
think of magnetic flux as analogous to the flow 
of water, which has a particular direction at a 
particular point. At such a point the direction 
and speed of the water, combined with the den- 
sity of the water, can be used to define a flow 
density of the water. This flow density is a vec- 
tor, and could also be called the flux density of 
the water at that point. 

The flow of a stream of water is usually de- 
fined as the rate at which water passes through 
an imaginary plane or other surface. This 1s 
also the average flow density or flux density 
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multiplied by the area of the imaginary surface, 
provided that the surface is perpendicular to 
the flow density at every point. This can be ex- 
pressed as a surface integral, which involves the 
scalar product (dot product) of two vectors: the 
flux density and the surface element’s normal. 

Magnetic flux density is similarly defined in 
terms of magnetic flux. Using the symbol B 
for the vector magnetic flux density, the flux 
through a surface S is 


6= | [ B-1,ds 


where 1,, is the unit vector normal to the ele 
ment of surface area dS. 

The MKS units of magnetic flux density are 
seen to be volt-seconds per square meter, or 
webers per square meter which are also called 
teslas. 

Magnetic flux density B is a vector quan- 
tity which is defined at a so-called macroscopic 
point. A macroscopic point is not the same as 
a mathematical point, but should be thought 
of as a very small volume or surface area, yet 
large enough to accommodate a large number 
of atoms. 

Faraday’s law of electromagnetic induction 
can then be written in terms of the flux den- 


sity as: 
emf=—= ff Ba dS 
=e e 


This can further be expressed using the circu- 
lation of the electric field intensity E to replace 
the emf, i.e. the line integral of the electric 
field intensity, E, around the boundary of the 
surface S in the counter-clockwise direction 1s 
the emf. 


(2.2) 


(2.3) 


me E-a=-5 [[ B-1,4s (2.4) 
Cs § 


where the curve Cg is the boundary of the sur- 
face S. 


2.3 FLUX INTO A VOLUME 


In equation 2.4 the circulation of E is the emf 
around any closed path, and the surface over 
which B is integrated is any surface which is 
bounded by the circulation path. This leads 
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to an important law, which was formulated by 
Gauss together with a more general theorem of 
vector calculus. 


Consider a square box as shown in fig. 2.1. 
We take the boundary of side A as the path for 
the circulation of E, which encloses the surface 
A. The circulation is —d¢/dt passing through 
surface A. Notice also that the path of circula- 
tion also is the boundary of a surface composed 
of the box excluding side A. The circulation 
of E is also —d¢/dit through that complex sur- 
face. The conclusion to be drawn here is that 
the time rate of change of the flux into a vol- 
ume is identical to time rate of change of the 
flux out of the volume: there is no net time 
rate of change of flux into (or out of) a volume. 
Thus, in the absence of magnetic monopoles, 
the flux into a volume is equal to the flux out. 


figure 2.1 


It is useful to shrink the volume until it ap- 
proaches the macroscopic point in size. In do- 
ing this, it is a straight-forward process to show 
that the divergence of the flux density is zero. 


div B= 0 (2.5) 


exercise 2.1: Using cartesian coordinates, and 
taking limits as the dimensions of the box in 
figure 2.1 shrink to zero, using the definition of 
the partial differential and the fact that the net 
flux into the box is zero, demonstrate that the 
divergence of B is zero. 


Equation 2.5 tells us that magnetic flux is 
continuous, having no sources or sinks (in the 
absence of monopoles). This requires that mag- 
netic flux ”closes” on itself, and for that reason 
is called solenoidal. | 


2. MAGNETIC INDUCTION 


2.4 MOTIONAL EMF 


A conductor moving in a magnetic field is a rel- 
ativistic situation, and is equally viewed as a 
magnetic field moving with respect to the con- 
ductor. The charges within the conductor have 
a velocity relative to the magnetic field, which 
results in an effective electric field acting on the 
charges. If the charges are free to move, as in 
certain closed circuits, they will participate in 
a current flow. The force per unit charge is, by 
definition, the electric field intensity. The en- 
tire situation is summed up in the Lorentz force 
law: 


F=Q(E+vx Br) (2.6) 


here Q is the charge, E is the electric field in- 
tensity, v is the relative velocity of the charge 
with respect to the magnetic force field Br. 

In a situation where there is a segment of 
Wire moving in a magnetic field, and there is 
no provision for current flow, then an equilib- 
rium condition will arise according to equation 
2.6, where the net force will be zero, so that the 
electrostatic field arising from free charge distri- 


bution within the wire and the electromagnetic 


field due to the relative velocity of the charges 
in the magnetic field are balanced. In that case, 
the electric field must be: 
E=-vx Br (2.7) 
The emf generated over a length of conductor 


is found by the line integral of the electric field 
over the length: 


sj Se [o xBr)-dl (2.8) 

Equation 2.8 is fundamental to electrical 
power generation, and provides the relation be- 
tween emf and motion of a conductor in a mag- 
netic force field. In principal such a force field 
can be calculated from known current distri- 
butions and the resulting emf generated on a 
conductor moving in the force field also calcu- 
lated. Such a calculation will agree with an 
experimental measurement, showing that the 
force field is directly involved with the motional 
generation of emf. 

A loop of conductor can be placed in this 
calculated force field and the flux of the force 
field passing through the loop determined by 
integration. This flux will be a linear function 
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of the current. Thus the flux of the force field 
can be varied by varying the current. In doing 
so, the time rate of change of the flux due to the 
force field is found experimentally to be equal 
to the measured emfon such a loop. So, for the 
case of a gaseous media the magnetic force field 
Br and the magnetic flux density B are found 
to be equal. 


The magnetic force field Br acts on electric 
currents and on magnetic poles, but the results 
are measurable only in a non-solid environment. 
On the other hand, magnetic flux in a solid ma- 
terial is measurable as the emf in a conductor 
wrapped around the solid can be measured. As 
seen in equation 2.5, the divergence of the mag- 
netic flux density is zero, which means that the 
normal component of magnetic flux is continu- 
ous in crossing from one medium to another. 


Except for ferrofluids, a non-solid environ- 
ment is also non-magnetic. In non-magnetic 
environments experiments show that the mag- 
netic force field and the magnetic flux den- 
sity are identical. Inside solid magnetic envi- 
ronments measurements of forces and torques 
are impossible. However, one experiment in 
a ferrofluid indicates that the magnetic force 
field and the magnetic flux density are different 
in a magnetic material (environment)!. This 
experiment appears to support Maxwell’s un- 


‘derstanding of electromagnetic forces and emf. 


Maxwell seems to have espoused the notion that 
the magnetic force field depends only on cur- 
rents and magnetic poles. Some later workers 
have found this an unacceptable idea, and have 
insisted that the medium or environment plays 
a key role. These dissident workers have suc- 
ceeded in establishing their opinion in the most 
commonly used version of the MKS unit sys- 
tem, which is known as the Sommerfeld Sys- 
tem. This will be discussed in some detail at a 
later point. 


In order to follow Maxwell’s understanding, 
we first establish the concept of the magnetic 
field intensity, H, which is exactly proportional 
to the force field Br. 


1 Experimental demonstration that the couple on a 
bar magnet depends on H, not B, R.W.Whitworth and 
H.V Stopes-Roe, Nature, v 234, pp31-33, Nov 1971. 
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2.5 MAGNETIC FIELD INTEN- 
SITY 


The magnetic field intensity H is a vector, and, 
following Maxwell’s ideas, is calculated from 
the combination of current distribution and 
magnetic pole distribution. It is proportional to 
the vector we have been calling the force field, 
with the definition: 


Br = 40H (2.9) 


The general expression for the calculation of H 
is obtained by modifying the Ampere Law and 
Coulomb Law expressions of chapter 1 (equa- 
tions 1.2 and 1.7). The current and its direc- 
tion are combined in the form of current density 
enclosed in a volume element as indicated in fig- 
ure 2.2, which then requires an integration over 
all current density. this results in an expres- 
sion which is usually called the law of Biot and 
Savart: 


Jx1,y 
Hm = ff fara ave! 


where J is the current density (source point), 
r,s is the distance from the source point to the 
field point (the field point is that where the field 
is being calculated), and 1,, is the unit vector 
from the source point toward the field point. 
Because magnetic materials are typically mag- 
netized using electric currents, this part of the 
magnetic field H,, will be called the magnetiz- 
ing field. 


ERY toh 


Ca 


(2.10) 


fig. 2.2 Biot-Savart field source element 


Coulomb’s Law is modified to deal with pole 
density rather than isolated poles, which also 
results in a contribution to the field intensity 
requiring integration over all pole density: 


H.= [| | pm™! Vol 
Vol 4nr?, 


(2.11) 
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where pr, is the magnetic pole density (source 
point), 1,, is again the unit vector from the 
source point to the field point, and r,, is the 
distance from the source point to the field point. 
Because inside a magnetic material this part of 
the field is in in opposition to the magnetiza- 
tion, Hg will be called the demagnetizing field. 


fig. 2.3 pole density field source element 


Then the complete expression for magnetic 
field intensity is the sum of the magnetizing field 
and the demagnetizing field H,, + Ha: 


If TXT a Vol + 
Vol 4nr?, 


‘ote 
[[[faitera ean 
2.6 SUMMARY 
Faraday’s law of electromagnetic induction is: 
_ 9 
emf= ah 


In the MKS unit system, magnetic flux is in 
webers while magnetic flux density is in teslas. 

The divergence of the magnetic flux density 
is zero: div B = 0 

It has been shown that for the case of a 
gaseous media the magnetic force field Br and 
the magnetic flux density B are equal. 

As magnetic forces have been shown to be in- 
dependent of the medium (as was envisioned by 
Maxwell), the magnetic field intensity is defined 
as proportional to the quantity we have called 
the force field, as in equation 2.9: Br = pH. 


Chapter 3 


MAGNETIZATION AND POLE DENSITY 


3.1 INTRODUCTION 


The magnetic moment of a dipole, arising either 


from a bar magnet or a current loop has been 
defined in chapter 1, and found to have MKS 
units of amp-meters’. All atoms have magnetic 
moments which are attributed to the orbiting 
electrons and to the spinning of those orbiting 
electrons. In most materials the average mag- 
netic moment tends to align with an applied 
magnetic field. This behavior is called para- 
magnetism. In some materials, the magnetic 
moments tend to align in opposition to an ap- 
plied magnetic field. This behavior is called 
diamagnetism. Both paramagnetism and dia- 
magnetism are weak effects. 


Some materials have strong magnetic effects. 
To some degree, all such materials exhibit spon- 
taneous magnetization, the property observed 
in permanent magnets, in that they behave 
magnetically in the absence of an applied mag- 
netic field. These materials are broadly classi- 
fied as magnetic. Other terms such as ferromag- 
netic, ferrimagnetic and anti-ferromagnetic, are 
used describe specific materials, but we will use 
two different classifications of magnetic materi- 
als, one of which is soft, and the other as hard. 


The atomic magnetic moments in soft mag- 
netic materials tend to align with an applied 
magnetic field rather easily. Also, soft mag- 
netic materials have little remnant magnetiza- 
tion when the field is removed. Such materials 
are used in magnetic recording heads. 


Hard magnetic materials are more difficult to 
magnetize, but retain the alignment of mag- 
netic moments after the magnetizing field is re- 
moved. Such materials are used in magnetic 
recording materials. 


10 


3.2 FARADAY VOLTAGE AND 
FLUX MEASUREMENT 


The most common method of measuring mag- 
netic flux is via the Faraday Law relation of 
equation 2.1: 

d¢ 


dt 

which involves the measurement of the voltage 
generated in a loop of wire. The ends of the 
wire are not quite touching, so that the voltage 
is the instantaneous value of the time rate of 
change of the flux passing through the loop. To 
obtain the flux, that voltage is electronically 
integrated as indicated in figure 3.1. 


emf= — 


figure 3.1 


The magnetic flux density is a quantity which 
is closely related to the magnetization of a mag- 
netic material, and the flux density is approx- 
imated by dividing the flux by the area of the 
loop (which actually gives the average flux den- 
sity): 


Bavy = -4 [ems at (3.1) 


In some cases, in order to increase the voltage 
available, several turns of wire are formed into 


3. MAGNETIZATION AND POLE DENSITY 


a coil of N turns. The resulting voltage is then 
seen to be: 


d¢ 
Veo! = —emf = N— (3.2) 
Then the average flux density can be calculated 


in terms of the voltage from the N-turn coil: 


Bavg = ns al —— dt 


(> (3.3) 


3.3 MAGNETIZATION 


Magnetization ss defined as the average mag- 
netic moment per unit volume at a macroscopic 
point. Recall that a macroscopic point is a vec- 
tor point, one which can contain a large number 
of atoms. The magnetization is a vector quan- 
tity, which is symbolized with M in the MKS 
unit system: 


M= jj >.m 


~ Voino Vol (3.4) 


* [amps/m] 


where the units of the individual magnetic mo- 
ments (m,) are amp-meters?. 

The units of magnetization, M, are here de- 
fined to be the same as magnetic field intensity 
H, namely amps per meter. Within a magnetic 
material, as can be demonstrated experimen- 
tally using Faraday’s law, the vector sum of the 
magnetization and the magnetic field intensity 
determine the flux density: 


B = po(M + H) 


It is recalled that H consists of two parts, the 
magnetizing field due to currents, and the de- 
magnetizing field due to magnetic poles: H = 
H,, + Hy. Then we can write: 
B= Ho(M +H, + Hz) [teslas] (3.5) 
The demagnetizing field arises from magnetic 
moments within the magnetic material. This 
can be understood from the properties of the 
magnetic flux density, B, which must be con- 
tinuous on a macroscopic field basis. 


3.4 MAGNETIC POLE DENSITY 


A fundamental decision that Maxwell made in 
the formulation of electromagnetism was in the 
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definition of magnetic flux density. This is ex- 
pressed as in eq. 2.5: 


divB=0 
Applying this to equation 3.5, we obtain that: 
div(u.(M + Hm +Ha)) =0 


where p, is a constant. In the quasi-static sit- 
uation (i.e. no macroscopic eddy currents) the 
divergence of the magnetizing field (which has 
its sources outside of the magnetic material), is 
zero in the magnetic material., thus: 


div M + div Hy = 0 (3.6) 


We can use this result to calculate the demag- 
netizing field arising from a known distribution 
of magnetization. To develop this ability, we 
first integrate equation 3.6 over a volume, sep- 
arating the two parts: 


// div Hg dVol= -|f div M dVol 
Vol Vol 


The divergence theorem can be applied to the 
left side of this equation to obtain: 


[J div Hg dvol= | H, -1,dS, 
Vol Sy 


where S, is the surface of the volume and 1, 
is the outward pointing normal to the surface. 
Then, substituting this result into the previous 
equation: 


If Hz -1,dS, = — // div M dVol 
Vol 


(3.7) 
Which is a form of Gauss’s law. This shows 


that the divergence of the magnetization plays 
the part of a pole density: 


where py, is the magnetic pole density. Then 
equation 3.7 can be written: 


/ Hy -1ndS, = / / pm dVol (3.9) 
S» Vol 


We can further develop this by taking the 
case where there is simply a small volume, 
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AVol, and we take the surface for integration as 
a sphere centered on that small volume. Then 
on the surface the flux due to the small vol- 
ume is perpendicular to the surface, and of the 
same value at any point on the surface. The 
surface integral then becomes simply 4xr?Hg, 
with a vector direction from the small volume 
to the surface point. Then equation 3.9 can be 
re-written to take in all pole density as: 


mls 
H.= [ff Prost aVol 
Vol dary, 


which matches the demagnetizing field portion 
of equation 2.12. 


(3.10) 


3.5 FIELD MEASUREMENT 


A fundamental relation for magnetic materials 
is the B-H or M-H chaaracteristic which is typ- 
ically a hysteresis loop. Measurements for a 
B-H loop typically use the integration scheme 
describe with figure 3.1 to determine the flux. 
The measurement of the magnetic field in a 
magnetic material is impractical, because force 
measurements are not generally possible. The 
strategy employed is to calculate the magne- 
tizing field from the current configuration, and 
either calculate the demagnetizing field or else 
use a sample configuration where the demag- 
netizing field is negligible, the latter practice 
being preferred. Here we consider three current 
configurations to give known field values. In 
each of these systems, a Faraday voltage pickup 
coil will have mutual inductive linkage with the 
magnetizing field source, so that such a volt- 
age will be induced by a varying magnetizing 
field even in the absence of a magnetic mate- 
rial sample. This effect requires cancellation by 
some instumentation means, or correction by 
computation. 


3.5.1 THE HELMHOLTZ COIL PAIR 


A Helmholtz coil pair is a system of two identi- 
cal coils, each having N turns and carrying the 
same current I. The coils are on a common axis 
separated by the radius, r. This results in an 
axial field which is quite uniform on the axis 
at the midpoint between the two coils. This is 
illustrated in figure 3.2. 
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figure 3.2 


This system is used for small samples of ma- 
terial which require relatively small magnetiz- 
ing fields. The field at the center of the system 
1S: 


r?-NI 
[r2 + d2/4]3/2 


= 0.7155 [a/m] 

(3.11) 
The region of uniformity is a volume of about 
10% of the volume of the system. 


Hasis = 


3.5.2 THE LONG SOLENOID 


A cylindrically wound coil of wire is called a 
solenoid. For a long, uniformly wound solenoid 
of N turns and length L meters, carrying a 
current I amps, the field at the center of the 
solenoid is: 


NI 
Hy ctenoid = —— [a/m] (3.12) 


L 

The solenoid field, excluding the ends, is 

quite uniform. The end regions excluded are 

lengths of about a diameter. Such systems can 

be used with relatively large samples of mag- 

netic material which do not require very large 
magnetizing fields. 


3.5.8 THE ELECTROMAGNET 


For materials requiring very large magnetizing 
fields, electromagnets are commonly used. The 
electromagnet is a structure of (magnetically) 
soft iron, constructed so that there is a gap 
where a uniform field is produced. The gap has 
two cylindrical poles with flat faces separated 
by a distance of g meters. The structure is en- 
ergized by current J flowing in coils of N total 
turns. An electromagnet is sketched in figure 
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3.3. At distances of at least the gap-length (g) 
from the pole edges, the field is quite uniform: 


(3.13) 


NI 
A yap = a [a/m] 


SOFT \IRon Yo KE 
figure 3.3 


3.6 MKS UNITS 


The treatment of magnetic moments is differ- 
ent in two competing MKS unit system. These 
two systems are further described below. The 
first has magnetization (M) in the same units 
as H, as has been done here. The second uses 
a variable called Intensity of magnetization (I) 
and its units are the same as B. The defining 
equations are: 


Sommerfeld: B = p.(H + M) (3.14) 


(3.15) 


In linear magnetic materials, also called soft 
magnetic materials, the magnetization is pro- 
portional to H for fields less than some partic- 
ular level: 


Kennelly: B= p,-H+I1 


M = xH 
for example, and x (greek chi) is called the mag- 


netic suceptibility. For such materials, the per- 
meability is p: 


B= Hobr = po(1 + x) 


and p, is the relative permeability. 

The two competing MKS unit systems also 
differ in that the magnetic moment for a cur- 
rent loop of current ¢ and area A is defined dif- 
ferently for each: 


Sommerfeld: m = ?Al,, 
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Kennelly: j = potAl,, 


The definition of torque on a magnetic mo- 
ment is therefore different in the two systems. 
Sommerfeld has the torque proportional to the 
magnetic flux density B, whereas we have de- 
fined it as proportional to the force field pH 
which is in accordance with the Kennelly sys- 
tem. 


3.7 UNIT SYSTEMS 


A variety of unit systems have been developed, 
and several survive in use today. The System 
International, S.I., which is an MKS system, 
was devised by G. Giorgi in 1902. Electrical 
Engineering, has long been defined in this unit 
system. 

Two competing MKS systems are in use for 
magnetic units. They both use the same def- 
initions and units for magnetic field intensity, 
flux density and permeability: Magnetic field 
intensity: H [a/m], and Magnetic flux density: 
B [tesla]. And po = 4” x 10-7 [henry/m] 

The system most commonly used in text 
books on electromagnetic fields was invented by 
A. Sommerfeld, and was adopted by one of the 
most successful texts in the middle 20th cen- 
tury, Electromagnetic Theory by J.A. Stratton. 
Later authors tended to follow this text, so that 
the Sommerfeld system has become a de facto 
standard. The units used in the Sommerfeld 
System are: 


Magnetic moment: m [ampere-m?] 
Magnetization intensity: M [ampere/m] 


The defining equations are: 
B= pOo(H + M) 


T=mxB 


The competing MKS system is the Giorgi- 
Kennelly system, generally following Giorgi’s 
original MKS system, but was so strongly sup- 
ported by A.E. Kennelly that it is known as the 
Kennelly System. The units used in the Ken- 
nelly System are: 


Magnetic moment: j [weber-m] 
Magnetization intensity: I (weber/m7] 
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The defining equations are: 
B=p.H+I 
T =j x H 


Unfortunately there is an inconsistency in the 
two systems. The symbols j for magnetic mo- 
ment and I for magnetic intensity are employed 
in the Kennelly system and are used here to 
pursue the inconsistency as follows: 

As B and H are the same in both units, the 
first defining equations give that: 


po(H+M)=p4.H +1 
so that I= pp>M 


The magnetic moment of a volume is defined 
in both systems as the volume integral of the 
magnetization. Thus, with a particular vol- 
ume, assuming uniform magnetization within 
the volume, the magnetic moment in each sys- 
tem (the same volume, only the unit systems 
are different): : 


(3.16) 


M vol 
I vol 


For Sommerfeld: m = 
For Kennelly: } = 


In the presence of an applied field, the torque 
must be independent of the unit system: 


Sommerfeld: T = 
Kennelly: T 


m x B=(M x B) vol 
jx H=(1Ix H) vol 


Then, as B = pop,rH in both systems, we 
have that I= pop,yM 

This differs from the result from equation 
3.16: I = poM. These differ by the relative 
permeability of the media, and cannot both be 
correct. The problem is to determine where the 
error lies. 

In the paper, Experimental demonstration 
that the Couple on a Bar Magnet depends on H, 
not B, in Nature, vol 234, November 1971, pp 
31-33, R.W. Whitworth and H.V. Stopes-Roe 
of the University of Birmingham in England, 
reported the results of an experiment, which 
apparently proved the error was in the Som- 
merfeld system. This error is in the expres- 
sion for the torque on the magnetic moment, 
so that the magnetization intensity relation is 
I = poM, and the torque must be corrected in 
the Sommerfeld system: 
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Sommerfeld correction: T = m x poH 


The other surviving metric unit system is a 
cgs system called electromagnetic units (emu), 
which was one of the incomplete systems in- 
troduced by J.C. Maxwell in 1863. Yet another 
system which he introduced in 1863 was also an 
incomplete cgs system named the electrostatic 
unit system (esu) which is no longer used. The 
primary use of the emu system is in magnet- 
ics and magnetic materials, primarily because 
the International Electrotechnical Commission 
(IEC) has not yet set S.I. standards magnetic 
moment and magnetization. The emu system 
has the following magnetic units which are in 
common use: 


Magnetic field intensity: H [oersteds, oe] 
Magnetic flux density: B [gauss, g] 
Magnetic moment: m [ergs/oersted, or emu] 
Magnetization intensity: M [emu/cm? ] 
Specific magnetization: o [emu/gram] 


The defining equations are: 


B= 
T 


H + 47M 
mxH 


Thus it can be seen that the permeability of 
free space is unity in the emu system: p, = 1. 


3.8 MAGNETIC UNIT CONVER- 
SIONS 
Magnetic Field Intensity: H 
1 oersted =250/x amperes/m 


Magnetic Flux Density: B 
10,000 Gauss = 1 tesla 


Magnetic Intensity: I (Kennelly) 
4x x 104 emu/cm® = 1 tesla 


Magnetization: M (Sommerfeld) 
4m emu/cm® = (250/7) a/m 


Chapter 4 
LINEAR MAGNETIC 


4.1 LINEAR CHARACTERISTICS 
Magnetic materials exhibit a hysteresis loop in 
a form similar to that shown in figure 4.1. The 
usual parameters are the coercive force H,, the 
saturation flux density B, and the remnant 
magnetization B,. A fourth parameter, the sat- 
uration field H, is added for our convenience 
here. A so-called soft magnetic material has 
B, that is a small] fraction of B,. When its 
operation is limited to the unsaturated region 
({H| < H,) the characteristic can be modeled 
mathematically as a linear device with reason- 
able accuracy. In this operating range, ignoring 
the hysteresis for the moment, the relation be- 
tween the magnetic field intensity, H, and the 
magnetic flux density, B, is reasonably approx- 
imated by the linear relation: B = pH. (The 
hysteresis results in a power loss, which will be 
dealt with at a later point.) Referring to figure 
4.1, the permeability p is the slope of the unsat- 
urated region and is approximated by B,/He¢. 


figure 4.1 magnetic characteristic 


A more general way to express the relation 
of B and H is to include the magnetization M. 
Here the Sommerfeld unit system is employed 
because it avoids p, in relating M and H: B= 


MATERIALS 


15 


Ho(H +M). For a linear material, M = yH, 
where y is a constant. Because the divergence 
of B is zero, and 


B= ».(1+ x)H = pH 


the divergence of H is also zero inside of a linear 
material. In general 


din H + M) = 0 80 divH = —divM 


4.2 


Inside of a linear material the divergence H is 
zero and the divergence of M is zero. At the 
surface of the material, the interface between 
the magnetic material and a non-magnetic re- 
gion, there may well be a discontinuity of the 
magnetization. This occurs whenever there is a 
component of magnetization perpendicular to 
the surface. The discontinuity arises because 
the normal component of magnetization exists 
inside the magnetic material, but not in the ex- 
terior nonmagnetic region. In such a case there 
will also be a divergence of the magnetic field 
H at the surface. 

One viewpoint is to consider the continuity of 
the magnetic flux density, which cannot change 
as it crosses the surface separating the magnetic 
material from the nonmagnetic region. Desig- 
nating the normal components as Hy; and M,, 
just inside the magnetic material and Hp, just 
outside or in the nonmagnetic region, we have 
the relation that 


B, = Bo( Ani + M,) = Bono 


so that H,,. — Hni = Mn. 

The normal component of magnetization at 
the surface is a source of magnetic field inten- 
sity, and is experienced in the form of magnetic 


SURFACE POLARITY 
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poles. The magnetization is a volume density of 
magnetic moments, and its normal component 
at a surface results in a surface pole density. 
This is entirely analogous to a surface electric 
charge density on an electric conductor. The 
mathematical] treatment of Gauss’ Law for elec- 
tric charge density is applicable here. For an 
element of surface area dS, with a normal com- 
ponent of magnetization M,, the field element 
at a point p will be: 


M,dS 
dH — MP Gar, 


(4.1) 
where 1,, is a unit vector directed from the sur- 
face element toward the point p, and r,p is the 
corresponding distance. In order to determine 
the total field at the point p it will be neces- 
sary to perform an integration over the entire 
surface. 

As we are first interested in the effects quite 
close to the surface, it is reasonable to take the 
surface as an infinite plane. In performing such 
an integration it is perhaps easiest to set up an 
element of surface which is made up of a ring 
centered below the point p, from which only 
a normal component of field will result, and 
which, when integrated over the infinite plane, 
will yield a result that H,, = M,//2. 


4.3 DEMAGNETIZING FIELD 


Recalling that H,. — Hn; = Mn, it is seen 
that with H,, = M,,/2, it is necessary that 
Hn;y = —M,/2, so that the surface poles cause 
fields in both the nonmagnetic region and in the 
magnetic material. The interior field is in the 
opposite direction to the normal component of 
magnetization, and for that reason is called a 
demagnetizing field, Hg. 


NONMAGNETIC REGION 
MAGNETIC MATERIAL 


1z M 


figure 4.2 Demagnetization 


The magnetization 1s magnetic moment per 
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unit volume, and as such is 8 spatial average 
over 8 macroscopic point. Each individual mo- 
ment within a macroscopic point is subject to 
torque due to the local magnetic field, which 
act to align the magnetic moment with the field. 
Thus, the demagnetizing field caused by a mag- 
netization component normal to a surface tends 
to cause the normal component of magnetiza- 
tion to decrease. This can be seen in figure 
4.2 where the magnetization is shown at an an- 
gle a to the surface normal, together with the 
normal component and the resulting demagne- 
tizing field Hg. The torque on a magnetic mo- 
ment is the vectorcross-product, so the torque 
per unit volume on the magnetization is simi- 
lar: o.M x H. The cartesian coordinate sys- 
tem shown corresponds with: M = 1,M sina+ 
1,M cosa, and Hg = —1,(M/2)cosa, so the 
torque is: 


2 


[oM x Hz = 4 


. 2 


sin @ COS @ 


2 
—l1, ~ sin 2a 


It is seen that the torque is maximum at an 
angle of 45°, and has minima at zero and 90°. 
The minimum torque at a = 0 can be shown to 
be an unstable equilibrium point (as any dis- 
turbance in that angle will tend to increase the 
angle), and the minimum at a = 90° 1s a stable 
equilibrium for the magnetization. Thus, in the 
absence of any other torques on the magnetiza- 
tion, it would assume an orientation parallel to 
the surface. 


4.3.1 DEMAGNETIZATION 
APPLIED FIELD 


WITH AN 


For the case where there is an applied field act- 
ing on the magnetic material, we can break 
that field into a component tangent to the sur- 
face and a component normal to the surface: 
H, = 1;4;+1,Hn where H; = Hasin®@ and 
H,, = H,cos@ where 6 is the angle the applied 
field makes with the normal to the surface as 
indicated in figure 4.3. 

The total interior field is the vector sum of 
the applied field and the resulting demagnetiz- 
ing field: H = Ha + Hg. The magnetization 
will lie in the direction of H and will be propor- 
tional: M = yH. As was previously developed, 
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Ha 


NONMAGNETIC REGION 
MAGNETIC MATERIAL 


Int 
1; 
figure 4.3 Magnetization Near a Boundary 


the demagnetizing field is Hg = 


—$1nMp. 
Then the total field can be written: 


H = 1,7, + 1n(An _— 5Mn) 


As My = x(Hn — 3Mn), 


H=1,H,sin@+1,H,(1—- = ) cos 6 


x 


As M lies parallel to H, the angle the magne- 
tization makes with the surface normal is found: 


+X 
2 


For a susceptibility of y = 100, which is low for 
typical magnetic materials, for an applied field 
at 1°, the magnetization would be at an angle 
of 89 degrees to the normal. 

For the following cases of magnetic circuits 
without nonmagnetic gaps, it is a reasonable 
approximation to assume the magnetization 
near a surface will be parallel to the surface. 


2 
tana = tan 0 


4.3.2 DEMAGNETIZATION AT A GAP 


For magnetic circuits with nonmagnetic gaps it 
is necessary to reconsider the foregoing mate- 
rial. Here we consider a relatively thin nonmag- 
netic region sandwiched between two magnetic 
regions which have similar magnetic properties. 

The surfaces considered have large dimen- 
sions compared with thickness of the gap, and 
effects due to surface poles near the outer edges 
of the gap are considered far enough away to be 
negligible. The region being investigated is rep- 
resented in figure 4.4. The nonmagnetic region 
has a thickness g, and the magnetization within 
the two magnetic regions is assumed to be the 
same (or at least the components normal to the 
surface are both the same). 

As the magnetic material and the nonmag- 
netic material are both linear, the principle of 
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figure 4.4 Magnetization Near a Gap 


M, | M 


superposition is applicable. The surface be- 
tween the gap and the lower magnetic material 
has a positive surface pole density because the 
normal component of magnetization is directed 
toward the surface. This will result in normal 
field components pointing away from the sur- 
face, and having a magnitude of 5 My, indicated 
as H,., (pointing upward, extending through 
the gap and through the upper magnetic ma- 
teria] as well) and Hg, (pointing downward in 
the lower magnetic material). 

The surface between the gap and the upper 
magnetic material has a negative surface pole 
density because the normal component of mag- 
netization is directed away from the surface. 
This causes normal field components pointing 
toward that surface, and having a magnitude 
of 5Mn, indicated as H,., (pointing upward, 
but extending downward through the gap and 
into the lower magnetic material as well) and 
Hq, (pointing downward but extending upward 
through the upper magnetic material). 

Applying superposition, it is seen that these 
normal field components add within the non- 
magnetic material, and subtract within the 
magnetic materials, so there is no demagnetiz- 
ing field associated with the gap. 

The foregoing is an approximation, but is rea- 
sonably accurate in gap regions at least two 
gap-lengths (2g) away from any lateral edges 
of the gap. 


4.4 SCALAR MAGNETIC PO- 
TENTIAL 


For linear operation of soft magnetic materials 
with negligible hysteresis, the relation between 


LINEAR MAGNETIC MATERIALS 


the magnetic flux density, B and the magnetic 
field intensity, H, is linear, and the divergence 
of B is always zero, so that the divergence of H 
within a linear magnetic material is also zero. 
Any vector whose divergence is zero can be ex- 
pressed as the gradient of a scalar potential. 
Thus it is useful to express the magnetic field in- 
tensity within a linear magnetic material as the 
gradient of a magnetic scalar potential which 
we shall refer to as VW, which in the MKS unit 
system has the dimension of Amperes. Its gra- 
dient then has the dimensions of Amps per me- 
ter, and the relation is: 


H = —grad ¥ (4.2) 


As the divergence of H is zero in the linear mag- 
netic material, it is also true that divgrad ¥ = 
0, which is Laplace’s Equation. The operation 
div grad W is called the Laplacian of W, and is 
written in shorthand as V?W This result is that 
the Laplacian of the magnetic scalar potential 
is zero. In cartesian coordinates it is written: 
eu ay 


3p t Oy? +37 =0 (4.3) 


EXAMPLE 4.4.1. 


A block of linear magnetic material has a 
Square cross-sectional area and a length £, 
which is along the x-axis so that one end is 
at x = 0 and the other at z = @. Assuming 
the scalar magnetic potential is 100z within the 
block: 

(a) Determine the potentials at the two ends 
of the block. 

(b) Determine whether Laplace’s Equation is 
satisfied within the block, and if so 

(c) Determine the magnetic field intensity in- 
side the block. 


SOLUTION OF EXAMPLE 4.4.1. 

(a) With Y = 100z, the potential at z = 0 is 
zero, and atz = ¢@, ¥= 1008. 

(b) With W = 100z, the first partial deriva- 
tive with respect to x 1s a constant, so the sec- 
ond partial derivative with respect to x 1s zero. 
As no other derivatives exist, Laplace’s Equa- 
tion is satisfied. . 

(c) The gradient is: 


bv OW OW 
grad 1, Or yoy + LL 
= 1001, 
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Here is is seen that the field within the bar 
is H = —gradW = —1001,, which is perpendic- 
ular to the equipotential ends of the bar, and 
directed from the higher to the lower potential. 


END OF EXAMPLE 4.4.1. 


Laplace’s Equation in circular cylindrical co- 
ordinates is written: 


1ev &y 


EXAMPLE 4.4.2 

A uniform ring of magnetic material is cen- 
tered on the z-axis at the origin. It is assumed 
that the magnetic field is due to a current along 
the z-axis which will result in the field being cir- 
cumferential within the ring (and elsewhere as 
a matter of fact). Prove that the equipotential 
surfaces will be planes radiating from the the 
Z-axis. 


SOLUTION OF EXAMPLE 4.4.2 


The easiest approach is to assume the answer 
and demonstrate that the field is circumferen- 
tial. To do this, it is necessary to recognize 
that a plane radiating from the z-axis will in- 
clude a particular radius vector, p1,, at a par- 
ticular angle, ¢. It is assumed that the equipo- 
tential surfaces are a linear function of ¢ alone: 
W = Ko+Ai¢. The Laplacian of this potential 
is obviously zero, and we determine the field 
from its negative gradient. 

In the cylindrical coordinate system the gra- 
dient is: 

Oy 1. Ov OV 


gradW = ay + aay) + 1,5- (4.5) 


so the field is calculated: 


1 1 0( Ko + Ki¢) _ Ky 
o> oo — apa 


H = -— lg 


d¢ p 
and it is seen to be only circumferential, which 
agrees with the assumption. Again it is seen 
that the equipotentials are perpendicular to the 
field vectors, and the field is directed in the di- 
rection of decreasing potential. 


END OF EXAMPLE 4.4.2 


At any point in a linear material the magnetic 
field, H and the equipotential surface passing 
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through that point are perpendicular, as a re- 
sult of the relation H = —grad Y. This can be 
written on a small scale as HA£ = —AY, with 
the understanding that H and Aé@ are aligned, 
so that this is equivalent to a vector dot prod- 
uct. This relation is illustrated in figure 4.5. 


Aft Wi+ AY 


Wi 


figure 4.5 magnetic element 


The magnetic element in figure 4.5 is suff- 
ciently small that there is no significant vari- 
ation of H over the area A, and the lateral 
boundaries are parallel to H. As B = pH, 
the foregoing statements also apply to B. From 
there it is a simple step to find that the element 
flux, A®, which is B- A = ypHA enters the 
element through the surface marked A at po- 
tential YW; and leaves the element through the 
other end, which end is indicated by ¥, + AV. 
Then the relation between the potential differ- 
ence and the flux can be expressed as: 


Ag 


HAf=—Ab=-AY 4, 
TA (4.6) 


Then the relation of the flux, which enters 
the magnetic element at the higher equipoten- 
tial surface and exits at the lower equipotential 
surface, is geometrically related to the drop in 
equipotential. It is to be noted that the flux 
does not pass through the lateral surfaces be- 
cause those surfaces are defined as parallel to 


the field H. 


The relation of the element of flux to the 
drop in magnetic scalar potential is geometri- 
cal. This relation is often expressed in terms 
of the parameter reluctance (R) or its inverse, 
permeance (P). For the magnetic element this 
can be written as: 


(4.7) 
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4.5 MAGNETIC CIRCUITS 


The concept of the magnetic circuit is useful 
in dealing with high permeability linear mag- 
netic materials. It is based on Ampere’s Law of 
Circuital Magnetism, and upon the concept of 
the continuity of magnetic flux. Ampere‘s Law 
states that the magnetomotive force, F, around 
a closed path is equal to the current enclosed 
by the path. This is a definition based on the 
circulation integral of H: 


Fee L tenclosed by c =4 H - dl (4.8) 
¢ 

which is the counterpart of H = —grad V. The 
circulation integral can be broken into as many 
segments as necessary. For amenable geome- 
tries, it is possible to replace the integrations 
by a series of products of field and length on 
paths where the field is uniform: 


f H -dl= H,Al, + H2Al,+---+ H,Al, 
c 


It is convenient to conceive of each of the 
terms H;Al; drops of of magnetomotive force 
between two equipotential surfaces, that is, 
mathematical surfaces which each have a con- 
stant magnetic scalar potential. The length of 
path between the equipotential surfaces must 
be short enough that the field intensity, H, is 
reasonably constant along the path, and the 
path must be aligned with the direction of H. 
When these conditions are met, we can call the 
resulting segment of the magnetic circuit a leg. 

The continuity of magnetic flux is another 
way of saying the divergence of the magnetic 
flux density, B, is zero. This requires that 
what-ever flux enters a volume must also leave 
the volume: fluxoy; = flux;, 

Then the properties of a leg are as follows: 


e A leg is a length of magnetic material be 
tween two equipotential surfaces. 


e The magnetization within a leg is parallel 
to the lateral surfaces of the leg. 


e The magnetic flux which passes through 
one equipotential surface in the leg is the 
same as that passing through any other 
equipotential surface within the leg. Thus, 
the same flux which enters one end of the 
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leg exits at the other. That is to say that 
there is negligible leakage of flux from the 
sides of the leg. 


The leg of a magnetic circuit is similar to a 
magnetic element with the difference that the 
field H need not be uniform over the equipo- 
tential surfaces of the leg. When considering a 
segment of a magnetic structure as a candidate 
for a leg, it is to be remembered that the mag- 
netization prefers to align parallel to a lateral 
surface, which is no difficulty when the segment 
is straight. However, when a lateral surface is 
curved or contains a corner, it is necessary that 
the magnetization somehow conform to the sur- 
face. 

The way that the magnetization and mag- 
netic field conform to the lateral] surfaces is by 
the creation of just enough magnetic surface 
pole density to cause the internal field and the 
magnetization to nearly conform to the lateral 
surfaces. By not quite conforming to the sur- 
face, there will be a normal component of mag- 
netization, which gives rise to the surface pole 
density as discussed in section 4.3. Within the 
linear magnetic material M and H are colinear 
(i.e. they have identical orientations at any in- 
ternal point). For this reason, the assumption 
of no Jeakage flux from a magnetic leg is only an 
approximation. There must be sufficient leak- 
age flux to cause the magnetization to conform 
to the lateral surfaces. | 


figure 4.6 Toroid Core 


The simplest example of a leg is taken from 
the toroid having N turns of wire threaded 
through the aperture, so the magnetic field 
within the toroid is entirely in the circumfer- 
ential direction, and can be expressed as the 
circumferential component: 


(4.9) 


The toroid core has an inner radius of a and 
an outer radius of 6, with a thickness of A, as 
shown in fig. 4.6. 
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It is assumed that there is a total mmf of 
Ni amperes threading the aperture of the core, 
so the total potential drop around the toroid is 
also Ni. 

An arbitrary section of the toroid can be 
taken as a leg, as indicated in figure 4.7. The 
equipotential surfaces indicated as VW, and W, 
are radia] surfaces which would pass through 
the axis of symmetry of the toroid if extended 
that far. Using the cylindrical coordinate sys- 
tem, the potential on an arbitrary equipotential 
surface can be written as W(¢) = Ko + Kid, 
where ¢ is the coordinate angle, not to be con- 
fused with the flux ®. | 


VW; 


figure 4.7 magnetic circuit leg 


It can be shown that with potentials as given 
above, the magnetic field is perpendicular to the 
potentials and thus circumferential in direction, 
as is to be expected from Ampere’s Law of Cir- 
cuital Magnetism. 


Chapter 5 


RECORDING HEADS 


FIELDS 


5.1 INTRODUCTION 

In magnetic recording, a strong, localized field 
is needed to change the direction of magneti- 
zation in semi-hard (nonlinear) magnetic ma- 
terials. A piece of copper wire could provide 
enough field when soft materials such as iron 
wire were used for recording in the early times, 
and as Valdemar Poulson used in his original 
invention in 1898. 

A current carrying wire produces a field 
somewhat less than the long strarght wire field 
used to explain Ampere’s law. For the long 
straight wire carrying a current I, the field is 
the current divided by 4x divided by the dis- 
tance from the axis of the wire. American Wire 
Gage (AWG) copper size #40 has a radius of 
0.0508 mm, which means that the maximum 
field (at the surface of the wire) is about 3000 
A/m for 1 ampere of current flowing in the wire. 
The current capacity of this wire size is about 
20 mA, so the maximum practical field it could 
provide is about 60 A/m, or about 0.8 oersteds. 
AWG #30 copper wire has a radius of .152 mm 
and has a current capacity of about 100 mA, 
to provide a maximum field of about 100 A/m 
or 1.3 oe. at the surface of the wire. These 
are strong enough to magnetize iron, but iron 
is not a practical recording material, because it 
is magnetically soft and can easily be demag- 
netized. Practical recording materials require 
fields of from 300 to 1500 oe. which cannot be 
achieved with a simple copper wire. 

A magnetic circuit can be used to amplify the 
effective current by the number of turns, and 
can concentrate the flux in a very local region. 
However, to provide a recording field, the flux 
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must be released from the core and directed into 
the recording media. The simplest method of 
releasing the flux from its conduit is to put a 
nonmagnetic gap into the circuit, breaking it 
open and allowing the flux to spill out. 

Magnetic flux is continuous, and is contained 
by magnetic materials which behave as conduits 
of flux, bending the flux gently around curves 
in the surface. The mechanism for the shaping 
of the field to fit the shape of the magnetic ma- 
terial is by surface poles, which act as further 
sources of magnetic fields. Surface poles are 
caused by the abrupt change of magnetization 
at the interface between a magnetic material, 
which has magnetization, and the nonmagnetic 
material which does not. This is most readily 
seen by studying the divergence of the magnetic 
flux density, which is zero: 

B = p.(H + M) (5.1) 

As the divergence of B is zero, the divergence 
of H+ M must also be zero. This means that: 
div H = —div M. The divergence of magneti- 
zation is therefore a source of magnetic field. 

If there is a magnetic pole density pm, then 
we can use Gauss’ Law and Gauss’s Theorem to 
establish that with div H = py, the field due to 
such a pole density is related through Gauss’s 
Theorem as: 


/ : H-1,dS = ff i H)dVol (5.2) 


Which is Gauss’s Theorem, which can be ap- 
plied to any conservative vector field. This is to 
be understood as follows: The normal compo- 
nent of the magnetic field intensity is integrated 
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over the entire surface of a volume, and the net 
result is the same as if the divergence were in- 
tegrated over the volume bounded by that sur- 
face. Thus S,.) of eq. 5.2 is the surface of vol 
on the right side. Using the fact that the di- 
vergence of H is pm, eq. 5.2 can be written 


as: 
/ Hinds = ff / pmaVol (5.3) 
Svol vol 


Equation 5.3 can be further manipulated to 
give the algorithm for calculating an element of 
H at a field point due to an element of charge 
density at a source point. This is: 


| a Pmlsy 
dH, = tar? dVol 


(5.4) 


where r is the distance from the source point 
to the field point, and 1,, is the unit vector 
pointing from the source point toward the field 
point. The subscript of p is used on H, to indi- 
cate that this is a component of H that is due 
to magnetic pole density. This magnetic pole 
density is due to the divergence of magnetiza- 
tion: 

Pm = —divM (5.5) 


As was previously discussed, a linear magnetic 
materia] has M = xH, so that there is no di- 
vergence of M within linear materials, but can 
occur at the surfaces. This is the primary rea- 
son that magnetic circuits have the ability to 
contain flux. When the magnetization would 
otherwise have a component perpendicular to a 
surface, surface poles are created by the discon- 
tinuity of magnetization which are then sources 
of H, which create a torque on the magneti- 
zation which tends to make the magnetization 
parallel to the surface. 

The flux containment of magnetic circuits 
provides another benefit, as the fields due to 
the surface poles will compensate for lack of 
symmetry in the current configurations. 


5.2 RELUCTANCE OF 
NONMAGNETIC GAPS IN 
MAGNETIC CIRCUITS 


Gaps in magnetic circuits are of great tech- 
nological interest, because magnetic circuits are 
able to contain and concentrate magnetic flux. 
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In electrical machinery very strong fields are ob- 
tained in air gaps and are used to provide the 
emf for generators and torque for motors. 

The analysis of the fringe field around a gap 
will be undertaken later, we first must under- 
stand the effects of gaps on the magnetic cir- 
cuit itself. From extensive experience in the 
measurement of reluctance of air gaps in elec- 
trical machines, where gaps which are short 
when compared with the dimensions of the s0- 
called pole-faces, the reluctance is found from 
the length of the gap, g, and an adjusted area 
of the faces. The pole-faces are surfaces so de- 
signed that they would have been equipotentials 
had the gap not been introduced. 


fig. 5.1. Gap pole faces. 


For a rectangular cross-sectioned magnetic 
circuit as shown in figure 5.1, they have dimen- 
sions of width w and throat height h, so the 
pole-face area is hw. The gap length (between 
the pole faces) is g. The adjustment neces- 
sary to reasonably approximate the reluctance 
is to increas the effective area by adding the gap 
length to both w and h before calculating the 
effective area of the pole faces. Thus, following 
from equation 4.7 in the previous chapter, the 
reluctance of this gap is: 


secant , 

Rew = Tht owra) 
In ferrite recording heads, to obtain very flat 
pole faces, the head is made in two parts, to 
allow the two pole faces can be machined and 
polished. These two halves of the head are then 
put together to create a magnetic circuit having 
two gaps, one called the front gap, and the other 
the back gap. A simpified head geometry is 
shown in figure 5.2 
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Lent gap 
flux is useful in certain geometries where the 
thickness is constant. This makes use of the 
reluctance relation given in eq. 4.7, which is 
| | _ repeated here for convenience: 


| | | For a rectangular leg, the cross-sectional area 
is the product of its width and height, A = 
wxh. 


ay LR i es 
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fig. 5.2. simplified head geometry. | desig es 
a VBES 


For this particular example, for purposes of 
illustration, the height of the back gap (h,) is 
3.25 times the height of the front gap (hy). The 
reluctance of the front gap is: 


Ry, = ae ee nc 
4 powhy(1+9/hy)(1 + g/w) 


and for the back gap: 


Te 
powhe(1 + g/hs)(1 + g/w) 


The total reluctance of the magnetic cir- 
cuit including the nonmagnetic gaps (typically 
called air-gaps, although the material is not 
usually air) is: 


Riotal = Reore + Ryy = Rog (5.7) 


Then it is necessary to determine the reluc- 


tance of the magnetic material portion of the 


circuit, referred to as the core. 


5.2.1 MAGNETIC CORE RELUCTANCE 


To deal with the magnetic material in the 
magnetic circuit of figure 5.2, we eliminate the 
air gap for the moment, to generate a form as 
shown in figure 5.3. The gaps are in series with 
this part, and so will be added in later. 

A non-mathematical technique for the deter- 
mination of the relation between the MMF and 


Cut 


fig. 5.3. circuit with gap removed. 


The magnetic circuit can be understood to 
be composed of tubes of flux, which close upon 
themselves in enclosing the MMF. The poten- 
tial drop around each flux-tube circuit is equal 
to the MMF. The element of figure 5.3 is part 
of one tube of flux, indicated as tube A. No flux 
leaves a tube, though the tube may change in 
cross-sectional area in making the circuit. The 
quantity which is the same throughout the cir- 
cuit of a tube is the element’s flux Ad. The 
potential drop from one end of the element to 
the other is: 


Al 


A® (5.8) 


where the negative sign indicates that the po 
tential is actually dropping in the direction of 
the flux. 


The magnetic circuit geometry is subdivided 
into equipotential surfaces and a number of flux 
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tubes which are perpendicular to the equipo 
tential surfaces. The thickness of the circuit is 
uniform at w, and elements of a tube are con- 
structed so that the length to width ratio (AZ 
to Ah) is 1.0 on all elements. In doing this, the 
reluctance of each element is the same as any 
other. Thus the total] reluctance of the tube is 
simply the number of cells times the reluctance 
per element: 


] 
Rude —le- 
pw 


If the flux tubes are chosen so that they each 
contain the same flux, then the total flux will 
be the number of tubes (n;) multiplied by the 
flux in one tube. As reluctance is mmf divided 
by flux, the total reluctance will then be the 
reluctance of one tube divided by the number 
of tubes: 


(5.9) 


In figure 5.3 the inner and outer surfaces are 
circles, and they are the inner and outer sur- 
faces of the total flux tube. We can divide 


this into smaller tubes by using other circu- 


lar shapes as indicated on the right hand side, 
where the inner tube was laid down first by 
drawing a circle through a point on the top cut 
a quarter of the way from the inner to the outer 
surface, and a point on the bottom cut a quarter 
of the way from the inner to the outer surface. 

Next, the cells were drawn in on the inner flux 
tube, making them as nearly square as possible, 
keeping the rule that flux tubes and equipo 
tential surfaces must be perpendicualr to one 
another. 

The next flux tube was constructed by ex- 
tending the equipotentials of the first tube, and 
then fitting the best circle to cause those ex- 
tensions and the circle to make the most nearly 
square cells possible. In doing this we assure 
that the flux tubes contain the same flux. 

The third tube was constructed using the 
same technique as the second. 

A full fourth tube would not fit in, so the 
equipotentials were extended to the outer sur- 
face. Then one of the oblong cells was divided 
up into the best squares to determine how wide 
it was. Because there are four squares by two 
squares in this oblong, the width is one half the 
length, so the outer tube contains only half the 
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flux of the others. This gives a total number of 
tubes as n; = 3.5. 

In counting the number of cells in a tube, on 
the half drawn there are 15 and a fraction. To 
determine that fraction, a convenient shortened 
cell is divided into the best squares and we find 
it is 4 wide by 2 long, so its length is a half a 
cell. Thus there are 15.5 cells in half the tube 
circuit for a total of 31 cells in a tube of flux. 

The total reluctance of the magnetic material 
portion of the circuit is then: | 


31 1 
Reireuit = =~ — 
3.5 pw 


A representative recording head has hy, = 
w © 259, (with hy, = 3.25h;,). Taking these 
values in this example, we have that: 


Riot 


= = ee ae eS 
wow \3.5p,  1.025hy,  1.025h,, 


: (= + 0.0474 
BoW \ Hr 


where p, is the relative permeability of 
the core material. Practical head materials 
have relative permeabilies from about 1000 to 
10,000. At the low end, we can see from eq. 5.9, 
the core will contribute about 16% of the total 
reluctance, and at the high end it would con- 
tribute about 2%. With a typical relative per- 
meability of 5000, the core contributes about 
4% of the total reluctance. The use of more 
massive cores can reduce this to 1 or 2% of the 
total reluctance. 


5.3 HEAD EFFICIENCY 


The purpose of the recording head is to pro- 
vide a strong local magnetic field. This field is 
the fringe field of the front gap, which we would 
like to have as large as possible. The potential 
drop across the front gap divided by the gap 
length will be the field at the geometrical center 
of the pole faces, and this will be the same over 
most of the gap, beginning to drop off at about 
one gap-length distance from the pole edges. 

The total potential drop around the magnetic 
circuit is the MMF due to currents threading 
the aperture. The potential drop of each part 
of the circuit will be the flux multiplied by the 
reluctance of that part. In order to maximize 
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the potential drop of the front gap, we would 
like to make the reluctance of the front gap large 
compared with the other reluctances in the cir- 
cuit. Eq. 5.7 shows the reluctances of the core, 
the front gap and the back gap. The relative 
permeability and size of the core determines its 
reluctance. The throat height of the back gap 
is the only factor which makes it different from 
the front gap in the simple head of figure 5.3. 
To minimize the effect of the back gap, it is nec- 
essary to make hp, much larger than hy,. Typi- 
cally the ratio of hy, /hy, is more than 10:1, and 
as much as 20:1. Thus, to make a very efficient 
head, the back gap throad height hj, must be 
much greater then hy,, and the core material 
must have a very high relative permeability. 

The head efficiency n is the ratio of the front 
gap reluctance to the total reluctance of the 
magnetic circuit: 


100Ry, 


a a RE ERR 5.10 
Ne eR Rea et) 


The deep gap field in the front gap is the po 
tential drop across the gap divided by the gap 
length: H, = Wyap/g, and the potential drop 
across the front gap is the MMF multiplied by 
the fractional head efficiency: 


Wrap = 9H, = —- MMF 


0 (5.11) 


5.4 KARLQVIST’S FRINGE 
FIELD 


In our example with a relative permeability 
of 5000, we found the core had only 4% of the 
potential drop. With higher permeabilities the 
potential drop around the core becomes even 
less. It is a short step to assume infinite per- 
meability of the core, and thereby make the 
core parts of the magnetic circuit equipoten- 
tial volumes. In practice the potential along 
the core changes very slowly as compared with 
the rate of change across the gaps. In addi- 
tion, the length of the gaps is small compared 
with the length of the pole tips for ferrite heads. 
The model used by O. Karlqvist was infinitely 
long pole pieces separated by g, each with infi- 
nite width wand infinite throat height hy. The 
geometry of his model is shown in figure 5.4. 
The left pole is assumed to be at a potential of 
gH,/2, and the right pole at -gH,/2, making 
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the midplane of the gap to be the zero equipo- 
tential, and is taken as the plane z = 0(with the 
positive direction to the right in figure 5.4). The 
outer surface of the core is taken to be the plane 
y = 0, ignoring any curvature, and the positive 
y-direction is outward from the core (up in fig- 
ure 5.4). The other coordinate is z, which has 
its positive direction out of the plane of the pa- 
per. The plane z = 0 is half-way across the 
head width (w). We shall carry out the analy- 
sis on the plane z = 0, assuming no variations 
in the z- direction. The problem then reduces 
to one of 2-dimensional! potential theory. 


U=gE,/ 2 


fig. 5.4. Karlqvist head geometry. 


In 2-dimensional potential theory, the poten- 
tial in the y > 0 half-plane can be calculated 
if the potential on the y = 0 axis (z-axis) is 
known. Except for along the top of the gap, 
called the gap region, the potential is known. 
Deep within the gap, the potential can be found 
to be W = —H,z’. For this analysis we assume 
that same potential distribution holds along 
y = 0, in the gap region, from z’ = —g/2 to 
z' = g/2. The potential for y > 0 is computed 
from: 


"poo ! 
; ydz 


1 
V(z, y) oe x [ 

(5.12) 
where (z,y) is the field point, and (z’,0) is 
the source point and the integration is over the 
source point on z’ from negative infinity to pos- 
itive infinity. The integration is performed and 


=~ OO 
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the result is: 


(z — £) arctan =—4 
— (z+ 4) arctan =+4 


+4)n ts 503 
(5.13) 


The field components, H, and Hy, are found 
by taking the negative gradient of W, which is 
a rather tedious process, but which eventually 
results in the following: 


KARLQVIST FRINGE FIELD EQUATIONS 


fae ee - (5.14) 
— arctan : 


H,, (z-4)?4+y 
an (z+ fh)? +y? 
In these equations, lengths are commonly nor- 
malized to half the gap length g/2, and upper 
case variables are used for normalized lengths: 
X = 22/g, and Y = 2y/g. 

The field components are normalized to H,, 
and lower case h’s are used: hy = H,/H, and 
hy = H,/H,. Then these equations in normal- 
ized form are: - 


W(z,y)= = 


Ay= (5.15) 


NORMALIZED KARLQVIST EQUATIONS © 


hx =< (arctan 2 eran 
(5.16) 
— 1, (X-1)4+Y? 
hy = on IDX ips (5.17) 


These expressions give components which have 
constant values which are circles in the z-plane, 
(actually cylinders parallel to the z-axis in three 
dimensions). 

For the perpendicular component hy = con- 
stant, the cylinders are centered on the y = 0 
plane, and have field directions upward for neg- 
ative values of z and downward for positive val- 
ues of z as indicated in figure 5.5. Cylinders 
with equal radii have the same magnitude of 
constant field. 


For the perpendicular component hy = con- 
stant: 
cylinder axis is lime: X = —cothahy, Y = 0 
cylinder radius of R? = esch*zhy 


The z-component is called the longitudinal 
component because the relative motion between 
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fig. 5.5 Karlqvist cylinders of Hy. 


the head and medium is in the z-direction. The 
cylinders for the longitudinal component hy = 
constant are shown in figure 5.6. These cylin- 
ders are tangent to the corners of the pole- 
pleces. 


fig. 5.6. Karlqvist cylinders of H,. 


For the longitudinal component hy = constant: 
cylinder axis is line: X = 0, Y = cot thx 
cylinder radius of R? = csc? rhx 


We are primarily interested in the longitu- 
dinal component because most recording me 
dia is essentially longitudinal, i.e. where the 
recorded magnetic moments are essentially in 
the x-direction. There may yet be interest in 
perpendicular media, where the magnetization 
is in the y-direction. 

The Karlqvist Head function is an approxi- 
mation, but has been found to be quite ade 
quate for many purposes. In some cases a s1m- 
pler model is useful in understanding the com- 
plex process of magnetic recording. 


5.5 FAR-FIELD MODEL 


When the recording medium is at least one 
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gap length away from the gap, the gap-length 
can be reasonably ignored in calculating the 
field. The Karlqvist’s result is not easily modi- 
fied to obtain this result, so one must begin with 
eq. 5.12 and set the z-axis potentials to 9H, /2 
from negative infinity to z = 0, and —gH,/2 
from z = 0 to positive infinity. This eliminates 
one integration and yields: 


H 
Ww = 22 arctan Z 
1 z 


Then, taking negative partial derivatives to 
obtain the field components we obtain that: 


FAR FIELD EQUATIONS 


25 gH, y 
_ _9Hy 2 | 
H, = ape ee (5.19) 


The locus of points where H, is constant, from 
eq. 5.18 is where 


2,,2_ 99 = 
a TH, ° 


As the equation of a circle is (zr — z-)? + (y— 
y-)? =r, where the point (z¢, yc) is the center 
of the circle, and r is the radius, equation 5.18 
can be put into the same form by adding the 
term (gH,/27H;)* to both sides (completing 
the square). Thus the locus of points where 
H, is of constant magnitude is a cylinder with 
its axis through the point (0,gHg9/27H~) and 
parallel to the z-axis and having a radius of 
(g|H,|/2x|Hz|). This is shown in figure 5.7. 


‘P 


fig. 5.7 Far field H; cylinders. 
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Similarly, equation 5.19 can be put into the 
circular form for constant Hy with a cylin- 
drica] locus with the axis through the point 
(—9H,/2xH, ,0) and parallel to the z-axis with 
a radius of (9|H,|/27|Hy,|), as indicated in fig- 
ure 5.8. 


fig. 5.8. Far field H, cylinders. 


Chapter 6 


_ THE RECORDING PROCESS 


The recording process is complex, involving 
head fields, media characteristics and the de- 
magnetizing fields which arise because of the 
magnetic record. The demagnetizing effects al- 


ter the recorded pattern, complicating the pro-. 


cess to a point where it is nearly incomprehen- 
sible to the student. To simplify the procedure 
of learning about recording, the demagnetizing 
field effects are postponed to the next chapter. 
Thus the development in this chapter results in 
an idealization of recording which is to be taken 
as the limiting case. Actual recording can only 
be worse than what we find in this chapter. 


6.1 CHARACTERISTICS OF 
RECORDING MEDIA. 


The magnetic materials employed for record- 
ing are generally classified as semi-hard mag- 
netic materials. The linear materials we have 
discussed for magnetic circuits are soft mag- 
netic materials, in that they are well character- 
ized by permeability. Permanent magnets are 
hard materials, and are characterized by an en- 
ergy product, which will not be covered in this 
course. 

Recording materials are classified primarily 
by their coercive force, H., and a squareness 
ratio. Figure 6.1 shows a typical saturation 
M —H characteristic of this type of material, 
with the coercive force H,, saturation magne- 
tization M, and saturation remnant magneti- 
zation M,, indicated. The squareness ratio is 
M,,/M,. For this characteristic can be either 
static, where the field is increased an increment 
and the magnetization measured at that field in 
e vibrating sample magnetometer (VSM), and 
then the field incremented and the measure- 
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ment repeated, or in a dynamic state, using a 
of a B — H looper under the assumption that 
M = B/po. 


Figure 6.1. Saturation hysteresis loop 
for a magnetic recording medium. 


The saturation M —H characteristic, or M— 
H loop as it is usually called, is of some value 
for recording, however it only approximates the 
information needed to study the recording pro- 
cess. In the recording process the recording 
medium moves past the head, experiencing a 
field while it is near the head, and then passes 
beyond the head field to a region of only am- 
bient fields. Thus it receives a field and then 
relaxes to a zero field condition with what- 
ever remnant magnetization it received from 
the recording field. Thus a characteristic for 
recording materials in which the remnant state 
which occurs after a field has been applied and 
then removed is of greater interest. 

Data and video recording use a technique 
which is called non-return-to-zero (NRZ) or sat- 
uration recording. This is largely used in instru- 
mentation recording and in digital audio record- 
ing as well. In this technique, information 1s 
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stored in the form of reversals of the magneti- 
zation from one direction to the other, t.e. from 
the positive longitudinal direction to the nega- 
tive longitudinal direction (the z-direction with 
respect to the head coordinates in figure 5.4). 
These reversals of magnetization are commonly 
called flux changes or flux reversals, and infor- 
mation is coded using these flux changes. In 
this process recording is accomplished by driv- 
ing the head winding current first in one direc- 
tion and then reversing its direction. For each 
direction the current must be sufficient to satu- 
rate the local recording medium. So recording 
takes place by reversing the direction of a sat- 
urating field so as to reverse the direction of 
magnetization, so that the recording medium 
will experience a full saturating field except 
during the short time interval required to re- 
verse the current direction. For this reason, 
the sort of an M — H characteristic needed is 
one where the material begins at one satura- 
tion level and is driven with a certain field in 
the reversed direction. the characterstic of in- 
terest in recording is then a quasi-static one, 
which is shown on the following page in figure 
6.3, and is called the remnant magnetization 
characteristic. This characteristic is not a loop 
as would be obtained with the dynamic M —H 
characteristic, but is a two branched curve, one 
branch at negative remnance curving upward to 
the point where the field just gives full remnant 
magneitzation. The other begins at positve 
remnance and curves downward to the nega- 
tive field that just reaches negative remnance. 
These two branches do not intersect. 


The remnant magnetization characteristic is 
obtained by a series of measurements, each be- 
ginning by driving the magnetization firmly to 
negative saturation. Then a positive field is ap- 
plied for a short period of time and then re- 
duced to zero, so that the magnetization re- 
turns to an intermediate remnant state. The 
measurement system required is a vibrating 
sample magnetometer (VSM). The applied field 
is measured when applied, and the remnant 
magnetization is then measured after the field 
is reduced to zero. This process is indicated in 
figure 6.2 where the dynamic saturation M—H 
characteristic of figure 6.1 is shown as a dashed 
curve. The remnant magnetization is plotted 
against the field for a series of these tests, us- 
ing increasing values of field for each succeeding 
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test, until the maximum remnant magnetiza- 
tion M,, 1s obtained. 7 


Figure 6.2. Generation of the remnant 
magnetization characteristic 


The parameters of interest to us from the 
remnant magnetization characteristic are the 
the saturation remnant magnetization M,,, the 
threshold field Hz and the switching or satura- 
tion field Hs, as indicated in figure 6.3. The 
threshold field Hy is the amount of field nec- 
essary to change magnetization from —M,, to 
—.9M,,, and the saturation field Hg is the field 
necessary to change the magnetization from 
—M,, to +.9M,,. 


An ideal recording medium would have Hs = 
Hr 


Figure 6.3. Remnant magnetization 
characteristic 
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6.2 LONGITUDINAL 
RECORDING. 


At the present time, most magnetic record- 
ing is tn-plane, that is the recording medium 
is in the form of tape or a disk which has a 
thin layer of recording medium coated or de- 
posited on the surface, and the recording takes 
place essentially with the magnetization lying 
in the plane of the medium. For our purposes 
the principal component of field which affects 
in-plane magnetization is the longitudinal com- 
ponent of the head field (z-component in equa- 
tions 5.14 and 5.16), so we shall concentrate on 
the action there. While there are experimental 
systems which use materials which magnetize 
only in directions normal to the plane, and as 
such utilize the perpendicular (y—direction in 
equation 4.15) field, we shall not consider per- 
pendicular recording at this time. 

Figure 6.4 is a schematic representation of 
the head-medium interface, with parameters of 
head-medium separation d, media thickness 6 
and gap length g. The recording media is at- 
tached to the substrate which is moving to the 
right past the gap with a velocity v, which a 
constant. The head has a width w perpendic- 
ular to the paper, and it is assumed there is 
no field variation across the width. For most 
purposes we shall] use normalized values for the 
parameters d, 6 and v, normalized to the half 
gap length g/2. 


Le x ane th ° / of Faane 
a ae Dae ae 
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Figure 6.4. Head-medium interface 


NORMALIZED PARAMETERS 
head-media separation: D = 2d/g 
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medium thickness: A = 26/9 
velocity: V = 2v/g 


For saturation recording, the applied field 
must be sufficient to saturate the medium di- 
rectly above the gap. We consider the cylinder 
of constant H, (as discussed in chapter 5) which 
passes through that point, and adjust the gap 
field H, so that the constant value of H, on 
that locus is Hs. This insures that all points 
closer to the gap are saturated, as all points 
within that cylinder experience a field greater 
than Hs, neglecting the demagnetizing field, as 
illustrated in figure 5.6. The highest reach of 
a cylinder is at x=0, and the height that the 
Hz = Hs cylinder must reach is d+ 6. Then 
setting z=0,y=d+6,and H, = Hs, we can 


_ solve for the minimum value of H, to insure 


saturation recording. In this case we wish to 
use regular rather than normalized units, as we 
are setting H,, the field used for normalization. 


Applying equation 4.14a, we obtain: 


THs 


oie 2 arctan[g/2(d + 6)] oe) 


It is seen that the distance parameters in 
equation 6.1 become normalized in the process, 
so thatnormalized parameters can almost al- 
ways be used. In general the field at the back 
of the medium (X =0, Y = D+ A) should ex- 
cede Hs, so that equation 6.1 is the minimum 
for saturation recording, and 


aHs 


2 —— ; 
Hy 2 2 arctan 1/(D+ A) (6.2) 


All the magnetization within the constant H; 
cylinder of H, = Hs is guaranteed to switch to 
the direction of the gap field. That is because 
all of the interior of a cylinder of constant A; 
has H, greater than that on the cylinder. Thus, 
when equation 6.2 is satisfied it is known that 
all media magnetization within the cylinder is 
in the direction of the gap field. 

Let us assume that the gap field satisfies 
equation 6.2, and that the field has been on for 
some time. The medium is moving to the right 
with a velocity v, so that for some distance to 
the right of the gap the media is saturated in 
the direction of H,. In figure 6.5 the situation 
is shown where the gap field is in the negative 
z-direction an has been in that direction for a 
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long enough time so that the the magnetiza- 
tion to the right of the cylinder marked —Hss, 
indicating that it is the cylinder of H, = —Hss, 
is switched to the negative direction as indi- 
cated. Further to the left we see a cylinder 
marked —H7, indicating it to be the cylinder 
for H, = —Hr, which is just enough to begin 
switching some of the magnetization to the neg- 
ative z-direction, if the magnetization at that 
cylinder were known to be saturated in the pos- 
itive z-direction. The region between these two 
cylinders would be capable of being partially 
switched to the negative z-direction, however 
we know nothing of the previous history there, 
and so will leave everything to the left of the 
H, = Hs with a question mark. 


Figure 6.5. Condition before 
head current reversal 


With the condition of figure 6.5 stabilized, 
the gap field is instantaneously reversed by re- 
versing the head current. The new cylinders in 
figure 6.6 are seen to overlap the previously sat- 
urated region in the medium near the recording 
gap. Here the previous condition of the media 
to the right of the gap is well known, it was 
previously saturated to the left. Thus we can 
predict that the region to the right of the gap 
between the cylinders marked Hs and H7 are 
partially reversed, so that region can be con- 
sidered to be a transition region between the 
earlier situation and the present condition. 


It can be seen that the transition length is 
determined by the two material parameters Hs 
and Hr, together with the normalized (to the 
half gap-length) parameters D and A. The 
cylinder radii of the the loci of H, = Hs and 
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Figure 6.6. Condition at the time 
of head current reversal 


H, = Hr have the ratio (see equation 5.14): 


(Rr _ sinths , Hs 
Rs  sinzhr Hr 


(6.3) 


The approximation is within 5% when hs is 
less than 0.55. 


6.2.1 Example 


For a particular recording system, the normail- 
ized parameters are D = 1 and A = 5. The 
recording media has Hy = 280 oe. and Hs = 
580 oe. We will set the gap field according to 
equation 6.2, and then determine the transi- 
tion length at both the top and bottom of the 
medium. From equation 6.2 we obtain that H, 
must be 5,517 oe., then hs = 0.330 and hz = 
0.159 Then from equation 6.3, Rr/Rs = 2.059, 
while Hs/Hr = 2.071. 

In order to solve for X from equation 4.15a, 
note that it can be manipulated to the form: 


] 2Y 
hz = = arctan X?4Y2-1 (6.4) 
so that 
X?=1-Y?42Ycotxh, (6.5) 


Then at Y = D=1, Xs = 2.416, and X7 = 
3.532, so the normalized transition width at the 
bottom of the medium is 1.115, which corre 
sponds to 0.558g. 

AtyYy = D+A =6, Xs = 0, and X7 = 6.311, 
which is also the transition length at the top of 
the medium, 3.165 g. 
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6.2.2 Example 


Using the same system as example 6.2.1, but 
increasing H, by 20% to 6,620 oe. hs becomes 
0.275 and hy becomes 0.133. Then at Y = D 
we obtain Xs = 2.662 and X7 = 3.866 for a 
normalized transition length of 1.204 or 0.602g. 
AtY = D+A=6, Xs = 2.744 and X7 = 
7.395 for a normalized transition length at the 
top of the medium of 4.651 or 2.326¢. 


end of example 


Figure 6.7. Transition regions for 
examples: (a) 6.2.1, (b) 6.2.2 


The transition regions for these two exam- 
ples are shown in figure 6.7 as the shaded ar- 
eas. While they are quite similar in shapes, 
the second can be argued to be better, in that 
it is more concentrated than the first. On the 
other hand, it may be argued that first is su- 
perior because the transition length nearest the 
head is shorter. As other considerations, such 
as demagnetization and the reproduction pro- 
cess, are of importance, it is not completely ac- 
curate to base an optimum on either of these 
criteria. However, experiments find an opti- 
mum recording current in the range of these 
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two values, t.e. between placing the H, = Hs 
cylinder so that it reaches beyond the back of 
the medium, but with the axis of the cylinder 
not beyond the middle of the medium. 

From the exercises 6.2.1 and 6.2.2 it can be 
seen that the ratio of Hs to Hr is of great im- 
portance in determining the transition length. 
To illustrate this we shall consider the same ge- 
ometry with a different Hr. 


6.2.3 Example 
The medium has Hs = 580 and Hr = 350 oe. 


_ Here we set the center of the H, = Hs cylinder 


axis at Y = D+ A/2, i.e. the middle of the 
medium, requiring that H, be 6,547 oe (very 
similar to example 6.2.2) and find hs = 0.278 
and hy = 0.168. Then at Y = D we have X5 = 
2.647 and X7 = 3.434, for a normalized transi- 
tion length of 0.789, or 0.395g. 

At Y = D+A we find Xs5 = 2.655 and X7 = 
5.980 for a normalized transition length of 3.325 
or 1.662¢. 


end of example 


Comparing the results of examples 6.2.2 and 
6.2.3, we find with a 20% reduction of Hs/Hr 
the bottom transition width is decreased 34% 
and the top by 29%. 

It should be clear that the ideal recording 
medium would have Hy very nearly Hs. 


Figure 6.8. Example 6.2.3 


6.3. DYNAMIC SWITCHING 


In the previous section we considered tran- 
sitions for the case where the head current re- 
versed instantaneously. That could also be a 
reasonable approximation whenever the media 
moves only slightly during the time that it takes 
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for the head field to be reversed. In this section 
we shall examine the situation where the me- 
dia moves appreciably during the time the head 
field is reversing. 

Here we shall make use of the fact that the 
gap field changes from one saturation direc- 
tion to the other in some continious way, nec- 
essarily because of the energy stored in the 
magnetic field, which cannot change instanta- 
neously. This can be thought of as the effect of 
the inductance of the head winding. Referring 
to figure 6.9, we will examine the field reversal 
process from the medium’s point of view. The 
medium sees the head traveling past to the left, 
carrying it’s field with it. We begin our exami- 
nation when the head just begins to switch it’s 
field from the negative z-direction to positive. 
The right-hand circle in figure 6.9(a) indicates 
the constant H, = —Hs cylinder at that time. 
As time goes on the head moves to the left while 
the field is first reduced from its original nega- 
tive z- direction magnitude to zero, and then re- 
versed until reaching the same magni(Etude but 
in the positive z-direction. The center point in 
figure 6.9(a) is the point where the field has 
reached zero. The left-hand circle indicates the 
fully reversed H, = +H cylinder. 
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(dD) magnetization pattern recorded 
Figure 6.9 


Figure 6.9(a) indicates the size of the H, = 
Hs cylinder for a discrete number of time sam- 
ples, as to show more than that would further 
obscure the point at hand. Figure 6.9(b) shows 
the resulting recorded pattern of magnetization 
after the head passed out of the picture. Here 
again we see the pattern for the discrete sam- 
ples of figure 6.9(a). Each portion of the bound- 
aries of this discrete pattern is an arc of a circle 
corresponding to a segment of the H, = Hs 
cylinder for a particular gap field and time in- 
stant. However, with the field changing in a 
smooth manner, any locus of constant H, is 
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also changing continuously. Thus there would 


- be a continuum of circles in figure 6.9(a), re- 


sulting in a smooth curves for the boundaries 
of figure 6.9(b). 

The medium is moving in the +z-direction 
with constant velocity v. Viewed from the 
medium reference frame, the head is moving to 
the left with that velocity. When the switching 
begins, the region of the medium which is ini- 
tially in the H, = —Hg cylinder is magnetized 
in the original direction, and the region imme 
diately to the left has been partially switched 
to to that direction. As the field falls from the 
Original direction and the head moves to the 
left, any point that has been saturated in the 
negative z-direction will remain in that state 
unless and until it experiences a reversed field 
of +Hr. In figure 6.9 we can see that the 
middle region of the medium, marked with a 
question mark, did not receive a field of Hs in 
either direction during this switching process. 
The previous recording history is unknown, so 
there is no way of determining the magnetiza- 
tion in that region, Data recording writes new 
data over the old record, so that the presence of 
uncontrolled regions such as the center region 
of figure 6.9 must be avoided. If they are al- 
lowed to occur, they will bring excessive noise 
to the play-back signal. 

To eliminate noisy regions such as the cen- 
ter parts of figure 6.9, we must insure that the 
pattern generated by the growing H, = Hs 
cylinder overlaps the pattern of the previously 
shrunk H, = —Hs pattern. That is to say, 
in the process of flux reversal the new direc- 
tion of magnetization must overlap and reverse 
the direction of part of the magnetization just 
previously recorded. (This is sometimes called 
over-write, but that term has another mean- 
ing which causes some confusion.) During this 
field reversal process it is necessary to assure 
that each point experience a field of magnitude 
Hs at least once. Then those points which 
originally had experienced a —Hs field will be 
switched to some extent when they experience 
& maximum reversing field greater than +H7, 
and will be fully switched if they experience a 
maximum reversing field of +Hs5 or more. This 
will insure uniform transition patterns and sub- 
sequently uniform play-back signal waveforms. 
A particular point in this active region which 
has previously experienced a field of —Hs and 
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then experiences a reversing field of at least H7, 
its final remnant state (demagnetization being 
ignored) is determed by the remnant magneti- 
zation characteristic and the peak of reversed 
field experienced. 


Here we examine a point which is experienc- 
ing its highest reversing field. It will next move 
to a lower field and be finished with the re- 
versal process. Referring again to figure 6.9, 
as the H, = —Hs cylinder shrinks it leaves a 
boundary where the magnetization to its night 
is saturated in the negative z-direction, and the 
magnetization to its left is undefined. Similarly 
as the H, = +Hsg cylinder grows it also es- 
tablishes a boundary where the magnetization 
is saturated to its left and undefined to the 
right. We wish to remove the undefined region 
by overlapping the H, = +Hs boundary oblit- 
erating the H, = —Hs boundary, to establish 
a predictable transition region. 


The boundary is a sequence of points which 
have experienced a field of Hs as their peak re- 
versing field. In two dimensions the cylinders 
appear as circles, and each point on a bound- 
ary obtained from a unique circle. Two nearby 
boundary points have come from two circles, 
one occurring after the other. If we place these 
two points very near to each other, it should be 
clear that a line drawn between them will be 
tangent to both circles. Figure 6.10 shows the 
geometry of an expanding circle moving to the 
left (as seen from the medium). In a time dt, 
the circle moves a distance vdt to the left. The 
boundary will be the loci of all such points as 
the switching procedes. 


Figure 6.10. construction for 
the dynamics of switching 


The parameters used here are the radius r, 
the center of the circle at y = c, and the velocity 
v. The results obtained can be applied to either 
the Karlqvist or the far-field head model. From 
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similar triangles we can obtain that: 
= © = vdt—dz or 


zdzr + (y—c)dy = evdt 


(6.6) 


The equation for the circle is: 27+ (y—c)? = 
r?, which we differentiate to obtain: 


2zdz + 2(y —c)(dy 
zdr+(y-—c)dy = 


dc) = 2rdr or 
ydc + rdr — cd6.7) 
As the left hand members of equations 6.6 


and 6.7 are the same, the right-hand sides are 
equal, so that: 


z= 1 de + nat — oo 
ae ae a 
As r and c are not independent, it is worth 


while to discuss two cases. The simpler is the 
far-field approximation, so it is examined first. 


(6.8) 


6.4 DYNAMIC SWITCHING FOR 
THE FAR FIELD CASE 


For this case r = c, so that equation 6.8 re- 
duces to: 


(6.9) 


Substituting equation 6.8 into the equation 
for the circle we obtain that: 


2r 
14+ (22)? 


v at 


y= (6.10) 


from which we can substitute into equation 
6.8 to obtain z. However, this development has 
been in the head coordinate system, and we 
need the result in the media coordinates which 
are moving to the right of the head system with 
a velocity v. We use the subscript m to desig- 
nate the medium coordinates, so that yn = y 
and z,, = z — vt,so that: 


; 2r 
Spee Lem 6.11 
aR 


From the end of chapter j we have that r= 
gH,/2xH,;, and for this example we look first 
at H, = Hs, and secondly at H,; = Hr to 
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locate the transition region. For H, = Hs, we 
find that 


_ 98, 
ts = on Hs (6.13) 
drs _ g dH, 
“dt Hs dt 4) 
and for H, = H7, we find 
gf, 
TT onHy (6.15) 
drr _ g dH, 
or © osu a 18) 


6.4.1 Far Field Example 


The simplest case to consider is where the field 
changes from its negative saturation value to 
positive saturation as a linear function of time: 
H, = mvt so the loci where a particular value 
of H, has been the largest reversing field expe- 
rienced has radius varying as r = kvtR and 


dr/dt = kuR, where k = mg/2xH,. From 
equation 6.12 we obtain that: 
= Imy yl 
er (k + 7p) (6.17) 
substituting equation 6.17 into 6.12: 
Pe = yank =) (6.18) 


which is the equation of a straight line in the 
medium coordinate system. 

We recall that the recording of predictable 
transition regions requires that the boundary 
locus of the reversed H, = +Hs pattern must 
overlap the previous H, = —Hs pattern. we 
examine equation 6.18 for the situation of H, = 
—Hgs and find that, with ks = mg/27Hs: 


1 


— _ Ym = 
Im = — "5 (ks Es) (6.19) 
and that for H, = +H s: 
Ym ] 
= I aaiiateeahieetiall cee RS 20 


so that as the field begins reversal with the 
H, = —Hs pattern a straight line of a par- 
ticular slope, and as the field reverses the H, = 
+Hs5 pattern is also a straight line with the 
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same magnitude of slope, but of the opposite 
sign. 

The situation where the H, = +Hs line 
just touches the H, = —Hs line is where the 
two slopes are infinite, t.e. when Z,, is zero: 
ks = 1/ks, or k2 = 1. Which requires that 
the field must change with a minimum slope of 
m= 2nH5/Q or that 


H, = 2rHsv 


t (6.21) 


It is of interest to note that, while this de- 
velopment has taken place with un-normalized 
units, the normalization will yield that the min- 
imum slope is given by: 

Hy 


5: 


(6.22) 
where V = 2u/g 

While a linearly changing field is only an ap- 
proximation of practical field switching, this ex- 
ample serves as a guideline to the maximum 
field switching time for predictable transition 
regions. 


6.4.2. Dynamic Switching for the 
Karlqvist Case 


Using the results obtained in equation 6.8: 


_i GC OP. gee 
ee a a de 


For the Karlqvist fringe field approximation 
we found the radius and center of the constant 
H, on page 26 in normalized form: R= esch, 
and C = coth,. Then equation 6.8 can be writ- 
ten in normalized form: 


1 (dC dR dC 
X=o (vs +R - of) (6.23) 


It can be readily shown that: 


dR dC 
ar ae 


Thus we can use equation 6.9 in this case as 
well. Here it is rewritten in normalized form: 
YdR 
V dt 
Substituting equation 6.24 into the normal- 
ized circle equation and using the facts that Y 


KARLQUIST: X = (6.24) 
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must always be positive and R? = C? +1, we 
obtain the following expression for Y: 


C+ /R? + (342)? 
KARLQUIST: Y = Oty Rit (ola) 
1+ (94 


(6.25) 

With C = cscrwh, and R= cot xh, we have 
Gh Se (csc* —— 

dt H? H,° dt 


There is presently no known closed form so- 
lution for the boundaries of the transition re- 
gion for the Karlqvist fringe field case, but nu- 
merical solutions have been performed which 
show similar behavior to the far-field case, ex- 
cept of course in the region very near the gap. 
In that region the linearly changing field indi- 
cates that a greater slope is needed as the field 
passes through zero on its reversal, and less at 
the beginning and end of the switching inter- 
val. Fortunately these are the usual and natu- 
ral properties of magnetic fields in the transient 
state. 
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Chapter 7 


DEMAGNETIZATION 


In chapter 6 we saw how a transition region 
would octur in the absence of demagnetizing 
fields. The transition is a region of chang- 
ing Magnetization, which means that there it 
is a region of finite divergence of magnetiza- 
tion, which in turn means it is a region contain- 
ing magnetic pole density. The pole density is 
the source of demagnetizing field intensity Hg, 
which acts to spread out the transition as much 
as possible. 

The hysteresis characteristic of the recording 
medium limits this spreading of the transition 
region. For any particular level of magnetiza- 
tion, a minimum field is necessary to change 
the magnetization. If the combination of the 
applied field (magnetizing field) and the demag- 
netizing field exceed that minimum. the magne- 
tization will change. 


While the demagnetizing field can change the 
transition region after recording has been com- 
pleted (the applied field no longer present), the 
principle effect occurs during the writing pro- 
cess where the total field is tiie sum of the ap- 
plied field and the demagnetizing field. This 
will result in a transition region shape that is 
altered from those predicted in chapter 6. The 
transition it will be wider (in the z-direction), 
and may even be shifted depending on the dis- 
tances to adjacent transitions. 


Analysis of the recording process, including 
the demagnetizing effects of the recorded mag- 
netization, 1s an iterative problem requiring a 
large high speed computer to carry out self con- 
ststent calculations. This is often done using fi- 
nite element mathematics. a process beyond the 
scope of this course. Insight is obtained only 
by using simplifications to allow the nature or 
characteristics of the process to be expressed 
in analytic form. The general rule is to sim- 


ye 
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plify until only the barest necessary complexity 
remains and so obtain approximate relations. 
After that, complicate as much as is practical 
to gain understanding of secondary effects. In 
this way we can also obtain results useful in 
system design and analysis. 


7.1 THE UNIFORM 
TRANSITION AND ITS 
DEMAGNETIZING FIELD 


The demagnetizing field is due to the diver- 
gence of the magnetization, and can be calcu- 
lated using equation 3.10, which is written here 
in rectangular coordinates: 


H, = j | [pede sdusds, . 


1) 
4nrs, 


where 


rey = 1e(zy — 2s) + 1y(yy — ys) +12 (27 - 25) 


(zy — 25)? + (yy — ys)? + (27 — 25)? 


and l,y =Vresy/Tsy. 

The demagnetizing field is to be calculated at 
the middle of the recorded track, where 2, = 0. 
This requires integration of equation 7.1 over 
all space, which is effectively the space of the 
recorded track in the magnetic medium. To 
begin with, the magnetization distribution is 
not specified, but some specific examples are 
considered after some preliminaries where the 
y and z effects are taken into account by some 
simplifying assumptions. 

In the spirit of simplicity the following as- 
sumptions are made: 


sf 
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1. assume the magnetization is entirely in the 
z-direction, so the divergence of M is sim- 
plified to dM, /dz,. 


2. As the r-components are of primary inter- 
est, the yand z components in rg, are elim 
inated. : 


3. Assume there is no variations in the z di- 
rection and that the track width is much 
greater than any other dimensions, and 
calculations are made at zy = 0. The inte- 
gral limits on the = go to infinity, and with 
no z dependency of the magnetization, the 
integration! over = is carried out to obtain 
equation 7.2, which is the z component of 
the demagnetizing field Hyg: 


"+ — 2,)dr,dy 
Hog = [ [Se ea Be 79 
° 2a (zy — 2s) + (yy — Ys)? a 


7.2. ABRUPT OR STEP 
TRANSITION 


The simplest imaginable transition is one 
which has the magnetization reversing in a zero- 
width transition, which can be written as a step 
function for the case where the magnetization 
is negative to the left of the transition (7.e. the 
magnetization points in the minus z-direction) 
and reverses to positive to the right of the tran- 
sition, as indicated in figure 7.1. 


Mea) Meg 20 Sg) 


where Zz is the location of the transition plane 
in source space. 


The pole density is found from the negative 
of dM,/dz;: 


p(zs) = —2\,,4(r2, — Zo) (7.3) 
where 6(z, — Z,) is a mathematical delta func- 
tion and is distinguished from the medium 
thickness, 6, by the fact that it always has an 
argument [i.e. (r, — z,)}] whereas the medium 
thickness alwavs occurs as an entity [z.e. as 6 
without decorations}. 


I see appendix A. equation ] 
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Figure 7.1. Abrupt magnetization transition 


Then substituting the result of equation 7.3 
into equation 7.2 we obtain: 


Hr = 


- |[ & (zy —2,)dz5dy, 


oa (xy oe Bele Tye = ys)? 


The infinite integral over z, 1s a special type 
of integral involving the delta function, which 
is called the sifting function. This is evaluated 
by replacing the integration by the value of the 
integrand at the point where the delta function 
is non-zero”, t.e. where z, = Zo. This results 
in: 


Hig = — [ Mrs (2s — Zo)dys 
~5/2 ™ (Zy—2o)? + (yy — ys)? 
(7.4) 
The final step is integration over the thick- 
ness of the medium, using equation 2 of ap- 
pendix A to obtain: 


This function has its maximum magnitude at 
ry—-2,= 0: 


|Haz(maz)| = M,, for (-—6/2 <y< §/2) 


At the points yy = +6/2 the maximum magni- 
tude is only half this value. 


2see appendix A section A.2 
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The maximum of the function is found at the 
yj = 0 to be: 


Mrs an72 _ 8/2 


T Zp— Zo 


(7.6) 


At y = +6/2 which includes the near side of the 
recording media, where the transitions seen in 
chapter 6 were the narrowest. The demagnetiz- 
ing field at yy = is: 


\ a ee: 
tan 


ne Tj —TZo 


Hryg=- 


which has its maximum magnitude at ry = Z>: 


|Haz(maz)| ad ele at (Zo, —6/2) 
These demagnetizing field functions at the mid- 
dle and surface of the media are shown in figure 
7.2 with two intermediate heights. It can be 
seen that the surface of the media is a degen- 
erate case, as all other points on the transition 
have a demagnetizing field equal to AM,,. 


Figure 7.2. Abrupt transition demagnetization 
(a) at y=0, (b) at y= +6/4. 
(c) at y = 36/8, (d) at y = +6/2. 
These results provide upper limits to the de- 
magnetizing fields that are possible in a mag- 
netic medium with a transition involving 
the reversal of magnetization of the level Af,,. 
The maximum possible demagnetizing field is 
equal to the maximum magnetization. 


7.38 ARCTANGENT TRANSITION 


The most useful magnetization distribution 
for an isolated transition has been an arctan- 
gent. This model of the transition is written 


2M es wet P ie 
ae {AQ 0 


(7.8) 


a 


R.B. YARBROUGH 


© 9 JuLy 1990 39 


where a is called the a-parameter. Here z, is 
again the arbitrary center for the transition in 
the medium coordinate system, the point where 
M, = 0. At the point where z, —z, = +a 
we find that M, = +0.5M,,. At extreme dis- 
tances from z, the magnetization approaches 


the appropriate saturated value. The deriva- 
tive is taken: 
dM, 2M,; a 7.9) 
a AR es Rs OE Een ee { 
dz, m (%;—-2,)* +a? eee 


which can now be inserted into equation 7.2 
for integration, where the range of z, becomes 
—~O Lr, SC. 


6/2 
moe bff 
ze ~§/2 tm (Zs—Zo)2+a2 


a mies aa, 
————_—_+___—_—_—_—_... dy, 
(xy —Z;)°+ (yy we ys)? 


(7.10) 


This is the form of the convolution integral, 
as explained in Appendix A’. The method of 
evaluating such an infinite integral uses the 
Fourier Transform as outlined in Appendix A, 
section 4. Following that procedure, the follow- 
ing steps are made: 


a 
(rz, —2,)° +a? 


filzs) = 


Using identity I4 in table Al, and transform T6 
in table A2, it is found that: 


Fy(k) = me 2Fok oak 
Next is identified: 


folu-z) = fo(zy—72;) 


ee 
fo(zy — 25) (zy — 25)? + (yy — ys)? 
so that 
rene (z,) 


z? 7 (vy ie ys)? 


Because of the nature of the integral (equa- 
tion 7.10) it is not clear whether the term yy —ys 
or the term y, — yy should be used. Keeping in 
mind that the y-variable of integration is actu- 
ally y,, equation 7.10 can be written in either 


3Appendix A, section 3 
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of two ways, which ordinarily makes no differ- 
ence. However in this procedure it does make a 
difference: 


a §/2 

dates 
ic —z,)dz,d(y, — y;) 
(zy — 25)? + (ys — yy)? 


So that using transform Td from table A2, we 
obtain 


Hzq = 


(7.11) 


Fo(k) = —yre¥-¥2)* for d(ys — yy) (7.12) 
Alternatively equation 7.10 can be written as: 
6/2 co 
hee. = =e 2 f / St eee! a he 
6/2 Joc (2s —Zo)* +4 


i. — Z,)dzr,d(yy — Ys) 
(ry — 2)? + (ys — ys)? 


which results in a transform using T5 again: 


(7.13) 


Fo(k) = gre 4* for diye —ys) (7.14) 


When the entire evaluation procedure is com- 
pleted, it is necessary that the resulting demag- 
netizing field be symmetrical with respect to y. 
As neither equation 7.12 nor equation 7.14 re- 
sult in such symmetry, it is necessary to use 
both to achieve that symmetry. Therefore the 
transformed relation becomes: 


Me of 
Aza = = . 
27° J_$/2 


(ys~-¥e) di seiiae ) a 
J7eE Yf Ys 
| gre dy a) | — 


mei Tok pak y 


This can be manipulated to obtain: 


grelotys—~¥d(yy — v4) | (7 36) 
—yre(atyemyys Ed(yy oan ys) ° 


Using transform TS of Table A2 and iden- 
tity I4 from table Al, the inverse transform 1s 
obtained yielding: 

_ Ms [Pl  (wy-te)dla + ys-ys) 
27 Jogj2 (Zf-T0)? + (G+ Ys-Ys)? 

Mr, [°l? (zy-to)dla + ys-¥s) 

22 Jogja (Zy-Po)? = (at sry}? 
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This is evaluated using equation 2 of ap- 
pendix A: 


ay | -~§/2 
ty = de] eee 
Qe | —tan-} Styrts/2 
Pim-Zo 
1 o+6/2-yy 
. M;; | tan~ hey 7.17 ) 
l — —_ 
25 —tan~ SS 
which can be regrouped to obtain: 
as a-6/2+yF 
Hsq = Mrs tan bie 
° 7 | +tan7? 22f@=y 
te~Zo 


+tan- 1 546/249) fas 


ty~Zo 


M,, | tan7} atéa-y 


Equation 7.18 has the proper symmetry with 
respect to both zy and y,, and in the case of 
a — 0 the demagnetizing field also approaches 
that of the abrupt transition (equation 7.5). 
The highest field occurs on the plane y = 0, 
where the field is: 


Me. = a — 6 2 
——— {tan / 


Arg = 


— tan 


This function has its maxima at (zy—z,)* = 
(a — 6/2)(a + 6/2), which has the value of: 


om a—6/2 

y M,, | ‘aa o+6/2 
dz(maz) = 7 : 

7 | —tan7? ,/sH2 

a-—6/2 


But this can be simplified using the identity: 
tanAstanB 


ni fe ed 
tan(i= B) 1 stan AtanB 
to obtain 
Mrs -1 6 
H ar) = —— tan fae fe 
dz(maz) Tr 4a? = 62 \ 


The demagnetization field is shown in figure 
7.3 for an a parameter of 1.5 6. 

For a situation where the demagnetizing field 
alone set the transition, the a parameter would 
be adjusted:so that the maximum field would be 
approximately equal to the coercive force He. 
This is calculated from equation 7.19: 


a §/2 
~ sin(tHe/ Mrs) 
This can be seen to be a lower limit on the a 
parameter, although effects due to the high per- 


meability of the head could reduce this length 
somewhat. 


(7.20) 
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(. 


Figure 7.3. Arctangent transition demagnetization 
(a) at yy = 0, (b) at y= 46/4. 
(c) at y = 36/8, (d) at y= +6/4. 


7.4 MATERIAL 
CHARACTERISTICS 


The arctangent transition is described by the 
a parameter, for which a minimum value can be 
estimated by setting the maximum demagnetiz- 
ing field Hrgimer) to A, as in eq. 7.20. This 
may be optimistic, because the actual magne- 
tization distribution is determined by applied 
head field H:, and the magnetization charac- 
teristic In conjunction with the demagnetizing 
field. 

An accurate solution to the magnetization 
distribution is an iterative process which re- 
quires numerical] techniques and a computer. 
However, useful information on the a parame- 
ter can be obtained by some simplifications and 
intuitive application of known relations. 

For convenience, it is assumed that the arc- 
tangent transition being written is in the op- 
posite direction from what was discussed previ- 
ously, so the transition is being recorded with 
a positive gap field (in the positive z-direction) 
over a previously saturated region of at negative 
remanance. Then the resulting transition will 
be the negative of equation 7.8, with a value of 
+.©M,,/2 at (zr, — 2.) = ta. 

When the recorded transition is established 
and the transition has moved far away from the 
head, the middle of the transition will have no 
demagnetizing field. However, during the re- 
versal process, the field at that point will have 
exceeded the coercive force (H,) and at that 
time the magnetization and field will be at the 
point indicated in figure 7.4 as point I on the 
dynamic hysteresis curve. When the applied 
field decreases because of movement away from 
the head, the applied field relaxes to zero along 
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the path indicated by yx,. 

In the dynamic process of writing the tran- 
sition, the slope of the M-H characteristic, re- 
ferring to figure 7.4, is: 


OM Mrs _ Mrs 
6H H.-H, H.-S") 


so that S* = H,/H, <1. 


Figure 7.4 Magnetization characteristic 


Again referring to figure 7.4, the middle of 
the transition (z, = z,) will occur at the point 
where the net field has been H, during the 
recording process. The relaxation path in fig- 
ure 7.4 1s from point I to the origin, which has 
the equation: Af=y,H 

The dynamic characteristic in this condition 
can be described as a straight line: 


M = xa(H —H-.) where: 
OM M,;; 


= —— = Tel 
X¢= 9H ~ H-S*) oa) 
From these relations we obtain that: 
Xd—-Xr 


7.5 FIELD GRADIENT 


During the recording process the characteris- 
tic of the material and the total magnetic field 
interacts so that a particular magnetization dis- 
tribution results in a demagnetizing field which 
alters the total magnetic field, resulting in a 
changed distribution. This can be treated by a 
series of calculations of the magnetization dis- 
tribution, which then gives rise to a demagne- 
tizing field which adds to the field, so that anew 
distribution must be calculated. This is an it- 
erative process which is call the self-consistent 
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calculation method, and consumes much com- 
puter time. Here we simplify by recognizing 
that the derivative of the magnetization must 
follow the following condition: 


dM, _ 0M, dH 
dr OH az 
but H = Azgt+ Aza 


where Hyg is the negative of equation 7.18, and: 


oe tan-! znitt : : 
y 


From the negative of equation 7.9 at the mid- 
dle of the transition (z, = 2.) gives a slope of: 


(7.23) 


dM QM, 


dr TQ 


At zs = Z, the demagnetizing field is zero, so 
the applied field must be set to H;. However 
it is desirable to make the slope of the applied 
field, dH,,/dz as large as possible so as to ob- 
tain as small an a parameter as possible. 

Equation 7.24 is first differentiated, and then 
the second derivative obtained and set to zero 
to find the point of the maximum slope, which 
is the maximum of dH;,/dz: 


dHee _ Hy | gee | (7.23) 
dz T —- ——_ + 
(me g/2)?+y? 


Taking the second derivative with respect to z 
and setting the result to zero, the following ex- 
pression is obtained: 


r+q9/2 _ zr—g/2 
(r+9/2)?+y? (z-g/2)? +y? 


solving this for the position of maximum ap- 


plied field slope: 
“| - W@d-lal 
—- =\/j/—-| +|—-| +1 
Fa 3 bs g/2 

] y : 

3 Fal 


which will become the center of the transition, 
Zo 

At the center of the transition the field is also 
to be set to H, as given in equation 7.22. From 
this the value of the gap field can be determined 


(7.26) 
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depending on which point in the media is se- 
lected to receive H,. 

The calculations are to be done at the mid- 
dle of the recording media, which in the media 
coordinates is at yy = y, = 0, and in the head 
coordinates is at y = d+6/2. Then taking 
the negative of equation 7.9 at z, = z, for the 
left side of equation 7.23, using equation 7.21 to 
obtain 0M/OH, dH,,/dz from equation 7.25 at 
z=27, and y=d+ 6/2, and finally obtaining 
OH-.q/Ozr from the derivative of equation 7.18 
with y, = 0 and z, = z, to obtain the following 
relation: 


Dies M;s 


ma ~=H.(1—S") : 


d+6/2 
(rotg/2)2*+(d+6/2)? 
d+6/2 


~ (Fong/2)*+(d+6/2)" 


a Afes 1 ma 1 
7 e a—6/2 


As 
v 


(7.27) 


where H, is adjusted according to equations 
7.22 and 7.24 to give: 


H 
Xd— Xr 
-—1 Zot+g/2 
H, ban” ” "33372 2 
— —tan7! ro-9/2 ( 28) 


Equation 7.27 can be evaluated for any specific 
geometry. The reader is encouraged to evaluate 
several examples. 
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A INTEGRALS AND 
TRANSFORMS 


AA.1 STANDARD INTEGRALS 


The integrals here can be found in most ta- 
bles of integrals. They can be verified by differ- 
entiation. 


| du ] u 
[oS oe (1) 


du l u 
Cae 
udu 1 , 2 . A 
i = 5 In(u + Uv ) (.3) 


Ai.2 THE SIFTING FUNCTION 


The delta function, 6(u—), is also known as 
the sifting functzon because it has the property 
of selecting out the value of a function at the 
point u = 6 whenever it multiplies the function 
and they are integrated over any interval that 
includes the point 0. 


[ub ujdu = £0 fora<b<c (4) 


A.3 THE CONVOLUTION INTEGRAL 


The form of the convolution integral is as fol- 
lows: 


Peale | “fOpe-oe. © 


This form of integral occurs in the analysis of 
demagnetizing fields and of play-back voltages 
in magnetic recording systems. Two interesting 
forms of the convolution integral are: 


fe) =f fledhley - eae, 


f(vt) 


[fla \p(et ~ 2)az 


In principle these are readily evaluated us- 
ing Fourier Transforms. However, these partic- 
ular functions (and others of possible interest 
in magnetic recording) are not usually found in 
most tables of transformations. For that rea- 
son, a special table of such transforms are de- 
veloped in the following. 
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f1.4 FOURIER TRANSFORMS 


Transform pairs for various functions are de- 
rived using the fundamental definition of the 
transform and inverse transform. The term 
1/27 must appear either on the transform or 
inverse transform, or be split between them 
as 1/,/27, and the placement is not consistent 
from one table of transforms to the next, so 
it is important to note its placement when us- 
ing a particular table of transform pairs. Here 
we place it in the definition of the transform of 
F(k) to f(z), which is usually called the Jnverse 


Transform: 


co 
jee: F(k)eit*dk (6) 
2n k=—0o 
where f(z) is a function of z, F(k) is its Fourier 
transformation into k-space, and j = /—1, so 
that we are dealing with complex variables in 
k-space. 
The transformation from z-space to k-space 
is usually called the Transform: 


oe 
F(k)= [f(z (1) 
rt=— 6 

The derivations and proofs of the properties of 
Fourier transform pairs is beyond the scope of 
this work. Here some of those properties are 
provided as aset of fundamental transformation 
identities, which are given in the table Al. 


TABLE Al TRANSFORM IDENTITIES 
I f(z) 


2 Af(z) 


13 filz) + fa(z) 


14 f(z+c) 


15 Wiz) 


dz 
16 fF, f(ujdu 


17) fr f(z) fo(u — z)dz 


The procedure in evaluating convolution in- 
tegrals is to begin with identity I7 which has 
the convolution integral on the left side. We 
must next find the transformations of f,(2) and 
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fo(z): Fi(k) and Fo(k) respectively using table 
A2, the table of transform pairs. 

The form f2(u— 2) is essential, i.e. the func- 
tion must be written so that the integration 
variable occurs in a subtractive relation to the 
independent variable of the function being eval- 
uated. Thus f. must be a function of vt — z or 
zy —2z,. However, when finding Fo, it is to be 
found from f(z) or fo(z,). 

The next step is to manipulate the product 
F,(k)Fo(k) to the form of some function of k 
on the right hand side of table A2, which then 


allows the evaluation of the function of vt or | 


Ty. 
F(k) 
ej ak 


mae 
jm[5(k + 3) 
-6(k - 8)] 
n[5(k + B) 
—+6(k — B)] 


jreluk 
ore] uk 


= pjuk 
pk 


x pjuk 
EE 


276(k — a) 
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Magnetoresistive (MR) Sensor 


VOLTAGE, — 


4. 
7 V= iA, — 1AR(H) 
Sense at . Bat te AR =) 
CUCPLLENT =e 0 Roo 


unbiased response biased response 


FILM lS AMISO T2Oe | C 


f 


- Re may whe Mil CL 


© M voteteg ww pve sine of 


sw ® e(H\ . We 2K 
Ho) Mes, 


° Ae So es oF bun OS LL Venn Sangin 
| ) ) | TAG 
r | ve Kei a aa ow { 


Des ee ape | i . A 4g 
é ees a ae 4 oe WWR UNT OY 


Advantages 
e Large signa! output 
e Velocity independent signal amplitude 


@ Gocc linear resolution with shielded 
MR sensor 


Basic Issues 


e Suopressicn cf Barkausen noise 


e Linezrizeatics cf MR response 
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Operating Magnetic Environments 


Impesed By requirements of:. 


e Line=r Fesolrrtion: 
— Presence of soft-meécnetic shields 
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” MR Element (Edge View) 
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The BARKFRAUSEN NOISE Problem 


© Domai am yvities im small MR elements 


° Relater Setween domain activities & noise in MR respocnss 


° Origin of camaings 
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Dispersion of the Anisotropy Easy-Axis 
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Techniques for_Generating Transverse Bias Field: 
(a) (b) (c) 


Bias 
RMR Conductor aN r Permanent 
Direcsoar af 


Direction of Current 
Bias Current 


Sturt Soft-Film Permanent. Mzsnet 
Biasing Biasing Biasins 
e Week effects e Strong effects e Selection cf 
permeability limited pecmanent magnet 
tilm important 
e Electrical shunting ¢ Thin spacer for e Good bias grefite 


uniform bras prafile 
e Nanp-uniform bias e¢ Selection of soft: 
| if . 


prone 
important 
e Enhancements e Saturatec scft-film 
operations 
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MAGNETIC RECORDING CHANNEL 
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Conliguration of Codes 


Error 


Modulation 


User | , Modulation ; User 
Data —? | Correction Code Channel |—> Code Correction ae 
Code 
Encoders . | - Decoders 
Modulation Code: Matches recorded signal characteristics to channel 


bandwidth, delection method, read/write electronics, 
timing and tracking servo requirements 


Error Correction Code: Detects and corrects data detection errors 


Digital Magnetic Recording Channel 


Magnetic 
Track: 


Data: sO oe F « FT a OF 3 NRZI 


Write 
Current: 


Readback 
Voltage: 


‘Detected — | | 
bas lk Bw Fc ts 


Causes of bit detection errors 

e Random noise 

e Intersymbols Interference ‘| Bit pattern 
- ® Loss of clock synchronization related 
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# =Current reversal 
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$ = Pulse peak spacing 
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Pulse Crowding 
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M(x) H(x) 


voltage 


(b) 
a 
7, Two views (a) or (b) 
INPUT WRITE READ . OUTPUT 


FIGURE 7.4 


. Recording channel wavelength response 
(transfer function) 


M’x(k) = jkMx(k) 


jk represents the differentiation on readback 


and 
1-e-ké sin(kg/2 
Hx(k.d) - (—_—}-e-Ke — 
thickness spacing gap 
factor loss 10Ss 
OT 


E(k) = jkwoNvws (1-246 -kd Sinks Mx(k) 


N = number turns 

Vv = velocity 

W = track width 

& = medium thickness 


HIGH DENSITY RECORDING 


os Narrow PW ,, 
fA Required 


PWe = \/ g + 4 (a+d) (a+d+ 6) 


0 = gap length 
a = transition parameter 
6 = media thickness 


_ 2M, 8 
H. 
H 


a 


¢ = Media coercivity 


; = media remanent flux density 


Source: Williams - Comstock (6) 
MIG - Sutton vans Applieo Magnetics Corporation 
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Linear Density: Longitudinal Recording 
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DISTORTIONLESS REPRODUCTION 


vi(t) CIRCUIT volt) 
OR 


CHANNEL Vo(w) 


TRANSFER FUNCTION 


H(w) = Vo(w)/Vi(w) (independent of signal) 


AMPLITUDE RESPONSE 


a FLAT OVER 
FOR A BANDWIDTH 
DISTORTIONLESS OF INTEREST 
CHANNEL 


FREQUENCY 


> 
| — STRAIGHT LINE 


e.g., Audio amplifier, telephone channel WITH ZERO INTERCEPT 
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EQUALIZER 


| RESPONSE va 
~~ J 
; ‘\ CHANNEL RESPONSE 
: \ / FLAT EQUALIZED 
/ \ RESPONSE 
ul / 
-Q 
: = 
f - \ 


FREQUENCY — 
EQUALIZATION OF CHANNEL RESPONSE 
(AUDIO RECORDERS, INSTRUMENTATION 
RECORDERS) 


Saturation Recording 
-© For Track Density Change 


— Signal « as 
tpl 


— Media Noise « 1 
J tpl 


or SNR decreases 3db as double tpi 
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THREE MAIN SOURCES OF NOISE _ 


\ TOTAL 
\ geen ——_- TAPE 
"SLOT" 4 
NOISE We y 
POWER | — 
(dB) > HEADS 
ELECTRONICS 
| dc FREQUENCY 
ORIGIN OF NOISES IS ALWAYS SOME 


UNCERTAINTY ABOUT A PHENOMENON: 


ELECTRONICS: ELECTRON (NOISE) STATISTICS IN FIRST 
STAGE REPRODUCE PRE-AMPLIFIER 


RESISTORS: ELECTRON STATISTICS 
HEADS: DOMAIN WALL/ANGLE STATISTICS 
TAPES: TAPE PARTICLE POSITION OR ORIENTATION 


STATISTICS 


JOHNSON NOISE 
RESISTORS < 


NYQUIST THEOREM 


v=V/4kTR Af) VOLTS 


k= BOLTZMAN’S CONSTANT 
(1.38 10°79 JOULES/DEGREE) 


T= ABSOLUTE TEMPERATURE 


R= RESISTANCE IN OHMS 


Af= BANDWIDTH OF MEASUREMENT 


“, 


CPS 


ELECTRICAL IMPEDANCE OF ATOROID 


CROSS SECTION 
AREA A 


d AV. PATH LENGTH, L 


PERMEABILITY ,s° = Hy—)uy 
eS 


INDUCTANCE, L* — [0.4 wv 107-8} N2 ab 


N2A N2A 
IMPEDANCE = JuL* = Jo NWA 4 gy ANTA 
j L Lp L tad | 


PURE INDUCTANCE RESISTANCE __ 
STORES ENERGY _ DISSIPATES ENERGY 


.. NOISELESS .. NOISY 


IF PARTICLE: 
POSITIONS ARE RANDOM 
DIRECTIONS OF MAGNETIZATION ARE RANDOM 


THEN, TO FIND THE TOTAL NOISE POWER, WE 
SIMPLY ADD THE POWERS 


IF E, ANDE, ARE TWO SIGNALS, 
ar ae ay ee 

(E,+E,) Er +E5 + 2E,E5 

IF E, ANDE, ARE RANDOM E,E, = 0 


UNCORRELATED 
INCOHERENT 


THE NOISE POWER IS THE SUM OF THE SQUARES 
(NOT THE SQUARE OF THE SUM) 


(sf =——: 


WHY IS IT THAT: 


SNRa nwa? 277 
10,000 — 40dB 
1,000 — 30dB 


THE REPRODUCE HEAD ONLY SENSES A VOLUME 
OF TAPE, 1/2 A WAVELENGTH LONG AND ABOUT 
1/5 \ DEEP, AND THE TRACKWIDTH WIDE: 
VOLUME aWa2 


.*. ALWAYS BETTER TO REDUCE W THAN REDUCE All 


MEDIA NOISE 


-e Noise in thin metallic media is associated main/y with 
written transitions. : | 


e Position jitter and amplitude variations contribute to total 
noise. 


e In general, maximum noise occurs at a density which is 
related to the transition width (A, ~Aa, Belk et al). 


e Two regimes, (i) noise initially increases linearly with 
density followed by (ii) “supra-linear” increase as 
transitions crowd together. | 


Track Divection, ym 
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CONTROLLING MEDIA NOISE (2 


com 


e Reducing the coupling between grains reduces the 
coercive squareness of the hysteresis loop, and 


decreases the irregularity in the transition boundary. 


1.0 
6 /a=0.75 
U/H,=0.4 
0.5 a=Q.i 


magnetization M,/M, 
© ] 
ro) 


—0.4 —0.2 0.0 0.2 0.4 
applied field H/H, 


(Zhu and Bertram, J. Appl. Phys. 63 (1988)) 
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Magnetic recording process 
| apecnia* 


Input eee of Output (voitage) 


recorded magneuzation) 


(a) 


Electncal filter 


TUw) 
TUw) = BUw) GUw) ° 


(b) 


Fig. 6.20. Input-output trensier function block diegrams. 


been defining H(z) for —o <2 < +, end since in our anslogy = (or= 
corresponds to ‘“‘time,’’ we see that our impulse response transier Iunction 


provides an iniinioateey OUuTDUL. 


Now we can moaiy the input-output reiziion snown in Fig. €.20c by 
insertion of an electrical Slier et the head outpur. Lf the voltage at ine 
input terminals of this Dlier is nov disturbed Dy Its presence, we can re- 


Craw the input-output transier relation, as snown in Fig. 6. 


G(sw) represents the irecuencr characteristic of the slewicinal ali 
Sltuzvion will epply, ior A. represents the irensier Iunction of 7 


200, nv wher eS 


metic Tecorcing process. and the ree parameters Of the mernedc 
heed can be viewed 2s incorporeied into the Gesign of ine inserted Slter. 
é tns iollow- 


Line input Af’(s) is then reistec to the output voltage e(=) DY i 


Ing transfer function 
Tz) = BGe)G Qe) 
Theoretically, if Gijz) could de chosen such thet 


Zun)EGe) = 1.0 


ee) = AL (=) 


| | es lee og 
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tone) Dendwidsh, then E(i) woule be identicsliy ecuel to 
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bi” \. __ FINAL SLIMMED PULSE 
’ =a+b+c 


eee TIME. 

a: C 

COSINE EQUALIZER tT OT 
| DELAY DELAY 
PLAYBACK T T 
SIGNAL | , : 

| TA 

-1/4 & ~1/4 


OUTPUT 


San 


PULSE SLIMMING 
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RLL Code Classification 
(after Franascek) 


for instance MEM: 


d km ner 
1 3 1 2 _ 


Number of word lengths 


(Variable word length code: 


only) 


Number of code cell per bit 
(=number of possible flux 
reversal positions per bits) 


‘maximum number of bits 
required for coding at 
any one time 


m/n: ratio data to code bits 


maximum number of code cells 
between flux transitions 


minimum number of code cells 
between flux transitions 


Source: N. Mackintosh 
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Comparison of RLL (d,k) Constraints 


een b----- - 
Data: 1 . Gg 1 
r— b---4 
1/2(0,1) ooh % Fe dl eo se Oe De Pe 
| cane aaa. 
1/2(1,3) cs Od oe Oa 8 oO 4 ta OD « 
oo freee P==-" 
1/2(2,7): eV «On Oe Tx Oe eo Oe Ox 
-}—____—__— b--—-4 


C3): «alt « © we 1+ + Oz Ba V 


| - (2,7) { Less pulse crowding +——4 
Sur ale ~ (1,7) Same/larger window #-----4 


(2,7) vs (1,7) | Involves tradeofis 
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aS 


(2,7) vs MEM (14,3) 


Fixed User Bit Density 


MFM. 1. O04, Trin = 2 - 
Ay «te Oe Wn fs * Ties’ 


(2,7) reduces intersymbol interference 


Fixed Minimum Transition Spacing Tala 
—_ min 
Peete 1 | | : 
MFM. . 1.0.1, , ' Density Ratio 
2x 1/2 = 1 bit/T,,, 
| | | b---4 . 
(2,7) ct « O « O & x | Density Ratio 


3.x 1/2 = 15 bit/Tmp 


(2,/) increases density ratio by 50% 
-with decreased detection window 
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CHANNEL CODE PARAMETERS 


T = DATA BIT TIME INTERVAL 
= +t OF DATA BITS IN GROUP 


n = # OF CODE BITS IN GROUP 
m — CODE RATE 

d = MINIMUM # OF ZEROS BETWEEN ONES 
k = MAXIMUM # OF ZEROS BETWEEN ONES 


Tin = MINIMUM TIME BETWEEN TRANSITIONS 
Trin = mMeUT 


Tmaz = MAXIMUM TIME BETWEEN TRANSITIONS 
T ros = —— 


W = CLOCKING WINDOW 
W="2T 


DR = DATA DENSITY/TRANSITION DENSITY 
Drs Tain = m(d+1) 


LOW FREQ ALL ONES/HIGH FREQ ALL ONES 
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Practical Example: Rate 2/3, (1,7) Code. 


00 101 

01 100 

1Q 00] 

a 
Basic Coding Table 


-  exaGaaeeeeinnnematentennc eae niCn AEN at INARA N OP AC: CE COLE OE rennet 


Data. Code 
00.00 101.000 
00.01 100.000 
10.00 001.000 


10.01 010.000 


Substituting Coding Table 


it | 


DOUBLE =| 
I 

rreuency ; - i 2 | ¥@ oo 1 V2 
H | 

MFM . 3 —_ 


ZM 1 2 1/2 


| ) 
6 
8 
5 


2/3 i ] NT/3 | 16T/3 
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SIGNAL DETECTION 


DataStorage'90 


SPURIOUS 
TRANSITIONS © 


t | 
| || | 
| | ZERO-CROSSING 


PLAYBACK 


SIGNAL 
___ [LOW PASS 
FILTER om 


|. 
FULL-WAVE 7 | - 
/ RECTIFIER | THRESHOLDING |—e 


ee ERD CR ONCE SY SNE 


ISOLATED PULSE 


DETECTED 
TRANSITIONS 


GATED PEAK DETECTION 


(sf —— 
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DATA DETECTION 
Window 


Tw 


a x 
eee gC el 
Detection : is 
threshold eed 
| a: 7 a i 
: \ Noisy 
; i signal 
4 re > 
4 x = 
- rms noise 


G, = Signal slope within Ty, 
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RITE CURRENT >——— 
WRITE DATA —>——_ | -PIG. B, #1 = ICout FIG. B, #2 = anki 
“SAD, with: N, L, -——— ———_—_, ——_—__ | 
, etc. >4 HEAD- | | AMPS | | NOISE | KF 
>4 ARM p——————>4 AGC -—>4 FILTER -—>———, 
™e : | R/W AMP! 1 & tI | i | 
| "CHIP" -——>4 ETC. F—->4 (LPF) ->— ! 
3 beeen a | | 
da | 
oe DISK ‘ wi th Mr p He : = ; @2c.. | 


< 
[nr 
| a I a aaa | aa a aa a i 
| | “GATING" | BI-DIRECTIONAL ie PIG. B, #4 = TGout 


Came ot CEQUALIZER: ero. “THRESHOLD=GAIED” = 


[ 
l 
| 
| 
I | | ' | PEAK QUALIFIER | | DE leCs 2D 
>—— |! >4 H-———=- >i (GATE GENERATOR) | —_ RAW-READ 
| | ene DATA 
| ! GATE DETECTOR or "PEAK QUALIFIER" PATE ee 
| FIG. B, #3 = DFout ———>1 HIGH= | 
| ] oe : | SPEED 
| | | "PEAK" | | DIFFER- | | BI-DIRECTIONAL ! i — 
| bee >} EQUAL-—- Fem—o4d ENTIATOR we 4 ZERO-CROSSING }———>- "AND" 
: | SZER. | | | PULSE | | CIR- 
wen > 6 (EQ. ) rat 6) fat ae GENERATOR | L CULT | 
a | (nsnsssnessensenneenmnancupstonasianamamenunemsennnattl {hesnannsassnsetsnnssnsssesfessceSNsesesneSeNSANRRNS hencrseaneeatensoacnennnnnsmennma, 
PEAK DETECTOR PATH 
» nee 4 
re «6 SCERARAT ED 
| FIG. B, #5 = RRDout [ULTRA-EIGH SPEED ! DATA 
>4 "AND" CIRCUIT FIG. B, 
———— nn PTS. B, #6 = Wouty & [#7 = SDout 


! RECORDING-CODE ———————————"————————>- PULSE STRETCHER >" 
s—>1 "DATA-CELL" or | CLOCKING SIGNALS ee eer 
| | "WINDOW" GEN. aa a aS 


I 

: 

{ 

| 

lL | Re FIG. B, #8 = CLKout ” 
I 

| 

| 


< 

Perce eee es rer nn aaa acacia erate : 

| PHASE | | SAMPLE & HOLD |! | FEEDBACK |! |! VOLTAGE- 
Ke—— >} DISCRIM- |! | CIRCUIT, & | | COMPEN- | | CONTROLLED 

| INATOR p——e >>] INTEGRATOR pam >of SATOR me >Y OSCILLATOR ~——c 
——>i (time-to- | 1 (current-to- | | 

1 current) ! ! voltage) | |(lead/lag) | | (VCO) 

a <iiiiealsaaanceeeeaaed I scessosinesiresisesaiiiaceemnisiatdl Iacansssseensesstenessnsnnnsnnasvunenetsaanamante! 


ee a ee gen ee 
<————1 PROGRAMMABLE (DYNAMIC) COUNTER |! 
ODE PATTERN LOOK AHEAD ! be ne © mm 
O DETERMINE "nn" o> «=(DIVIDE VCO FREQUENCY BY "n") | 


"VFO", "PLL", "PLO", or PHASE-LOCKED DATA SEPARATOR CLOCKING SYSTEM 


| FIGURE A. TYPICAL HIGH-PERFORMANCE DISK DRIVE RECORDING CHANNEL 


F. J. SORDELLO i9&&é 
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| tRADBACK TRANSITIONS for 2, 7, RLL 
SODE ->0 O 000000014100 000000010 0 0 0 i 0 0000 


Fl ICout 
INDIVIDUAL oil at 
ISOLATED Be 
PULSES /,4 \! ! a 
ae | ey, e lacs 
LPFout — | ——_ ae = 
eT NA LO 
; y, ; 
F3 DPout ; oe ; | - 
"SOLID" : oe 
} 
F2 LPFou \ 
"DASHED"! are 
#3 DFout | —— 
| "SOLID" 2 | 
| 
! 
#4 TGout Leeemiemetinattinttninenncernnnnenrcmed | haan a | ea | 4 ee 
| | 4 | 
#5 RRDowt— | $$$ 0 
een I ee BO ce cee ce cee cree ce Se ee OT ce re ee 
#6 Wout Mammen amd dd aed aed | Led sd | Ld et fmt 
#Z SDout eee 
3 | 


PRR MRR A RMR RAR AMAA eA eee 
Sa a ee ee Oe Oe ee Oe DO Oe Oe Oe ee ee ee 


CODEout 0/0: 1000000010 0} 0000000:10 of2 O O13 0000 
eto 


om?" RRADBACK SIGNALS AT VARIOUS STAGES WITEIN READ CHANNEL. 
{All intrinsic circuit delays have been eliminatec fer 
Clairity. For 24 Mbit/sec = 3 Mbyte/sec, one "Wout" 
cell (% complete cycle of "Wout") = 20.833 nanoseconc. 


F. J. SORDELLO 1988 
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% AMPLITUDE OF UNWANTED ADDITIVE SIGNAL 
BIT SHIFT ="-------------------------------------------------- = nsec 
SLOPE OF DFout SIGNAL NEAR "ZERO-CROSSING" POINT 


IT IS AS SIMPLE AS THAT! ANY UNWANTED SIGNAL THAT "LIFTS" OR "LOWERS" 
THE DESIRED RECORDED INFORMATION BEARING DFout SIGKAL WILL CAUSE THE 
"ZERO-CROSSING" OF THE COMPOSITE DFout SIGNAL TO BE SHIFTED IN TIME FROM 
THE IDEAL "ZERO-CROSSING" TIME. THE TIME WIDTH OF THE DATA TRANSITION- 
CLOCKING WINDOW, Wout, DETERMINES THE TOTAL AMOUNT OF "LIFTING" OR 
"LOWERING" THAT CAN BE TOLERATED BEFORE THE COMPOSITE DFout SIGNAL HAS 
ITS "ZERO-CROSSING" SHIPTED IN TIME BEYOND A PARTICULAR TRANSITION OR 
"BIT" CELL. FIGURE D SHOWS THAT THE 24 MEGABIT/SECOND RECORDING SYSTEM 
DESCRIBED CAN TOLERATE COMPOSITE DFout “ZERO-CROSSING" TIME“-SHIFTS ("BIT 
SHIFTS") UP TO + (20.8 nsec)/2 = + 10.4 nsec BEFORE OCCURRING OUTSIDE OF 
THE TIME BOUNDARIES OF Wout. FIGURE D ALSO SHOWS THAT THE + 10.4 nsec 
Bit SHeirlt WILL BE CAUSED BY TEE ADDITION OF UNWANTED SIGNALS EQUAL TO 
APPROXIMATELY + 80 ® OF THE BASE-TO-PEAK VALUE OF THE INFORMATION 
SEARING DFout SIGNAL. 


THE ASOVE NUMBERS ARE ASSOCIATED WITE A SINGLE ISOLATED READSACK 
TRANSITION. TEE AMPLITUDE OF THE BASE-TO-PEAK VALUE OF DFout CONTROLS 
TEE ASSOLUTE AMOUNT OF "Bit SEIFT" TEAT OCCURS FOR A GIVEN AMOUNT OF 
JINWANTED ADDITIVE SIGNAL. IGEN PORTIONS OF ISOLATED DFout: PULSES 
-OMBINE DUE TO TEETR PROXIMITY, TEE BASE-TO-PEAK VALUE OF IEE N=t 
RESULTING DFout SIGNAL MAY INCREASE BY AS MUCE AS 50%, YIELDING LESs 
A'SIT SHIFT" FOR A GIVEN AMOUNT OF UNWANTED ADDITIVE SIGNAL. SEE FIGURE. 
» ON PAGE 3, AND ESTIMATE THE POSITIVE PORTIONS OF TEE DFout SIGNAL 
~ORRESPONDING TO TEE SECOND, THIRD, FOURTE, AND FIFTE TRANSITIONS; AND 
THE NEGATIVE PORTIONS OF TEE FIFTH AND SIXTH TRANSITIONS. TEEN, COMPARE 
©HE ESTIMATED PORTIONS STATED ABOVE WITE THE NEGATIVE PORTIONS OF TES 
OFouT SIGNAL CORRESPONDING TO TEE FIRST, SECOND, TEIRD, AND FOURTH 
TRANSITIONS, AS WELL AS TEE POSITIVE PORTIONS OF THE FIRST AND SIXTH 
TRANSITIONS 


TEeres ARS CISs8 WAYS TEAT CHE "ZERO CROSSING" CAN OCCUR SHIFTED. MOST 
PROMINENT IS BY WRITING TEE ZERO MAGNETIZATION POINT IN THE MEDIUM IN AN 
JNWANTED, SHIFTED LOCATION. EEAD/MEDIUM COMBINATIONS WITH POOR OVERWRIT: 
{CAPABILITIES, COMBINED WITH FINITE WRITE CURRENT AND/OR EEAD GENERATED 
WRITE FLUX RISETIMES AND PREVIOUS RECORDED EISTORY IN THE MEDIUM MAY 


*AUSE TWO OR TEREE NANOSECONDS BIT SHIFT. 


HE MAJOR CAUSES OF DIGITAL MAGNETIC RECORDING CEANNEL SIT-SEIFT ARE 
EE FOLLOWING: 


ae NOLoe. 2hoVone bee tenes s (Nae) 4 
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on CVESWroos shove ee Seem once. (Ons) 
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4. ASYMMETRY INDUCED BIT-SHIFT (ASIB). 
5. SYSTEM NOISE INDUCED BIT-SHIFT (SNIB). 


6. MINOR MEDIA DEFECTS THAT PASS DEPECT SCANNING TESTING BE- 
CAUSE THE MINOR DEFECTS ARE WITHIN SPECIFICATIONS. (MMD) 


7. ADJACENT TRACK INDUCED BIT-SHIFT (ATIB). 


NIB IS CAUSED BY MEDIA, HEAD, AND ELECTRONIC CIRCUIT NOISE SOURCES. THE 
INSTANTANEOUS AMPLITUDE OF THIS NOISE IS COMMONLY ASSUMED TO BE 
DEPENDENT UPON THE NOISE SOURCE'S RMS VALUE, MAGNIFIED BY THE SIGMAS 
ASSOCIATED WITH THE PROBABILITIES OF GAUSSIAN DISTRIBUTION OF CHANCE. 

IF THE NOISE IS ASSUMED TO BE GAUSSIAN, IT LENDS ITSELF TO THE FOLLOWING 
PROBABILITIES: 


RMS VALUE = SIGMA 


SIZE OF NOISE "SPIKE" IN PROBABILITY OF OCCURRENCE 
TERMS OF NOISE RMS VALUE OF THAT NOISE "SPIKE" 

0 1/1 

1 1/3.2 

2 1/21.74 

3 1/370.4 

4 1/15,625 

5 1/1,724,137.9 

6 1/500,000,000 

6.47 =" 6.5 1/10,000,000,000 

7 1/400,000,000,000 

8 1/80,000,000,000,000,000 


FIGURE D SHOWS THAT IF THE 24 MEGABIT/SECOND EXAMPLE SYSTEM HAD SNR 
EQUAL TO 26db (BASE-TQ-PEAK SIGNAL DIVIDED BY RMS NOISE) AT THE OUTPUT 
OF THE DIFFERENTIATOR , THE COMPOSITE "ZERO CROSSING" WOULD SHIFT IN 
TIME WITH AN RMS JITTER EQUAL TO: 


26db = 20:1; THE RMS NOISE IS 5% OF THE BASE-TO-PEAK DFout. 


HENCE, THE "ZERO-CROSSING" WILL SHIFT: (5%)/(10%/nsec) = 0.5 nsec. (RMS) 


s | 
NORMALLY, SNR IS SPECIFIED AT THE HEAD'S OUTPUT OR THE HEAD'S IC 

PREAMP'S OUTPUT. IN ORDER TO OBTAIN THE SNR AT THE DIFFERENTIATOR'S 
OUTPUT, THE RMS NOISE OR THE SPECIFIED NOISE SPECTRA CHARACTERISTICS ARE 
OBTAINED FROM THE SPECIFIED SNR TEST CONDITIONS, THEN POWER SPECTRA 
MODIFIED BY THE INVOLVED READ CHANNEL BLOCKS SUCH AS AGC AMPLIFIER, 
NOISE FILTER, "PEAK" EQUALIZER, AND, OF COURSE, THE DIFFERENTIATOR. THE 
RMS NOISE PRESENT AT THE OUTPUT OF THE DIFFERENTIATOR OBTAINED 

CAN BE USED EITHER TO CALCULATE SNR OR DIRECTLY TO PREDICT NIB. 
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BIT-SHIFT BIT-SHIFT 
ON TRACK ON TRACK and DURING MINOR DEFECT 
10 (DFout slope divided by 0,3) 
(@ ANY TIME @ x10 @ ANY TINE @ x10") 
(6.5*sigma) (6.5*sigma) 
NIB 0.71 3420 1.01 4.64 
PIB 2.4 2.4 3.43 8.43 
OWIB I 0.8 0.8 0.8 0.8 
OWIB II 1-> 3 1 -> 3 1 =-> 3 1 -> 3 
ASIB 0.5 -> 2 0.5 => 2 0.5 -> 2 0.5 -> 2 
SNIB O.1 -> ? O.1 -> ? 0.14 -> ? 0.14 -> ? 
ATIB ~O “0 “0 “0 
5.51 8.05 6.88 9.71 
ERROR | 
RATE to7° 107° . 107° 107° 
MARGIN 4.91 2.37 3.54 0.71 
BIT-SHIFT BIT-SHIFT 
| OFF TRACK OFF TRACK and DURING MINOR DEFECT 
(DFout slope divided by 014°) (DFout slope divided by 05370-7275) 
| (@ ANY TIME @ x10 @ ANY TIME @ x10) 
(6.5*sigma) (6.5*sigma) 
NIB 0.95 4.33 1.35 6.19 
PIB 2.4 2.4 | 3.43 3.43 
OWIB I 0.8 0.8 0.8 0.8 
OWIB II 1-> 3 1-> 3 1-> 3 1 -> 3 
ASIB 0.5 -> 2 0.5 -> 2 0.5 => 2 0.5 -> 2 
SNIB 0.13 -> ? 0.13 -> ? 0.19 -> ? 0.19 -> ? 
ATIB 1.5 5 2.14 2.14 
7.28 10.66 9.41 14.25 
ERROR 
RATE <67° 10° 107° 10 
MARGIN 3.14 fe) 1.01 fe) 
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Visualizing the shifts in detection that occur due to noise and Peak Shift as a Statistica! Distribu- 

~ tion within the Data Window is the first step in the modeling process. As shown in Figure 2, the 
perspective gained by this model yields an immediate observation, the error rate of a drive is the 
total area of the distribution falling outside the Data Window. 
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FIGURE 2 


STATISTICAL DISTRIBUTION 


The next step in the modeling process is the choice of statistics. Assuming that the shifts in peak 
detection conform to a Gaussian model, makes the task of comparing noise and Peak Shift to the 
Data Window simple. According to the Gaussian model, the area outside the Data Window is di- 
rectly related to the ratio of Data window width to the Mean (1) and Standard Deviation (co) of the 
distribution. This is referred to as the intrinsic error rate T of the system. 


+ Katz and Campbell, °Effect of Bitshift on Error Rate in Magnetic Recording’, IEEE Transactions on 
Magnetics, Vol. Mag. 15, No. 3, pp 1050-1053, May 1979 
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CONCEPTS 


*% = SAMPLE THE SIGNAL WAVEFORM ONCE PER BIT PERIOD 
*% DELIBERATELY ALLOW PULSES TO INTERFERE 
*% EXAMINE SEQUENCE OF SAMPLES FOR DECISION 


~ PRINCIPLE 


Match waveform against a library of all possible waveforms and select 
the one that gives the minimum sum of squares difference between 
the sample and the stored waveforms. 


FOUR oy 


a 
é 


89S-Ts- wal’ 


er 


H. Kobayashi 


Application of Probabilistic Decoding to 
Digital Magnetic Recording Systems 7 : 


Abstract: A digital magnetic recording system is viewed in this paper as a linear system that inherently includes a correlative level & 
encoder. This encoder can be regarded as a linear finite-state machine like a convolutional encoder. The maximum likelihood decoding = 
method recently devised by Viterbi to decode convolutional codes is then applied to digital magnetic recording systems. The decoding + 
algorithm and its implementation are discussed in detail. is 


Expressions for the decoding error probability are obtained and confirmed by computer simulations. It is shown that a significant © 
improvement in the performance with respect to other methods is achievable by the maximum likelihood decoding method. For =~ 
example, under the Gaussian noise assumption the proposed technique can reduce raw error rates in the 10-2 to 10-‘ range by a factor * 


of 50 to 300. These results indicate that the maximum likelihood decoding method gains as much as 2.5 dB in signal-to-noise ratio a 


over the conventional bit-by-bit detection method. 


1. Introduction 

In an earlier paper [1] it was shown that a digital magnetic 
recording channel can be viewed as a partial-response 
channel. The partial-response signalling or the correlative 
level coding is a technique recently developed by Lender 
[2], Kretzmer [3], van Gerwen [4] and by others in data 
communication systems, in which a controlled amount 
of intersymbol interference is intentionally introduced to 
improve the information rate [5]. In a digital magnetic 
recording system, on the other hand, the differentiation 
Operation inherent in the read-back process generates, 
in effect, a correlative level coded sequence. Since the 
representation of a digital magnetic recording channel 
in terms of its equivalent partial-response channel is 
essential to the development of the present paper, a brief 
review of some earlier results [1, 6, 7] is given. 

In the ordinary digital magnetic recording system 
saturation recording is performed, i.e., two stable states 
of magnetization represent binary data to be stored. Let 
{a,! represent an information sequence of **0" and “1” 
lo be recorded on the magnetic surface. The magnetiza- 
tion pattern a1) recorded by the NRZ (Non-Return-to- 
Zero) method is representable as 


mr) = D> (Qa — Wut — AT) — Mr, 
ket} 

where i1) is a rectangular pulse of duration 7 seconds: 
1, OS1<T 

“n= | at (2) 
0, _elsewhere, 


H. KOBAYASHI 


and 1(r) is a unit step function: 


2 | 
w= | ee (3): 
| 0, r< 0. 


Here the amplitude of m(r) is normalized by its saturation $ 
levels, i.e., +1 and —1 represent two saturation levels 3 
corresponding to “1” and ‘0’ of the sequence {a,}. $4 
We assume in Eg. (1) that m(1) = —1 for 1 < 0, Le.,: : 
the magnetic surface has been magnetized to the —1‘@ 
saturation level before the arrival of data stream fg 3 

In the read-back process the relationship between the x 
Output voltage e(1) and magnetization pattern :m(r) is 3 
given by 


‘v 


e(f) = E mo * A(t), (4) # 


where * means convolution and /(7) represents the mag- 
netic head field distribution characterized by the response y 
due to a unit step function in m(r). Figure 1 illustrates © 
waveforms at various stages in the NRZ recording method. : 3 

On substituting Eq. (1) into (4) we obtain j 


e(t) = A(t) * bP (Qa, — 1)(80 — kT) 
— 6&1 —kT— T)} + x | 
= 2 >> x,hle — kT), (5) 3 


ékaQ 
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Figure 1 Waveforms at various stages in the NRZ record- 
ing system. 


Qy ~~ Goi, k2>1 


Qo, ‘ k = Q. 


(6) 


As seen from Eg. (6), the sequence |x,} is a three-level 
sequence of —1°s, 0’s and +1’s. Unlike the situation in 
data communication systems, the sequence {x,} per se 

venerated nor clearly observed in any part of the 
ie_.  utg system. What we actually observe is e(r), a 
linear function of the sequence {x,} as shown in Eg. (5). 
In other words, we consider for analytical convenience 
thit the magnetic recording channel contains some 
imaginary correlative level encoder as a part of the system. 
I ja, takes on “1” and “0° equally likely and is inde- 
pendent from bit-to-bit, the sequence {.x,} possesses the 
following statistical properties: 


Prjx, = 0} = 4, Pr {x, = —1] = Pr {x, = +1} =1 
(7) 
und 
| ee ae 
Elva) = 9—3, [A-—- N= 1 (8) 
| 0, elsewhere. 


Equation (8) shows that adjacent digits are highly cor- 
related and hence {x,} is called a correlative level coded 
sequence. In other words {.x,} is a sequence that contains 
redundancy. 

Let e(7) be passed into a linear filter (1), the output of 
which is denoted by r(s): 


U1) = eft) * f(r). (9) 


-SSuanr 1971" 


a,ma,@ a,_, 


x, a,~ 4). 
( b) Two-state machine 


Figure 2 Discrete system representations of (a) the NRZ 
recording system and (b) the NRZI recording system. 


If the total response function g(7) = A(r) * f(t) satisfies 


the condition 


B(kT) = 6: .0 k = 0, +1, +2, Bg (10) 
then the sampled value of the filter output is 
MKT) = 2x,. (11) 


Equation (10) is satisfied if the filter f(r) includes an 
equalizer so that the effect of intersymbol interference is 
removed. However, the sampled voltage cannot be exactly 
equal to 2x, because of the noise and the residue of inter- 
symbol interference. Therefore, what we actually observe 
at a sampling instant is represented by the following 
random variable );,: 


MEM A, (12) 


where z, represents the total disturbance. 

From Eas. (6) and (12) we obtain the block diagram of 
Fig. 2(a), which is a linear discrete system representation 
of a magnetic recording system. Here {a,! is a sequence 
of “1” and *‘0”, {.x,} is a sequence of —1°‘s, 0°s and + 1's, 
whereas {),{ is a random sequence that may take on any 
real number. In earlier papers [1, 6, 7) we described a de- 
cision scheme that quantizes {),} into a three-level se- 
quence {q,}. The data sequence }a,} can be estimated on 
the basis of this “hard” decision output {g,! by solving Eq. 
(6). However, an erroneous decision in }q,} would result 
in the propagation of error in decoding the data sequence 
{a,|. To avoid such error propagation, |a,} is transformed 
into another binary sequence ja;} by the following 
relation before being passed into a correlative level 
encoder: 


ai = a,@ a/_, mod 2 


el 


<—— 


(13a) 
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Figure 3 (a) A representation of the NRZI recording sys- 
tem as a partial-response system; (b) error detection and 
decoding method in the NRZI recording system. 


where 


Gg = ‘ bas b,-1, k 2 1 (9) 


bo, k= 0. 


he response function /A(r) virtually satisfies the fol- 
, condition at sampling instants, 


h(nT) = ‘ n= 0 (10) 
0, n <Q, 

then the sampled value of the output voltage is 

e(nT) = 2c,, (11) 


oan 


which is a three-level sequence, as is clear from Eq. (9). 


_ Thus, the recording channel can be regarded as a partial- 
response channel with a discrete transfer function 


G(D) = 1— D. (12) 


From our earlier discussion on partial-response channels, 
the corresponding precoder here should perform a trans- 
formation that is “linear” in the binary sense, and has 
a transfer function [1/(1 — D)]sca2. More specifically, 
the precoder input {a,} is related to its output {5,} by 


& = a+ b,-) mod 2 


mod 2. / (13) 


usually assumed to be zero. 


JULY 1970 


wer't~ S x. Sd \" 


It is easy to see that the precoder obtained above by 
viewing the magnetic channel as a partial-response chan- 
nel performs exactly the function of an NRZI encoder 
(Fig. 3(a)]. That is, a symbol in the precoded binary 
sequence {b,} is always reversed from the preceding sym- 
bol when a OnE in the binary input sequence [a,} is to 
be recorded. This new observation, nevertheless, enables 
one to employ the general detection method developed 
for partial-response channels and to take full advantage 
of the inherent redundancy in the three-level channel 
Output. 

A possible structure of the error detection method based 
on the theorem in Section 2 is diagrammatically shown 
by Fig. 3(b). Here the rectifier or “‘mod 2 detector” in 
the conventional system is replaced by the combination 
of the inverse filter with transfer function (1 — D)”* and 
the decoder with transfer function [1 — Daca >. 

It should be pointed out that, in this special case, the 
error criterion is equivalent to the observation that voltage 
pulses read from a saturation recording system must alter- 
nate in polarity. However, it is reassuring to know that 
all detectable errors can be detected this way. The special 
form used to implement the error detection logic may 
also be of some practical interest. 


4. A new high-density recording scheme 

The function A(1) of Eq. (1) is the output voltage response 
to an impulse in d[sm(r1)]/dr; i.e., to a step change in the 
direction of recording medium saturation. It is known 
that a Gaussian characteristic, A(f) = a exp (—dt’), can 


’ generally provide a good fit to experimentally obtained 


characteristic voltage pulses.’ It is seen from Fig. 4(a) 
that although the pulse shape A(t) has no overshoot, it 
is not suitable for high-density recording, since a bit 
interval 7 of a fairly Jarge value must be chosen to avoid 
excessive intersymbol interference. 

The pulse in Fig. 4(c) is the same as that of Fig. 4(b). 
However, the sampling rate is increased by 50 percent; 
i.e., sampling is done at every 7’ second, where 7’ = 
£T. If we define the new time axis 7’ by 


n= t+ $7’ (14) 

and a function f(t’) by 

f(r’) = hr + 47), (15) 

then, as can be seen from Fig. 4(c) 

tor} Me (16) 
~ 0, n ¥ 0,1, 


where the values at 2 = 0, ] are normalized. Hence, the 
transfer function is characterized by 


F(D) = >> f(nT’)D" = 1+ D. 


aa Se 
A\ 


(17) 
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Figure 4 The principle of a higher-density recording system. 


Figure 5 (a) Interleaved NRZI recording method; (b) error 
detecucn aaa added to the ales 


ACD) nm BID) = 1p [cto 


" €cp) 


(3) 


R(O): reset pulse 
72 (b) 
. ROBSAYASHE AND VU. 1. LANG 
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Ourput sequence {e,} i) 1 0 =-1 1 1 0 


Figure 6 Waveforms at various stages in the Interleaved 
NRZI recording system. 


Then the transfer function of the whole svstem is, from 
Eqs. (12) and (17), 


G,(D) = GoD) F(D) = 1 — D*. (18) 


The corresponding precoder has a transfer function 
(1/(l — D*)Luas. The input sequence {a,} and output 
sequence {6,} are related by 


b, = a, by: «mod 2. (19) 


We shall call this new system the “Interleaved NRZI” 
recording system. The block diacram of the systern 1s 
shown in Fig. 5. Fig. 6 shows the actual waveforms 
observed act various stages in the Interleaved NRZI 
system. 

The Interleaved NRZI systern described above requitcs 
the channel response to sausfy Eq. (16). In practice. 
however, such desired channel response can not always 
be maintained (especially if the sampling rate ts pushed 
higher), and. as a result, systern performance may dete- 
riorare. A common remedy is the introduction of “shaping” 
or “equalization™ in the frequency characteristic of the 
channel. A possible frequency characteristic H(u) [the 
Fourier transform of the channel response function /1(1)] 
assumes a cosine shape in tht magnitude, as shown in 
Fig. 7, while possessing a linear phase characteristic. 
Equalization or channel shaping can also be realized in 
the time domain. Automatic or adaptive equalization 
developed in data transmission systems*® can thus be 
applied to the magnetic recording channel. However, this 
Subject is beyond the scope of the present paper. 
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Figure 3 State transition diagrams of (2) the NRZ record- 
ing system and (b) the NRZI recording system. 


Figure 4 Trellis picture representations of the state tran- 
sition in (a) the NRZ recording system and (b) the NRZI 
recording system. | | 


eo “oo “<0 “oO “oo 


Stare 0 <= () 


Stare | —— C) 


“Omi” : Inpus. } a,; 
+ 1.0% 1 : Ourpus 5 2,{ 


and a1, = 0. (13b) 


This transformation is usually called precoding in data 
communication systerns [2]. It has been shown [1] that the 
so-called NRZI (Non-Return-to-Zero-Inverse) recording 
method is equivalent to a precoding operation followed 
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by the NRZ recording method. Figure 2(6) is an equivalent 
discrete system represencation of the NRZI system. 

An algebraic method of error detection proposed in 
References 1, 6 and 7 makes full use of the inherent 
redundancy of the three-level sequence | x,]. This algebraic 
approach has been further extended to the ease in which the 
receiver makes a “soft” decision, i.e., the nun.ber of 
quanuzation levels is increased frem three to five or seven 
including ambiguity levels (6, 7]. 

The present paper describes a completely different 


approach to decoding the magnetic recording output. 


This decoding method is a very simple scherne to realize 
the maximum likelihood decoding (MLD) rule and is a 
probabilistic decoding scherne rather than the algebraic 
one discussed earlier. It will be shown that a significant 
improvement ir, the performance is obtainable by the 
proposed decoding scheme. 


2. Maximum likelihood decoding 

In 1967 Viterbi [8] devised a new nonsequennial decoding 
algorithm for convolutional codes. Forney [9] showed 
that this algorithm is in fact the maximum likelihood 
decoding rule. Omura [10] discussed the algorithm in a 
state-space context and showed its equivalence to the 
dynamic programming. 

A correlative level encoder can be viewed as a simple 
type of linear finite-state machine over the real number 
field as opposed to a Galois field over which a convolu- 
tional encoder is defined [11]. Now that we know the 


‘equivalence between a magnetic recording channel and a 


correlative level encoder, it isnot difficult to show that 
the Viterbi decoding rule is applicable to our problern. 

We define s,, the state of the imaginary correlative level 
encoder by the latest encoder input, i.e., 5 = a, im the 
NRZ recording system and 5, = a in the NRZI recording 
system. A precoder defined by Eq. (13) ts also a two-state 
machine; hence we can combine the precoder and cor- 
relative encoder in the state representation of the NRZI 
System. Figures 3{a) and (b) show the state transiucn 
diagrams of the NRZ and NRZI recording systerns, 
respectively, where 1 and 0 in small circles represent two 
possible states. Each time the machine receives a new bit. 
“1” or “O™, a state transition takes place depending on the 
input and the current state. Numbers — 1,0 or = 1 attached 
tO arrows represent the encoder output }..1. Although 
the diagram of Fig. 3 completely describes our system. 
the description in terms of the trellis picture introduced 
by Forney (9] will provide a better understanding of the 
decoding rule to be discussed. 

The trvllis picture of Fig. 4 shows the transition of the 
encoder state as a function of time ¢. Here the input “1” 
or “0” and the corresponding output —1. 0 or —I1 are 
attached to each branch connecting two states. Staring 
from s, = O the encoder follows a particular path according 
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Figure 7 Probability of error vs SNR: (A) the maximum 
likelihood decoding. (B) the bit-by-bit detection with thres- 
holds at ~A/2 and A/2. (C) the bit-by-bit detection with 
the optimum thresholds. 


the channel symbol x, takes on A, 0, and —4A with prob- 
abilities, 4, 4 and 2, respectively. The optimum thresholds 
Ure fep, and —/,,,, where f,,, is the solution of the equation 


SD(tep/o) = 24[(teve = A)/c}, ° (40) 
which yields 

A o” 
tao = 5 of. yi In 2. (41) 


Thus the bit error rate with these optimum threshold 
Values is given by 


a i! o(-4 _ 4 In 2) 


Hegetm 

Curves B and C in Fig. 7 are plots of Eqs. (39) and (41), 
respectively. Their difference is fairly small. We can see, 
however, a substantial difference between Curves A and 
C. For example, at SNR = 13 dB Pg;¢ ™ 1.1 X 10° 
whereas P. utp ™ 1.8 X 107°, ie., improvement by a 
factor of 70. The performance improvement is even higher 
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SNR (dB) 


Figure 8 Simulation results of the maximum likelihood de- 
coding method and the bit-by-bit detection method. 


for a higher SNR: the decrease in the error probability 
by a factor of several hundred is possible beyond SNR = 
14 dB. In terms of SNR, the maximum likelihood decoding 
‘method gains as much as 2.5 dB in the range of raw error 
rate 10°* to 107*. 

Now we shall report the computer simulation results 
and confirm the analytical results obtained above. The 
discrete channel models of Figs. 2(a) and (b) are assurned. 
The data sequence {a,} was generated through the random 
number generator program, and the noise sequence {:,} 
was generated by transforming the random variable with 
a uniform distribution through a polynomial approxima- 
tion formula of the mapping 7 ‘(-) [12]. It will be worth 
mentioning here that most existing subroutine programs 
under the name “Gaussian Random Generator” are not 
appropriate to this type of simulation, in which a high 
accuracy is required at the tail of the probability density 
distribution. 

The simulation results are plotted in Fig. 8, where 
decoding error rates for the NRZ and NRZI recording 
methods are marked by X and CO, respectively. The size 
N of data is 10° for SNR = 10 to 11.5 dB, and 10° for 
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(from laws of scaling are relentlessly driven towards zero spacing) 


On Tape: | 
Assume contact and set limit on number of passes 
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tpi 
— Media Noise « is 
J tpl 


or SNR decreases 3db as double tpt 


“sTtT=- — 


; | A.Hoagland 3/90 


Track Width 


@® (W = Track width) (s = signal output) 
© ca We ~7Alfci) 
1 


© tpi art] | 
® track density versus linear density . 


bits/in? « (fci)(tpi) « (fci) x s x Md (fei) 
holding s constant, maximize bits/in* in terms of bpi and tpi 


2 
d(bits/in’)_-g of fy=-b 
‘Cl 1d 
(ignoring noise) operate where a further increase in 
_ linear density causes a greater signal loss than an 


increase in track density would. 
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Read head pickup of old/adjacent track information 
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HEAD POSITIONING 


Two requirements on head positioning 
e Track following 
e Track seeking 


Track following Track seeking 
e Sources of track position variation e For minimum access time 
spindle bearing runout Use maximum power 
disk deflections max accel then max decel 
e Mechanical disturbances bang-bang servo | 
windage design involves choosing 


carriage bearing friction switch point 


vibration and shock 
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- The carriage is actuated by a linear (voice 
coil") motor. The carriage position is sensed 
from a magnetic servo pattern. Most commonly, 
the servo pattern is on a dedicated surface, 
(servo surface) and read by a servo head. The 

magnetic servo pattern is prerecorded during 


a drive manufacturing. 
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Voice Coil Actuator || 
Force and Acceleration ; 


F =ma (9.3) 
Ki = ma (9.4) 
Vio. | 
} = —— (9.5) 
R, +R, 
ae an, 
i eS eee CS 
M M R, +R, 
"BANG BANG" 


BEGIN TRACK 
FOLLOWING | 


| Source: A. Wagner 
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HEAD POSITIONING 


Equations - where mass alone considered 
: a=F/M 


a = acceleration, F' = Force and M = mass of actuator assem- 
bly. 
For voice coil actuator 


F=KkKI 


I = maximum current available from power supply. K = Bl. 
(B equal magnetic flux density and / = coil length). 


Kk 
a= Vi 
and | 
v= ball or t= cv 
| M KI 
v = actuator velocity and ¢ = time. 
Then 
, = Kit 1M 
M2 2KI 
Let X, = stroke length and 7, = seek time 
“e Xs « Alp or Il « lf 
M™** K T?, 
As Power needed « I” 
M\? X?2 
Power « €a T4, 
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Structure of a Control System 


Desired Output Control - System 
Output Error Signal Output 


Xo(s) E(s) F(s) Physical | (8) 
System 
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sid Sensor 
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Analog Compensators 
Advantages 


@ Cheap parts 
@e Control Theory is well Known 


Disadvantages 
® Difficult to modify 
® Components change with time 


® Circuits cannot adapt to changes 
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Output ce Asse sae j — Signal Output 
Xols) | BU Fa) Fo] Xe 
| ysica 
DSP men 
Foca. ue eee “Mee Ee __I 
X(s 
(s) _ Sensor 
isTtT= = 
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Advantages of.Digital Control 


® Greater Flexibility 
- Rapid prototyping 
— built in diagnostics 
- dynamic decision making 


® Filter Characteristics do change 
with time or temperature 


. ® Learning Ability 
- On-line identification 
- Adaptive control 


— IST 
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Data Track Servo Servo Tracks/ 
Allocation Head Servo Pattern 


“Tribit” Servo Pattern 


ey —__e 8 +) eee ° Cemented 


Se @ (Reamer Te eID 


bd ee ewwinns mame cman 


Track A 
Track B 
Track A 
Track B 


Pulse Sequence: 
1. Sync pulse 

2. Track A pulse 
3. Track 8 pulse 


Head closer to 


Track A: 


Track A pulse 
>Track B pulse 


Head in the middle 
between Track A and 
Track B: 


Head closer to 


4 TRACK aren | 
TK ( DISTURBANCE 
| EXCITATION in Pree 


HEAD 
: — 
4 
POSITION COMPENSATOR (AcTUATOR) pee 
IRAWS DUCER dAMPLIFIERS SYSTEM 


Tenok Foltowine SERVO 


as { 25: aan 
St ee iia (ala Core 
\ 
 ~ERROR 
‘ 
‘ 
4 
é \ as ysLtociry 
: ‘, ‘, 
§ 6 : 
| : , - aa . _ 
i \ 4 ~ 
| \ va 
o| \ ee. 
i e , 
i 7 TO la - am 7 
Re 
/ 


f ~ posiTioN 


j 
rae 
DV otnnes 


~~. 


Hb. HHK0H 


o> wee @ ae ae oom 


Bee (Re mee een nese ieee mee 
0.060650006 


0 o: U + GEIR ree gy RED EEE EP NE 0 AE eR aD me REE «AD aI TTD - Cte ae SRT = a ue ae 


CYLINDER 005 004 003 002 001 000 -00: 
. INVERTED: IN-PHASE = P QUADRATURE = ~ 
! 


| 
oe 


+ 


\ : 
IN-PEASE = P 


eae. - aap! eer ee ee ee ——— 
(P>Q)*(P>Q)——— : 


San 4 
004 003 002 O01 000 -OC: 


POSITION TRANSDUCER "IN-PHASE", “INVERTED IK-PEASE", 
es "QUADRATURE", AND "INVERTED QUADRATURE” SIGNALS VERSES 
ay RADIAL DISK POSITION OR CYLINDER LOCATION; LINEAR 
REGION SELECTION LOGIC AND CYLINDER COUNTING IN COARS= 
MODE; AND, SELECTED LINEAR REGION FOR FINE MOVE. 
F. J. SORDELLO 1968 


ENCODER AND WEDGE/TRACK 


emer ems - — “> pen . ee ran eect tenn od > ae a eee on pti — ee ee ee 


ow! Ole ane aoe sane 3 aoe 7 'e 


+ PROS; 

- SIMPLE SERVO WRITER 

- LOW OVERALL COST | 

- LOW OVERHEAD | 

= TRANSPARENT TO CONTROLLER - FREE FORM FORMATTING 

+ CONS; : | 

- LOW BANDWIDTH FEEDBACK FROM DISK (ONCE AROUND) 

- CORRECTS FOR LONG TERM EFFECTS ONLY 

- HIGHER INERTIA OF ARM STACK (ACCESS TIME CONSIDERATION) 
- LIMITED TP! UPPER LIMIT 


To 
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DEDICATED SURFACE 


a mm ~~ —_ oe ee oe Papen ee ttt Nt “Ane timer lly no ae pa ww ew ewan 


+ PROS; , 

- FEEDBACK FROM DISK STAC 

- WIDE SERVO BANDWIDTH | 

- CORRECTS RUNOUT 

- BETTER SHOCK AND VIBRATION IMMUNITY 

- WELL DOCUMENTED APPROACH 

- MORE PRECISE SERVOWRITING | 

- TRANSPARENT TO CONTROLLER - FREE FORM FORMATTING 
+ CONS; 

- FEEDBACK FROM SERVO HEAD ONLY 

- NOT EFFECTIVE FOR DIFFERENTIAL DISK SLIP 

- SERVO TO DATA SURFACE RELATIONSHIP DESIGNED IN (COST OF MECHANICS) 

- EXPENSIVE SERVOWRITERS 

- % OF OVERHEAD INCREASES WITH SMALLER FORM FACTORS 


1ST— 
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- FEEDBACK FROM DATA HEADS 
- REMOVES THERMAL TILT VARIATIONS 


- LOWER OVERHEAD 
+ CONS; | 
- "NON ZERO" HEAD SWITCH TIME 
- LOWER BANDWIDTH SERVO SYSTEM 
_- HARD SECTORED - LESS TRACK FORMAT FLEXIBILITY 
INCOMPATIBLE WITH DEVICE LEVEL INTERFACES 


SwrtTquNné 119 


ammlitt<s wn 
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OPTICAL SCHEMATIC 


AN T.R. LEO LIGHT SOURCE 1S USED TO 
REFLECT THE IMAGE OF THE SERVO TRACKS 


THROUGH AN ASPHERICAL COLLECTION “TC 
1H ONTO A QUAD ELEMENT PHOTODETECTOR. 
FLOPTICAL 
DISKETTE 
at] COLLECTION LENS 
com 
‘Me = = PHOTODETECTOR 
STRAY LIGHT BAFFLE 


ROUTING HIRROR 
0 SLIDER 
| Zan STOP RAW HEA 


SN\ys 
a sfstuvere ¢ “Aes 


‘“ PEC geeeega 


10ST —— 
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INTERPRE TER 


FLEXIBLE DISK ORIVES 


Brier Technology 


A 


KEEPING AMERICA 
COMPET ITIVE 


fields a first: 


a 20M-byte, 3'/2-inch flexible drive 


Mike Seither, Sersor Editor | 


San Jose start-up Brer Technology 
Inc. hopes to turn a few heads at the 
COMDEX/Fall show in Las Vegas 
next month when it introduces a 
20M-bste. 3¥s-inch flexible disk 
drive. That's more than double the 
storage of the newest breed of high- 
capacity S'a-inch flexible disk drives. 

With an average access time of 35 
msec. Brner’s BT 3020 high-capacity 
product also rivals the performance of 
comparable-size Winchester disks 
and offers a removable media to boot. 
This is something that is attractive to 
' OEMs and system integrators work- 
ing in the engineering. scientific and 
military segments. where work- 
stations are shared and the best way 
to protect one’s data is to remove it. 


HOW BRIER PACKS A BT 3020 DISK 
BRIER'S MAGNETICALLY BURIED SERVO 


EMBEDOED OR INTERLEAVED SERVO 


MINI-MICRO SYSTEMS /Octooes 1987 


Brier’s first product uses a proprie- 
tary system that embeds continuous 
servo information into the magnetic 
oxide coating of a standard 3'4-inch 
high-density disk. The company 
claims that it has been able to divide 
the disk’s oxide media into two log- 
cal lavers. Servo information in the 
deeper laver. recorded at a verv low 
freqency. is used to keep the read/ 
write head on track. The upper laver 
is dedicated solelv to user data and is 


written at a much higher frequency 


than the servo information. 

By companson. most vendors sell- 
ing 5'%einch disk drives in the range of 
4M bytes to 10M betes embed servo 
information directly on the data sur- 
face. As a result. up to half the toral 
data area is given up to servo. 

Thanks to this buried-servo svsterm, 


Brer has been able to overcome the 
biggest technical hurdle—track densi- 
tv. The BT 3020 has 777 tracks per 
inch. compared with the 96 tpi on a 
regular §'4-inch disk. By increasing 
the bit density and using advanced 


Embedded servos (ieft) typically 
uS@ Dart of the Cata surface to ine 
entire depth of the oxide media to 
Store read/write head ‘racxing inior- 
mauon. Brier honzontailly divides tne 
oxide into two logical layers (rignt). 
Prerecorded, continuous servo intor- 
mation 1s written af a low frequency 
onto the lower. Data is written over 
the entire upper layer ata gh fre- 
quency. Grier’s Crive elecironics 
sends data one way and feac ‘TaCK- 
ing information another. 


AN FS SO OD OF ED aeRO 


MYLAR SUBSTRATE 


SOURCE. BARA TEQORLIGY OS, 


29 


Write 
element 
>Write gap 
i 
Read 
element 


===) Read gap 


2 
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7 -— ew pews 
Kes 


OUTPUT VOLTAGE 


TFI Head 


POSITION 


Read Element 


OUTPUT VOLTAGE 


MR Head 


POSITION 


HisTt= = 
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DISCRETE TRACKS - POSSIBILITIES FOR HIGH-DENSITY MAGNETIC pisxs? 


Ian L. Sanders and Steven E. Lambert 


PROCEEDINGS 


VLSI 

and 

. Computer Continuous 
Peripherals 


3° Annual 
European 
Computer Conference 


Magnetic 
coating 


Substrate 


Discrete 


Magnetic 
coating 


Substrate 


Underlayer 


0.10 3000 fc/mm 
a 
® 
eee 1000 fc/mm 
oe" 


0.06 500 fe/mm 
> 


RMS Noise 


0.0 1.0 2.0 3.0 
VTrack Width (pum)1/2 
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Media Noise (uVrms) 


Edge Noise Reduction in Discrete Tracks 


45 
a Li 
35 ‘ On-track 
30 
25 | 
a. 
20 / 
Written edge 

15 © .—— 8 

a — 
106 

@ 
5 

; Patterned edge 

A A Aa A 
7 1000 2000 3000 


Written Density (fr/mm) 


# 
tee 


LITHOGRAPHICALLY DEFINED SECTOR-SERVO 


/ : ae a i Soe SS . 
Mee LD A oi. 
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Track Density in Magnetic Storage 


Buried Special Medium. 
& Transducer 
ce Intermittant Control, Real Estate 
Escicarsonenavbacrsdecnte Geometric Stability of Head-Disk Assembly 
1000 2000 TP} 3000 4000 


lusv 


MSS '90, Monterey 


DATA STORAGE DEVICES 


MIST 
A.S. Hoagland '90 


Error Control and Correciion 


Raw BER — 


Bit Density ->> 


e For Density > 10° bpi¢ 
Rew BER 10“-10° 
=> ECC Overheed 10-20% 


(1s 'T—_ 
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The error probability P, e@qual to the BER, is also called 
the transition probability, where transition refers to a “i" 
changing into a “O", and vice versa. The probability of no 
errors is i-p, and n independent symbols are therefore received 


error free with a probability of (1i-—p)". We can summarize: 


Ptcorrect) = Clep)" 
P(l error) = pl(i-p)7r7 
P(2 errors) = p#(1-p)"°,2 
Srmtestntentctemeatenttee ape nfeenenfanirabe fen ena eed 
P(b errors) = p®(i-p)"-») (25.6) 
pee epee ep ethemechemes spe am sfpmeprndpmsnmogp> when penehtpn ene cepomarmn ep ep exp 


The actual probability will be greater, since one (or more) 


errors can result in different error patterns: 


ONE error: EX X X XxX 
= or X EX X X 
= “or xX EX X 
= or x Xx E X 
= oF x Xx XE. 


For a n-dimensional word there are n different patterns for 


asingle error, and the true error probability is therefore 
P(l error) = n-p(il—p)-~3 


Two errors will result in n-(n-1)/2 patterns: 


EE X x Xx EXE X 
EX EX X Xx EX X 

Ex XE xX Xx XE x 
Ex XXeE XX EXE 
X EE X X X X XEE 


lIST 
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Retransmission 


| request 


Error al | 
avoidance ee ear 
is oetection correction 
redundancy 
encoding 
Channe! 
errors 


Figure 7.1] The maior processes in an error-handling system. 


Error 


concealment 


hISsT 
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133 


RECORDING OCRSITY 
134 HARD YS. SOFT SECTORS 


Recording density (RO) fs the number of user bits which can be recorded 


| Hard Soft 
r unit length of track (e.g., inch). For our purposes ft 1s most conventent 
(factory sectored, (user sectored, 
take this unit length to be Tagn. the minimum recorded spot size. (This servo tracking, open-loop tracking, 
Medfa high density) low density) 
lantity {5 @ measure of the R/W channel bandwidth and fs usually difficult to Type src ty 
crease.) That fs, flexible High-performance floppfes 
(mylar) 


floppfes (e.g., 
Bernoul If Box! 


RD = number of user bits In interval Tay, 
: old 8" flopples 


« (m/n) > number of code bits In interval Tain 
rigid 


High-performance PC Winchesters 
it the Interval Taq (5 divided into d¢i detection windows of length T,; In (aluminum) 


Winchesters 


her words, there are d#} code bits fn an interval of length Tay,. Thus, 


Hard-sectored disks are formatted very precisely at the factory. They 
are capable of accurate tracking and thus achfeve high recording density. 


This high density, in turn, makes system more sensitive to nofse, IS! and 
timing jitter. 


ERRORS ON MAGNETIC DISKS 


° Surfaces of both rigid and flexible disks have few defects and 
very little dirt. Also, sfgneal-to-nolse ratio {s typically very 


high (e.9., 35db). ERROR CORRECTION METHODS FOR MAGNETIC DISKS 
© Ofsks are formatted and certified before use. During certification ' 
a test pattern 1s written at each sector and then read back. If it - Detect-and-Reread 
contains any errors, the sector 1s discarded. (Several spare 
sectors are provided. ) . Error detection code only. Reread until! correct read achieved. 
(System parameters like tracking offset or demodulator tim! ffset 
¢ As a result of the sbove two rhe ad rrors ay rare. Typical bit 9 ro timing offset may 
error rates gre in the range of 10°19 eo 10-! » giving sector error be adjusted on different reads.) Code must have large number of check 
rates of 10°’ to 10-9. bits Co achleve moderate relfability. Disk performance poor. 
© f€rrors are caused by: 
- dirt * mis-tracking 2. Reread-and-Correct 
> head skipping * demodulator timing jitter Re ee es 
° ran lta oxide areas * intersymbol interference Fire or Computer-generated (single-burst correcting) code. Same as 
* no 
Detect-and-Reread, except that after a fixed number of unsuccessful reads, 
© A typical error event Is caused by misreading s single flux correction 1s performed. Correction fs typical! f h 
change. There are three ways in which this can occur: P ypicatly performed by CPU wit 
- Drop tn - reading a flux change where there {s none. location and pattern provided by controller. Better performance. 


* Drop out - missing a flux change. 

> Shift - reading a flux change In the wrong timing Interval. : Weavcect-sddenecead 
. eee ee nerese 

© Errors are extended by the modulation code decoder. Result — 


typical error event 1s a short burst of errors (2-10 bits). Error correction ts performed tn controller. If a detectable but 


2 uncorrectable pattern occurs, reread ¢s Invoked. Typical code: 
e Most errors are soft fn that they do not repeat on multiple reads. JSC A 4 ¢ yp ode 


Interleaved RS2. Since rereading {$ rare, best performance. 


Binary System: 


| Module 2 
Addition Multiplication 
O 1 | O 1 
0; O 1 0; 00 
1/10 | 1; 01 
Vector 
transmitted word ( t ) 1011001 
received word ( T) 1001001 
error pattern ( e ) 0010000 


t=T+e 


thus once error vector determined can add to received word to correct 
and get transmitted word. 
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Parity Check 


Simplest method of error detection. 
Coding Rule: 
Add additional bit (named parity bit), so that 


the entire message including the parity bit 
contains an even number of 1s. 


-e.g.: odd parity message: 1101 f 
even parity message: 1111 © 


Parity bit 
Hamming Distance (Definition): 
Number of bits that differ between code words 


e.g.: Codeword 1: 1101 
Codeword 2: 1111 


The two messages differ in only one bit, the 
Hamming distance iS one. 


The addition of a parity bit increases the 
Hamming distance by one. 


A Hamming distance of two allows’ the 
deteclion of one error and no correction. 


‘1s TI———__ 
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ROW PARITY BITS 


4 4 


PARITY COLUMN PARITY BITS 
ERRORS 


eS —___ 
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k n-k 


A code word o: length 1 consists of an information block of 
&k symbols and a parity block of n—-# symbols; each symbol] com- 
prises 5 bits.. The number of possible words of 1m symbols is 2”°. The 
parity bits are fixed for each combination of the As information bits 
in accordance with established encoding rules. The number of code 
words is thus 2“. It follows that the fraction 2'*-")* of the number 
of possible words consists of code words. 


(1s TI— 
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CONSTRUCTION OF DECODING TABLE 
- BASIC CONCEPTS 
noe ee In a decoding table, al! 2" possible recelved words are grouped Into 2k 
columns, with each column headed by @ code word. important point: entries tn 
a column are selected to be the words most Ifkely to be recelved when the 


number of Information bits (k) - transmission of a code word conveys a corresponding code word 1s transmitted. 
EBT message, where k <n. 


code length (n) - number of bits per word (all words have same length) 


number of parity check bits - n-k Example: n*5, ke2 

code rate or efficiency (R) - k/n. Since k <n, it follows that R< J. Message 1 2 3 4 

systematic code - block code In which the first k bits of each code Code Words 11000 OO1T0 TOOTT 01101 
word are identical to the k-bit message 11001 oOoOt1! 10010 O1100 

encoding - associating a code word with a k-bit message vies : , : , : : : : ; ; : : ; : 

decoding - associating a k-bit message with a received word Sa 10000 o1110 t101! 00101 


01000 10110 o000TT F1101 


—_—— ee ee —<—_ a a —_—_——_ «a —_ a aad am sien 


tTrTI1O0 00000 1OTOT OFT0T? 
01010 10100 0000171 111171 


ENCODING TABLE 


f d thef 
URES NSS SMBS: NNSL ON MESS Ages 2 ne: TRENT Note: Since there are 2" entries in the table and 2* columns, there must 


associated code words. In principle encoding can be be 2nyak s gn-k rows. 
accomplished by looking up message in table, then select- 
ing corresponding code wor d. 
USE_OF DECODING TABLE 


Example In a system with four messages and binary 
Case 1: €rror Correction Only; Recetved word {s located In the table and 

is decoded into the code wrd at the head of its column. The decoder 
tsble is es shown, outputs the corresponding message. 


code words of length 5,one possible encoding 


Case 2: Error Detection Only: The decoder determines {f recelved word {s 
above or below the solid Vine. If (t 1s above, then it {s a code 
word and the corresponding message fs outputted. If it Is below, an 
error has been detected. 


Case 3: Combined Error Detection and Correction: The decoder determines If 
received word Is above or below the dashed Tine. If tt {'s above then 
{t 1s decoded Into the code word at the head of the colum. If ft 
is below an error has been detected. 


REAREST REIGHBOR DECODING 13 14 


ERROR CORRECTION CAPABILITIES ARE DETERMINED BY din 


Errors are normally unlikely; for errors which occur at random, the pepe Se ee 


ee eee ae ee 


eases as the weight Increases. . 
ibability of a multiple error pattern decreases | g renoe Curcec tion: 


other words, few errors are more likely than many. | Case I: bth = ?t ¢ > waranteed error-correction capability of code = t 
A good error-correcting code fs one which successfully corrects Case 2: dain = 2t 4 2 
‘terns having few errors and (f.e., Iikely ones) and (necessarily) fails Proof: let qt transmitted code word 
- high-wetght patterns. ie neighboring code word distance din away 
v = received word at distance e < t frome 
The decoding rule which {s most appropriate for random errors fs 7 t 
by triangle inequality 
irest_ neighbor decoding. Here the received word {s decoded into the code d(c, sc )< d(c, ov) ed(v.c ) 
n 7 n 
rd which {s nearest to ft ({n the sense of Hamming distance). Thus, the $O 
d(v,c > d(c..c ) - d(c,.v) 
coder always guesses that few errors, rather than many occurred. n) ~ ; t’ on ( t 
> Onin ; 
>2t+ed - tzted 
HAMMING WEIGHT AND DISTANCE Hence Ch is farther away and so nearest-neighbor 


decoding corrects any error pattern of weight < t. 
ep (Hawmeing) weight of a vector v fs 


w(v) = number of non-zero components tn v 


NEAREST NETGHBOR DECODING AND dain 


t (Hamming) distance between two vectors v, and v, 


d(v,, v>) * number of places in which v, and PI differ 


C, = transmitted code word 


t 
C., « code word at distance dain”! 


bed w(v, - v,)° 


e w(v, ® v,)* inthe binary case 
v = received vector 


e minimum (Hamming) distance ef a code fs the distance between the two pe ee: _ ~ 
so ANY ‘\X < 
Closest code words. / \ / | \. | 


: \ | \ 
a | 


Hasming distance ts primarily useful In dealing with \ C 24 n 
random errors. \ - or . / 


4 


at {s, Gain” smallest value of d(v,. vy) where 1 # j. 


ae “a NX we 


a ee 


*Addition (or subtraction) of vectors {s always performed by adding 
wr subtracting) corresponding terms. 


A code with minimum distance 
anin 
can correct any error pattern of 
d. ] 


Inin 
2 


or fewer errors in a codeword. 


11s T— 
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ERROR CORRECTION AND DETECTION 


Any error pattern having < e. errors can be corrected and any pattern having 


< ey errors (e, > e.) can be detected provided: 


ct eee - | 


Proof: Since ze. <« d, error correction works. For detection to work the 


recelved word must be further than e. from every other code word. 


d(v.c,) > d(c,.c,) - d(c, ov) 
>d -(d-1- e.) 


ea 
Thus, detection works too. 


A SPECIAL CASE 


if e,. * 0 we have the result: 


A code with minimum distance d has 
error detecting capability ey" d-1. 


Since error correction capability {fs eo. each code word {s centered 
in a sphere of radius Qe. These spheres are at least distance - d apart. 


correction , quaranteed detection false correction 
region ' region region 
' ' 
Gi ey 


If fewer than d- e. errors occur, the received word cannot “react! the 
sphere around C,- Hence any error patterns of weight d- = lor less will 
be detected so the error-detecting ability of the code fs 


ey" d-e.- 1 


CORRECTION VS. DETECTION 


At design time the designer decides how much distance wil! 
be used for error correction and how much for detection. Once the 


decoder {s bullt, this dectsfon {s frozen. 


Example: Suppose d=©7; then the 


following chofces are ee *4 
possible. —_——— 
0 6 
1 5 
2 q 
3 3 


CORRECTION OF ERRORS AND ERASURES 


Sometimes the read-head output signal is quantized into 
three rather than just two regfons. The third output {s then 
seen by the decoder as an erasure. 

We treat erasures as having Hamming weight 1/2. It takes 
two erasures (but only one error) to confuse a transmitted 1 
with a 0. 


one error two erasures 
7 a a 
de1/2 

d=1 x =) 


| de1/2 
cocci |, 0 _t 


Correction of efther errors or erasures {ts normally done 
by finding the code word which has the smallest distance to 
the received word; that is, nearest neighbor decoding. 


ERROR CORRECTION 


Hamming code (single error correction). 


_ Need check word that can point to single bit in error. 
n = word length 
kK = # information bits 

n-kK=# check bits 


n-k parity check bits have to indicate n + 1 conditions: 
l.e.; nN possible error locations 
+ error free condition 


Thus: 2 (tk) =n4+1 


or; =n=2 (m4) -4 


D-k ia k 
2 3 1 
3 7 4 
4 15 11 


lis Ti __ 
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TABLE 3.4-2 Encoding a 4-Bit Message and Locating Error 


Encoae 
] 


Q 


0 


th) 


i 


E 
l 
i 
x 


0 


Locate error 


Check ] 
Check 2 
Check 3 


Syndrome 


Correct (add 1 into the error p 


0 


] 


l 
] 


4 5 6 
“04 
00 3 
00 3 
13 5 
0°0 0 
2 3 6 
1 Q ] 
45 6 
001 
011 


0. O° J 


~) 


>— 


— ~} 


—_ ~)} ee ~]) 


1 


Position 
Message 


Encode 


Error 


Receive 


Fails —— 1 


Fails —— 1 


,-COnect ===>0 


3 —— position of error 


osition) 
Correct error 
Corrected message 


lis T—— 
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Code Sample after R. Hamming: 


Error correction code for four binary digits. We 
will need three parity checks, i.e. the message 
including parity bits is seven bits long. If there 
is no error, the syndrome will be 000, if there 
is an error, the syndrome gives position of 
error. 


Example message: 1101.For simplicity reasons, 
the parity bits are interlaced with the data. 
Dashes symbolize positions, where parity bits 
are to go. 


Bit Position 7 6 
Data 1 1 


© dl 
aS 
— WN 
NO 


18¢ parity check: 

Parity at position 

1,3,0,7. 

Parity bit position 1:°0" 1 10 - 1 - @ 


2°4 parity check: 

Parity at position 

2,5,6,/. 

Parity bit position 2:°1" 1 10- 1 17 0 


3° parity check: 

Parity at position 

4,9,6,7. 

Parity bit position 4:°"0" 11060 1 1 0 


Parity Check: 


Bit Position 7 65 4 3 2 1 

Message 11001 1 0 

Syndrome: 0 0 0 ~ no error 
Error: Bit 4 changed into 1 

Message 1101 1 7 0 

Syndrome: 1 O O- 4 decimal: 


Position 4 is position of error. To correct 
message, reverse bit 4. 


Parity check of corrected message: 


Message 11001 1 0 


Syndrome 0 0 0 

The syndrome of the correct message 

is OOO ~ message correct. 

Note: The Hamming error correction 


code has been’ chosen to 
demonstrate the principle of 
error correction. 


More sophisticated codes are used 
in disk drives. 


THE_REED- SOLOMON CODES 


These are cyclic codes whose symbols are binary m-tuples 
rather than bits. 


length = N > 2” 1 symbols or n * m(2™ - 1) bits 
1f use NW - K check symbols (m(N - K) check bits), then 


code can correct t = ((N - K)/2]} randomly located symbol 
errors. Theat is, d= @-K¢ }. ° 


For t-error-correction, need 2t check symbols so N-K = 2t 


nant tN 


In a t-error-correcting RS code, the following error 
patterns are correctable: 

1 burst of total length b, = (t-1)m *+ 1 bits 

2 bursts of total length bo * (t- 3)m * 3 bits 


{ bursts of total length b,* (t-24+ 1 3)m 4 21 - 1 bits 


Decoding 1s only slightly more complex than for BCH codes. 
In addition to the locations of the errors, their values 
(m-tuples) must also be determined. 


EXAMPLE 


Form» 8 we can construct the following RS codes 


(amrng atiers) 


(2040, 1976) 
(2040, 1912) 


(2040, 1784) 107 
where b, {s the total length of { correctable bursts. 


Mote: n, k, and b, are expressed as numbers of bits. 


THE_GENERATOR POLYNOMIAL OF AN_RS CODE 


eee oceeainecnmeeneneaenrennenmea meneemenat enmcetaaatibamnm ne casero nnataneiaeee 


For an RS code of length 2 - 1, symbols are binary 
m-tuples, f.e., elements of Gr(2"). The generator 
polynomial of an RS code with minimm distance d has the 


form 


gtx) © (xt alle ta) (neatly 


where Pa {s an element of GF(2"). 


Actually, any set of d-1! consecutive powers of a can *e used. 


Thus, we could also write 
q(x) . (x ry a iG r a) “ox ry ad?) 


[An RS (WK) code has d = WET] 


EXAMPLE OF AN RS CODE 


To construct a J-error-correcting RS code of length 
15 with 4-bit symbols, a total of d - 1 © 6 check symbols ere 
necessary. The generator polynomial of this (15,9) code fs 


a(x) = (x Halle a N(x eax ean Had(x & af) 


- x? ‘ o>, aa ‘ te P ane P an nl 0° 


Viewed as a binary code, this {s 8 (60,36) code which 
is capable of correcting any error bursts which are confined 
to three 4-bit symbols. 

As in the binary case, encoding can be accomplished by 
dividing by g(x). 
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codewords 


Figure 7.8 The vertical columns of this diagram are all codewords generated by the matrix of 
Figure 7.6, which can correct a single-bit error. If these words are recorded in the order shown, 
a burst error of up to four bits will result in one single-bit error in each codeword, which Is 
correctable. Interleave requires memory, and causes delay. Deinterleave requires the same. 


iis Tr—— 


A. Hoagland 3/390 


Coceworc A 


) Coceword E 


Sssition 
on iape 


Codeworc A 


PRUE TITTIES TITLED} | eoereree 


interleave 


sete = AN =A ZZ =A ZA 
2 B2 C2 Al B3 B BS Cs 


Ai B1 Ci A C3 AG 4 C4 As 


Burst Random 
error error 


(b) 


Figure 7.27 At (a). interleave controis the size of Dursi errors in individual code words. but at 
(b) the system falls down when a random error occurs adjucent to a burst. : 
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— 176 CODING GAIN 


To achleve a destred bit error probability (e.g. p 19°) 


Ee memimeadiiicnmemantanaditediiiend 


+ —— sige PLOT OF GECODED BIT- ERROR an uncoded system requires a certain value of EM (e.9..11.3 db). 
| PROBABILITY - VS CHANNEL 
we —— A coded system can achleve the same value of error prob- 


BIT CRROR PROBABILITY 


ability at a lower value of ELM: The difference between these 


two values of EMM is referred to as the gain of the code. 


Example: To achteve a bit error probability of 10"7, the 4-error 


correcting (255,223) code requires a channel! bit error rate of 


- 3.06 


p:- 10 . This value of p,in turn requires a signaleto-nvise 


ratio of E/N o” 7.5 db. Hence 


Coding Gain = 11.3 - 7.5 = 3.Adb 


CODING GAIN OF VARIOUS CODES 
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"Convolutional codes are self-orthogonal majority - logic decodable 
codes. 
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ECC IN DISK STORAGE DEVICES 


| DENSITY | DATA 
ipevice| pp: tpi | CODE ele 


| 250 | 100 20: HEZ | PARITY CHEER 


YEAR 


CHIP CAPABILITY 


‘yose | 23x | 520 50 | NRZ [CRC CHECK 
| 1020 +#4«©50 MULTIPLE ERROR DETECTION 
| 23xN | 2200 100 | FM 

| 
‘ascx | 38x | 4000- 200-) SINGLE BURST CORRECTION 
15000 800 | 

| 

| MODIFIED FIRE CODE 

_ |CORRECT./DETECT. = OF CHECK BYTES 
“ 187) MPM {11 BITS / 22 BITS 7 / RECORD 

:e73 | 3340 | 5636 300! MFM | 3 BITS / 21 BITS 6 / RECORD 
1s76 | 3250 |6425 475 | MFM | 4 BITS / 10 BITS 6 / RECORD 


MODIFIED REE D=SOLONMNON CODE 


9 BITS / 17 BITS 9 / BLOCk 
9 BITS / 17 BITS 12 / RECORD 
i? BITS ./ 33: 817s 12 / RECORD 
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ECC IN TAPE STORAGE PRODUCTS 


YEAR: | : 

EHIP, DEVICE DENSITY | ECC CAPABILITY | 

Iyo5N) T2X 1100—556 bp: IPARITY TRACK & LRC 

| | | | 

: { i 

| — eee = 

) | 

lyos2i 72S IB8OC bp: NEC (penile Seek CRE 

! IS TRACKS /i/2"  iONE TRACK CORRECTION ON REREA 

| | : 

i | | : | 
{ i 

1886 2200 11600 ‘spi PE {PARITY TRACK & PE ENCODING (CRC) 


9 TRACKS /1/2° i“ TRACK CORRECTION ON}TRE<FLY 


2878.) See0 16250 bpi CCR IGCR GROUP CODED RECORDING (ORC) 
I9 TRACKS/is2. {PARITY TRACK & 1 CHECK BYTE/7 BYTES 
| | TWO TRACK. CORRECTION ON=THRE< FLY 


4975 | 3850 (6888 bpi ZN INTERLEAVED SUBFIELD CODE 

MES |67 tpi WORD * 35 SECTIONS. SECTION. = 16 23725 
REDUNDANCY 2 SECTIONS/15 SECTIONS 
ON-THE-FLY CORRECTION OF 2 SECTIONS 


119BN | 3480 |-22000 bpi (0.3) IMULTI-TRACK CORRECTION 
| 18 TRACKS/i72z: |ADAPTIVE CROSS PARITY CODE 


MIST 


A. Hoagland 3/90 


DISK ARRAYS 
PARITY DISKS 


MIRRORING 
N+ 


2N 
(Also called shadowing or dual copy) Byte interleaved: 


@ Array viewed as single 


@ Complete duplication of data, 
logical disk 


i.e., disks paired 


e@ Parallel readout for high 
data rate - striping 
(Synchronized spindles) 


e Read from disk with shortest 
access (have twice as many paths 
to data) 


&® fwice the number of writes Block interleaved: 
@ Second write to update parity 
(therefore rotate parity among 
disks to minimize delays) 


Used in transaction processing 
where fault tolerance essential 


Notes: Disk drives indicate when they become inoperative 
Hot spares used for availability and minimum MTTDL 


MtIsT 
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Improving Reliability with Redundancy 


© Add redundant drives to handle failures 
Redundant 
Array of 
Inexpensive 
Disks 


e Redundancy offers 2 advantages: 
1) Data not lost: Reconstruct data onto new disk 
2) Continuous operation in presence of failure 


@ Several RAID organizations 
Mirroring / Shadowing 
ECC 
Parity 


Source: Industrial Liaison Program, Berkeley, CA - Mar.'90 
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Types of Disk Failures Disk Lifetime Distribution 


: 

Transient | NS : | 
=» retry or ECC modification = re 
Media a . = exponential lifetime good first guess 
=> ECC modification, factory remapping, => desparately seeking better data 

read after write verificaton . 

® also assume independent, but ... 

Electrical support hardware 
=> end-to-end checksum, complex interface bad batches 

protocol, internal buffer & path ECC human operators, maintenance 
Mechanical 


=> hope you can scavenge data 


Effects of Arrays on Failures 


Lifetime Parameters 


® higher data rates 
=> higher misdetection/miscorrection rates 


® temperature => better error detection 
= age = striped files 
| => more files exposed to each failure 
= power cycles > encourages failure tolerance 
® seek events ® more disks 
= more frequent failures 
® read/write events => automat isolation and easy repair 


® larger total capacity 
= more to backup 
= reduce frequency of backups 
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TYPES OF MAGNETIC RECORDING FORMATS 


LONGITUDINAL 


HELICAL 


e¢ 


ANMIPEX 


HELICAL SCAN HEAD ASSEMBLY 


-HEAD 


PRE-LOAD 
SPRING 


STATIONARY 
DRUM 


BEARING 


COUNTER. 
WEIGHT 


STATIONARY 
TRANSFORMER 
LEADS 


FERRITE ROTATING 


TRANSFORMER “\ 4 
LEADS — A 
| - 
STATIONARY AZ 

iN WA Ce 
TRANSFORMER ) Vesey 
LEADS gp 


VA ™ INTERCONNECT 
PRINTED 


QARITOIAIC? MCAT 


Read After Write 


2nd Stripe 1st Stripe 
2nd 360° Rolation 1st 360° Rotation 


Gay, Ge | 
nets 1st 180° of Rotation - Write 
mee 


2nd 180° of Rotalion - Read Back 
Check 


1st 180° of Rotation - Write 


2nd 180° of Rotation - Read Back 
Check 


WRITE HEAD =~ TAPE 
1ST 180° OF ROTATION 
1800 
SERVO HEAD | 
APM DRUM 
READ BACK CHECK HEAD (_)=-CAPSTAN 


2ND 180° OF ROTATION 


tl. 


Read After Write 


Normal Sequencing 


WLR, 
SON 


Rewrile Tew: 


MCR HEAD/STAPE FORMAT 


MAGNETORESTIS 


DATA READ 


TAPE 


STEPPER NOTOR PREAMP 


AND —.- 
BAND ASSENBLY Wi 


CARRAICE VOICE COIL | 
CUIDE RODS \ MOTOR 
OPTICAL DETECTOR 


MAGNETIC HEAD 
_ PREAMP 


PHOTODETECTOR --—~---—— 


LENS AND FPS Dg ya ne 
ROUTING MIRROR ta 
zomp . {he pane 
Fluptical == 
Diskette are LED 
3.5 20 MB 
FLOPTICAL _— 


DISKETTE DISKETTE 
SERVO TRACKS 
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OPTICAL RECORDING 
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ST "Data Storage Tutorial" 
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Optical Data Storage 


TABLE OF CONTENTS 


COMPARISON OF OPTICAL AND MAGNETIC RECORDING 
OPTICAL COMPONENTS 
- laser diodes 
- lenses, beam-shaping prisms, beam splitters, 
holographic lenses 
AUTOMATIC TRACKING AND FOCUSSING OF THE BEAM 
READ-ONLY DISKS: CD, CD-ROM 
RECORDABLE OPTICAL DISKS 
WRITE-ONCE DISKS 
REVERSIBLE DISKS: = phase-change 
- dye-polymer 
- thermo-magneto-optic (TMO) 
CURIE-POINT AND COMPENSATION-POINT RECORDING 
ORIGIN OF MAGNETO-OPTIC EFFECTS 
ENHANCEMENT OF KERR MAGNETO-OPTIC ROTATION 
PROPERTIES OF BARE EARTH - TRANSITION METAL AMORPHOUS FILMS 
SPUTTERING OF TMO FILMS 
DESIRABLE PROPERTIES OF TMO FILMS 
FACTORS LIMITING THE SIZE OF WRITTEN BITS 
COMMERCIAL TMO DISKS AND DRIVES 
FLEXIBLE OPTICAL DISKS AND TAPES 
ISSUES 


MAGNETIC RECORDING 


SIMPLE —- NO CRITICAL COMBINATION 
OF INGREDIENTS NEEDED 


RELIABLE 
INEXPENSIVE MEDIA 


NO PROCESSING NEEDED 


CREMOVABLE ) 


INFORMATION STORED IN A STABLE 
STATE OF THE MEDIUM . 


UPDATABLE 


INFINITELY ERASABLE AND REVERSIBLE 
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Areal Density versus Time 
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Device 


IBM 350 
IBM 1301 
IBM 2311 
IBM 2314 
IBM 3330 
IBM 3340 
IBM 3350 
‘STC 8650 
IBM 3370 
IBM 3380 
NTT PATTY" 
NTT PATTY" 
IBM 3380E 


NTT GEMMY 


IBM 3380(K) 


RIGIv DISKS 


Bits/in 


_ Tracks/in 

100 20 
500 50 
1100 100 
2200 100 
4040 192 
5600 300 
6425 476 

~— 6425 952 
12134 635 
~ 15000 801 
13970 1092 
25400 1800 
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Bits/in2 
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Thickness 


Spacing Gap 
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OPTICAL 


M@ Large Separation 


| MAGNETIC 


M@ Very Small Head-to-Disk Separation 
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Focussing 
— (Dy + Dg) — (Do + Dg) = 0 


Tracking 


Signal = 
D,; + Do + Dg + Dg 
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OPTICAL vs. MAGNETIC RECORDING 


ADDITIONAL DIFFERENCES 
. MAGNETIC RECORDINGS ARE ALWAYS "SWITCHED ON” 


(OPTICAL RECORDINGS - INACTIVE UNTIL SWITCHED ON; LIKE DNA 


- ONLY THE ILLUMINATED BIT CONTRIBUTES TO THE SIGNAL) 


» WHOLE RECORDED AREA OF MEDIUM EMITS FLUX THAT CAN BE PICKED UP BY THE READING HEAD. 


- MAGNETIC HEAD READS SIDEWAYS 2° SO ADSACENT 
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AS WELL AS DOWNWARDS 
- SENSITIVE AREA NOT LIMITED TO GAP 
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OPTICAL BITS ARE EQUIAXED 


BIT LENGTH 
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TRACK DENSITIES OF OPTICAL AND MAGNETIC DISKS 


VERSUS HEAD-TO-DISK SPACING 
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Focal distance 8 mm 
| ; NA 0.45 
_ Field diameter 1.0 mm (OPD at 400 wm from the centre is 
: 0.052 d) 

Field curvature . 6pm at the rim 

| Diameter 7.5 mm | 
Length 11 mm : 
Weight <2 g (mount included) | 


Focal distance § mm 
NA 0.45 
Field diameter 2.0 mm (OPD at 800 um from the centre is 
0.047 d) 
Field curvature 24 pm at the rim 
Asphericity (S1) +25 um (with respect to the best fit sphere) 
Asphericity (S2) +2 xm 
Diameter 7.5 mm 
Thickness 8.5 mm 
Weight <2 g (mount included) 
Focal distance 3.7 mm 
NA (object side) 0.10 
3 NA (image side) 0.45 
Field diameter 0.2 mm (opp at 80 um from the centre is 
0.045 i) 
Field curvature 2 pm at the rim 
Asphericity (S:) +11 um (with respect to the best fit sphere) 
Diameter 4.1mm 
Thickness 3.2 mm 
Optical throw 26.0 mm 
Weight <0.25 g (mount included) 


Corning Molded Glass Aspheric Lens 


—— 5.4 mm Diameter + .015—— 


: fess than 6 mm 
Operating wavelength: 780 nm 
Effective focal length: © 4.47 mm 
Entrance pupil diameter: 4.2 mm 


Magnification: infinite conjugate 

Working distance: 1.85 mm 

Field of view: 0.180 mm 

Mass: 0.248 

RMS wavefront aberration: jess than es dat full field 
Q= 48.62° 


NA = Sin «= 0:75 
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LASERS: Future 


Extending and improving current technology should 
allow the reliable, CW powers to be increased 
£rom 20-30 mw. to 50 mw. 

= If non-absorbing mirrors prove practical, it may 
be possible to make single-element laser 
Giodes that operate reliably at 100 mw. 

ps Phase-locked arrays have potential to exceed 100 mw - 


but nave Formidable Probiems. 


-- arrays not yet demonstrated - 
- Diffraction - Limited, 


- Stable, Singile-Beam in 
C.wW. operation 


-- long-term reliability is unknown 


Cdissgipation of heat more difficult than in 
single lasers) 


-- failure of one laser could alter beam's 
pattern. 


* New Ideas include:- Single laser diode with etched 
convex mirrors acting 4s 
wunstable resonator. 


- Coupling of individual lasers 
via external cavity so thact 
they operate as a single source- 
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Compact Disk CD , '20nm. } 
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Corning Molded Glass Aspheric Lens 


—— 5.4 mm Diameter + .015 —— 


Diameter: less than 6 mm 
Operating wavelength: 780 nm 
Effective focal length: 4.47 mm 
Entrance pupil diameter: 4.2mm 
Magnification: infinite conjugate 
Working distance: 1.85 mm 
Field of view: | 0.150 mm 
Mass: 0.24 g | 
RMS wavefront aberration: less than 0.05 J at full field 
O 
Q= 48-62 


NA = sinQ= 0-15 


cylinder lenses 


~want Cicwlor aot Ali Chicoh beams 
Les lenses are citolar 


use ANAMORPHIC SYSTEMS. 


— light incidenk on prism at Brendes Angle 
Ok light is Tans miked 

- — fesult is an ec pansion aT tke eam in 

° the ditedOn SNown , a CG factor 0} N 

| — uSing Two prisms , the eon emeoes 

a ction paralhe) to ts ori ainah Path, 


rocus Sensing MetHons | y 
Aa fput-22PWS  Deteks D, od T, ate eqpalley illuminated 
. 3 onley uhen fhe Knife - edog is at the 
Sock foink J the mirtor, ie. Az=O 


Mirror 


Ailherence Siqneh %S a non-linear Lundion of At 
pormennes 


~ 
“= 
“_ 


\ 


eee 
eeae*?* 
ae 
ove 
eee 
ee 
eee 
ese 
wnat? 
ee 
ae 
soe 
eee 
oon 
oe 
coe 
ove 
oe 
2 
voor? 
oes 
eetee 
eee ee 


Net Signal = A-B | 


/ 
“ 
v 
Pa 
. 
o 
7 
’ 
A 
f 
, 
7 
r) 
%e 
ern ener ete tetany renee. wee 
en BOCEOENOES 6000 SOSTEE SES TOECEHS ESS DOSS OOSO HOES UO VERDES LOR A . 
q 


SSM 
> 


esee tt? 
etee 
noeecone ooo" 
eae 
ltrs 
econ eee 


BrA 


eoee* 
eee 
ones POOe 


ft 
Waa 
oO 


e 
° 
og vo censvcceee peseevnevencoons 


a4 
rs 


Bi pr SN a 


_ 


TR vPRisy VARIATION 


Net Signal = A-B-C+D > 


F.tap- Outy Grrea, Dex 


RN » ) Pre-Embossed 
Aan, ¢ Data Pits — 


w 
= > BAS 
Uy) \ iP ) 


TY rehetetinad POEM BOG DPLP PEMB Gig . 
. (sti Hy Wi Z Sos, 


500A Aluminum Reflective Layers’ | i 


eH ap 
. ay Disk Substrate pe : 
| (PMMA) 


FY If 


Laser Read” “ 
Beam 


6 Spot size loregr thon pik width 
. pik edoes are bekecked 


625" 0- ROM 


(98 \ L Aeerevenen Ly 
~ ZOcm. diameter PMMA disk 


4 hr. video 


trex Compacr flupo Disk | cD 


—. \QOmm. diam, YorycARBovaTe 


om me “TS ening, ; O00 


CD-ROM 


= panenen diam, PC 


— SSO MB «a af data 


Track Density - Optical Storage 


P e To Limit Crosstalk 
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Typical Optical Recording Medi- and Their Writing Mechanism. 


1. Ablation * 
WH. CW. Te-alloys, organic dyes, 
Ag in gelatin 


2. Topography Change 
WAAL WWW . 
, Textured Ge, Si 


LL ily, Ton Au (or Pt) alloys/organic polymer 


Gneri 


3. Chemical Reaction 


Rh/Si, Se/Bi 


4. Particle Coalescence 
—) Discontinuous Au 


5. Phase Change 


TMA. “Cnalcogeride, Te0, 


| 6. Magneto-Optic 
e44as hangs Mn-Bi, Tb-Fe, Ga-Tb-Fe 
7. Seareqation 


TeQ, , Au/Ted, 


3. 


4. 


A 


. -To permit unlimited(> 


~ Basic Read/Write Medium Requirements 
: | 60 


. High laser writing sensitivity 


-Compatible with diode laser recording at 1-3MB/s 20) 


. Well defined threshold for recording 


Ins | 
100M) non-destructive reading 


High resolution and regular mark geometry 
~- To achieve adequate raw BER at high density 


*: 


High contrast between recorded mark and unrecorded track 


-To maximize read signal amplitude(SNR) 


. Long lifetime 


-To maintain low raw BER(<1/10K) for at least 10 years 


. No intermediate processing , postable (dakn can be added later ’ 


-To permit data verification during recording 


- Comparable with formattable substrate 


-For track servo, sector ID, etc. 
—To cerify uauritten mediq 


- Compatible with low cost, high volume manufacturing process 


\ 
CNR 


; > 
1 2,3 4 5,6 7 mW, 


¢oad write. 


Laser Power (mW) 


— Speci fea by laser manul- 


Weiti Hat Tower Density [em Wy 


— depends on laser powes 
and ON ophcs | 


Writing Enequ 5 
_ depends on laser power 
_. OOFics 


write- Once MEDIA. ——— 
ee mee —t, SSS q37°_°w~wqxi ISIS] Strongly absorptive 
“ YF recording medium 
50 
Zz 40 | | a | (a)| 
. s—~s_. Bs. > ae | . 
a 30 Approx. An 
t | Moderately absorptive 
20 . A 9 recording medium (n) 
| | h Poe ee me ALUMINUM ref| ector 
10 | 
ene 
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0 100 200 300 400 ~~ 500 | (b) | 
thickness (A) | | \* 
Notmal pe ace (nearly) ' less- dense ka mre. dence Approx. 4n dectrod ae S:m('c (OV d: tux Keer, A\ 


Strongly absorptive 
+-ASSSSSSI—<—recording medium 


—~— Transparent spacer 
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OrricAL Disk Sreuctures 


Air Gap 


Active Layers 


CAA CMB BRARAARAMRAABERMP ERB AAP DAMP PBA MAMMA ODRBMBMAD OOH OM 
CS  e 
Substrate 
(Glass or Plastic) 


(a) Air Sandwich 


Active Layers 


Overcoat 
LLL LLL LLL LLL LLL LLL LLL 
Substrate . 

(Glass, Plastic or Aluminum) 


(b) Encapsulated 


| Substrate. | 
3 (Glass or Plastic) | 


| Air Sand wich 


Used with ¢ 
= ablation 


— Bump Sorming 
= Topoqphy Change 


E ncaps ulated 

used with ; 

— (200- only 

— chemicoh Ceackon 
— AiScontauovs Golo 


— +hermo- MANNEKD ~ootic 


REVERSIBLE OPTICAT, DATA STORAGE 
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Advantages 


Lower cost of media 


Compatible with magnetic recording 


Multi-Function drives 


Cread-onily, 
reversible ) 


write-omce, 


Sd 


Disadvantage 


Erasure is not easily detectable 
tssues 


| 3 
Limited erasability CAy 10 cycles) may be 
acceptable for storage of images, audio, 
but not for data 


Erasure should be fast so that writing data 
rates match reading data rates 


G. Bate 
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Amorphous to Crystalline Phase Change in Te-Based Materials 
eg TeAsSe, TeGe, Te 


‘Crystalline 


Ann ——» Reverse-Mode Writing: 


Amorphised marks along 


a polycrystalline track 


. BELL 


Phase-Change Recording 
[TeO, 4] Ges Sng 


Twi-—---f-—-—-—-—------------------- 
Rapid Cooli 
Temperature — 
T@+---t+-----—---------- ———}—-_—_ 
Write. Erase 
Read Slow cooling 
C AA A A — C 


| O: (0: 0:0 
C=crystalline phase © © © high reflectivity, low transmissivity 
OeeECE® 


A=amorphous phase en : - low reflectivity, high transmissivity 
O 
=“ = G. BATE 


Ideal Optical Storage Material 


= Need a material: —writable with low-power laser 
—high reading SENSITIVITY 


—but STABLE with respect to TIME 
TEMPERATURE 
HUMIDITY 


STRAY FIELDS 


=» Medium must also be STABLE and REVERSIBLE 
after millions of write-read-write cycles 


Answer: Thermo-Magneto-Optics, TMO 


Writing and Reading Energies 


Write-Once Media 7 Reversible Media 
Nhp 
Energy to write 
Nh» (writing) nhi 
oe nhy (reading) Bias -/ 
7 | | _/ 
7 } 
ya 
‘‘Q”’ state ‘“1’’ state ‘<Q”’ state ‘4’ state 


e.g. Hole-burning e.g. Thermo-Magneto-Opti 


CyrRié Pwr 


THermo — MAsNETIC RECORDING 


@ 
Te Tq Te SIMA SeREALALIE 


Magnet 4 


Bias 7 H. eg eee oe 


i 
Ta Te Temperakure 
= 2S C [v 200°C | | 


HATER TYE 


@) 


rye dette hay ytt 


me 
=> 


TA 


(Og 


TMO RECORDING | 
(COMPENSATION our ) 


T=MxH 
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Properties of Rare-Earth, 
Transition-Metal Alloys - —=— 


Tb, Fey. 
x Tx 0.3 
e Kerr Rotation Determined by TM Magnetization 500 L “ i, —> / 
e 6 Ranges 0.25°-0.5° Depending on Composition " 
400 |-‘e | “N02 


Mt a ; 6 
emu/cc 300 a / degrees 
$s ze | 


@ 
Curie Compensation pedis Anisotropy Kerr 200 \ 0.1 
Temp - Temp : mae Constant | Rotation e 
| (K). {KY | gga itensy | (104i/m3) | (deg) | - 
7 100 


Te Teoma 


Te 

eee Tb23Fe77 400 300 7 
[Teaco | -- | 300 | * | | 033 ons 20 | 25 30 35 
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eaters |= [300] a ae Ss 
PTenaFesscorn | ~500| <0 8 | 10 | 038 
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T. M.O. Materials : Kerr Rotation, ©, 


ear Material 8, Wavelength Temperature 
1958 MnBi 0.7° visible 25°C 
1985  GdTeFeCo 0.45° 780 nm 25 
1987 ~ PtMnSb 1.9° 780 25 
1986 CeSb, _.Te, 1. 3° 1.2 ym 2K 


1986 CeSb 14° 2.4 2 


G.Bate 6/26/89 


Enhancemer. >f Kerr Rotation 


| nD | ry) P ) 
iG reduced ss enhanced 
Incident: on 
Plane polarized 
M-O 
Layer 


Single layer ; 
: GdTbFe, A = 800 nm 
Reflected: 7 Bilayer Ohta et al, Sharp (1983) 


Elliptically polarized Trilayer 


Quadrilayer 


Magneto-optic Ant-reflection Trilayer 
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ae 


«—Absorber—> 


«—Dielectric—>» 


«—Reflector— 


Destructive Interference Constructive Interference 
ry - 0 r. Enhanced 
High Absorption of Increased Playback SNR 


Record-Beam Energy 


ee Aree ne cman Plasma. 


| @) P,, * Argon Pressure 


@® V; # Target Voltage 
&) Va = Bias Voltage 
Target Composition 


Maanetton Soukkerin 


used to increase deposition 
(okes 
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Planar Maa nek TOA 
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| Alloy targets Fe Tb 
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Properties of Interest in T.M.O. Films 


Writing 
Te; Tcomp ,M Ss 
a (A), R(A) 
ae ae) 
Tcomp 
thickness 
specific heat 


therm. conductivity 
(film and substrate) 


spot diameter 


Storage 


H,, M, (20°C) 
AH¢, AM, (x, T) 
Ku > 2TT M,2 


magnetostatic vs. 


wall energy | 
corrosion 


Reading 


u* = n — ik (Temp.) 
[a(A), R(A)] 


thickness 
Of, Ox : Or, Ox R 
: OL 


T max 


grain structure 


interference layers 
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Domain Stability 


Demagnetizing Field (H,) 


— =xpancs Bits wells He : i 
mAVAVA TTT TT APT TD NN dy ttciemer a 


Wall Field (Hy) 
Contracts Bits 


@ Stability Criterion: Heee =IHp - Hywi<H, - 
a 


© Minimum Stable Domain: d_,, = wT 
Cc 


DOMAIN STABILITY 


(Nommatizes “Tree ("Seaton ) (rata) (Ft) 
AE (4) 2-2 + amm?. De (2)-omge 
Ar 2mrh r r h . 


FOR COERCIVE STABILITY: 


AE | : : 
Semana gy ea] >= MH 
Ar Fy, C =. 
OR: 
r ed (Small M, thin films) 


min “ 2M(He-H) 


TMO Drives 


Company Disk Capacity Access RPM Data 
Diameter (2 - sides) Time Rate 
MAXTOR 5.25" 1GB- 35 ms 2,200 10 Mb/sec. 
650 MB 
SONY | 5.25" 650 90 2,400 7.4 
CANNON 5.25" 512 80 3,000 9.1 
OLYMPUS/RICOH 5.25" 650 44 1,800 5.6 
SHARP 5.25" 650 67 2,400 7.4 
H-P 5.25" 650 95 2,400 5.1 
5.6 
VERBATIM 3.25 " 64 MB 30 2,400 2.1 
(1 - side) 


G. Knight 


OPTICAL DRIVE TECHNOLOGY TRENDS: 


e G.E. has demonstrated an 8-track laser device 


@ frequency - variable lasers are being developed for use with 
multiple-layer films | 


e seek-time now < 20 ms. 


@ rotational speed —— 3,600 rpm 
@ higher-level of electronic integration 


1995 | data rate > 3MB/sec. 


access time < 15 msec on3.5" drives 
20 msec on 5.25" drives 


capacity > 400 MB 0n3.5" disk 
> 2GBo0n 5.25" disk — 


cost <. $1/MB 


G.Bate 
Dec.1990 
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Data Storage in 2000—Trends in Data Storage 
Technologies 


MARK H. KRYDER 


MAGNETIC Disk DRIVE IN 2000 


Capacity: 10 Gbytes 


HE ipht 3.5-in disks 
— 125 000 bits /in 
7000 tracks /in 
0.7-Gbyte /side 
Zone Bit Recording 


o. aye * « ; 4 e 
og ? 1 ad Mt wn i 
* = ae as ea e.g e . 


Non- cemoyoble disk 


Data Rate: 100 Mbytes /s 


6 Mbytes /s—channel 
at 5000 r/min 
e | head /surface ¢ 


MAGNETOoopTIC Disk DRIVE IN 2000 


Capacity: 10 Gbytes 


7s One 3.5-in disk - 


~ 115 500 bits /in 


~ 


—/ _0.33-nim spot 


1.5 bits / transition 
63 500 tracks /in 

0.33-ym track 

0.07-~m guard band 
5-Gbyte /side 

Zone Bit Recording 


| 4 ‘i 
Pr ae Dhyy 


femovaule disk 


Data Rate: 100 Mbytes /s 


6 Mbytes /s—channel 
at 5000 r/min 
XK 8 lasers /head 4 


| head / surface 


LCSSUES >: 


7 Reversible versus write-omce or read-only 
~ Data rate, overwrite by magnetic Field 
modulation —- need fast erasibility 


= Bit density limited by MN and N.A. 
Magnetization determined by R.E. and T.M. 
Magneto-optic effects determined by T.-M. 


aa Optical Nnheads - integprated?, short / lasers, 
holographic lenses 


as CAV vs. CLV vs. banded CAV 


= Removability, vs. multiple disks/drive 

= Access time Chead mass ) 

a Continuous VS. sample servo 

= Multi-functional drives 

ae Hierarchy of storage must be transparent to users 


as Media cost 
bad Media life expectancy, reversibility 


= Media backward compatibility 


TMO versus Conventional Magnetic Recording 


TMO: MAGNETIC RECORDING: 

e much greater head—disk | ® mature, established tech. 
separation @ problems are known 

® no risk of head crashes e standards established 

® removable disks e small, light heads 

e tpi x 10 higher (—> faster access times) 

e lens can handle multiple e multiple disks /spindle 
signals e® high data rates 

e bits are passive until * deck cveraite 
interrogated — 


e stable media 


G. Bate 
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FLEXIBLE OPTICAL DISK SYSTEMS 


~ Bernoulli Plate 


Motor Hub Microhead 


« 


Static Disk 


Bernoulli Plate | Motor Hub Microhead — 


Spinning Disk 


[30 NM aAisk 
pia disks /cackei d0¢ LGB 
AO ms. AVETAGe accesS inne 


1:5 Bute [sec daka roke 


_ ot heioht 


© Bernowlli principle stabilizes 
disk — head separ On. 
~ No Junamic Joss ng regpired 
— Small , liar Weight head, high NA, 
e Iwo dicks Can be SPar Nn ONe 


Cartridoe 


IC| DIGITAL PAPER 


e FET suloetrate. 2S —7S pm Mick 
700 

e LR. alosorb' NQ » ove oclyner{agonm 

spot size e pre Format: a ol Aaka & traces 

Da em Dossing Ove - polymer layer 


6 WCU 19g DY pit— forming ro 
Melinex”’ Oue = polymer Laver, Am puts 


polyester 
substrate 


; lite expectancy . a years 


Metallic reflector 


Protective overcoat Dye/Polymer layer 


° ce cording crore of 10 Miz VSINA a LOmW laser 
. los. reel ol alt ind Ta Pe Wolds G00 GS 
(cf. magnekic tape 180 MB) 


IC! 1012 OPTICAL TAPE 


Laser Beam 


| Meiallic reflector 


Protective overcoat 


Dye/Polymer layer 


‘Melinex®' polyester 
film substrate 


PRODUCT DESCRIPTION 


TAPE PARAMETERS REEL PARAMETERS 
Length 880m Reel [.D. 

Width 35mm Reel O.D. 

Thickness 0.075mm Flange Spacing 

Elastic Modulus 4500N/mm Flange Thickness 
Substrate ‘Melinex®’ polyester film Reel Composition 
STORAGE AND OPERATING REQUIREMENTS 

Long Term Storage Operating 

Optimum Temperature 18 deg C Optimum Temperature 
Maximum Temperature 65 deg C Maximum Temperature 
Minimum Temperature 5 deg C Minimum Temperature 
Maximum RH 93% Maximum RH 


READ/WRITE CHARACTERISTICS 


Lasers Power: 10 milltwatt Wavelength: 830nm 
CNR 55daB using 830nm laser at 7 Meters per second 
‘Raw BER 1x10 -5 to 1x 10-6 


COMPATABILITY 
Creo 1003 Optical Tape Drive and equivalents 


76mm 
‘317mm 
37mm 
2mm 
Chemcor 


18 deg C 
32 deg C 
16 deg C 
80% 


addition of 1-8 at. 7 of [Pe], Ce. AL 


6 undec coating and Overcoating eq. Si 0, [ALN] 
| peels 


e addition of B00. % Th to Sid, 
ae. a or 
1) =41-9) rn <Q | 


| . oO S 
8. O04 —» 06 
Chromium (absoc ber) 


o ‘el > Goraek : high - , low H. i Ky 
Ss seneieetr eit meee, «3 low O, high HS // Ky 
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T>Teoomp 


Demag Fields 


A Domain Wallis Created 
(T > Tcomp ) 


Te] 


in Walls 


Second Pulse on the Written Domain» 
( T > Teomp) : 


A Domain Wall is Created (T>Teomp) 


The Two Domain Walls Annihilate 
by Wall Motion (T> Teomp) 


The Domain has been Erased ( T= Ta) 


HPS 


WORLDWIDE REVENUE WORLDWIDE SHIPMENTS 


om ALL OPTICAL DRIVES sy ALL OPTICAL DRIVES 
2.0 | 
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Storage Interfaces 
and Architectures 


IST. 


SCSI-2 
@ Overview 
@ New Features 
Storage Architecture 
@ Evolution 
@ Optimizing Data Transfers 


@ Disk Arrays 


e Disk Buffers and Caches 
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O Overview 

© New Features 
Storage Architecture 

@ Evolution 


© Optimizing Data Transfers 


© Disk Arrays 


© Disk Buffers and Caches 


TIST. 


Storage Architecture (1975) 


Host 


HOSE Hee eae Bridge 
omputer Adapter | 4 Controller 
Controller and 
Channel Adapter Control Unit Device 
Host Bus 1/O Channel Device Interface 
PDP-11 Q-BUS MASSBUS SMD 


iSBC Multibus SAS|I ST-506 


HIST 


raf Storage Architecture (1980) 


Single Board Controller 


; Host : 
H OS t "ee , . . : . oS 
es Bus o 
Ss Controller - 


Adapter 


Host Bus Device Interface 
LSI-11 Q-Bus SMD — 
iSBC Multibus | ST-506 
IBM PC PC BUS ST-506 
IBM PS/2 Micro Channel ESDI 


lIST 


Storage Architecture (1985) 


See Re te 8 "* 

see SEN 
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| Embedded 
2) Controller 


Host 
Computer 


Host Bus Me} Channel 
FC PC-AT Bus SCSI 
PS/2 Micro Channel SCSI 
sun = +~VME bus SCSI 
Vax BI bus SCSI 
IBM 4300 370 oemi 
| IP 1-3 
HPPI 


HIST 


Storage Architecture (1990) 


Embedded 
— f= | HostBus 
t Adapter 


Embedded 
Controller 


C O mM p U te r ESORR 


Host Bus 


@ Put the HBA and controller on the Disk Drive 


@ Lower cost, higher reliability 


@ Most popular System Bus is IBM PC-AT (ISA) bus 


® Motherboard manufacturers now including IDE connector 
with signals and power for attachment 


lIiIsT 


2 af mm Fe 


ST-506__ 


ESDI ———_ Are Device Level Interfaces Dead ? 


_ 


SMD-E 


LIST 


Existing storage architectures require DLI drives for expansion 


Most vendors continue to sell DLI drives 


¥ athena ttf MR 


Customers not designing DLI drives into new systems 


Primary emphasis away from DLI drives 
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SCSI-2 
° Overview 
© New Features 


Storage Architecture 


© Evolution 


© Optimizing Data Transfers 


© Disk Arrays 


© Disk Buffers and Caches 
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SCSI Vocabulary 


Initiator - Host, Host DMA channel, HBA 


Target - Controller/Control Unit and Device 
CDB - Command aiid shod Block 

LUN - Logical Unit Number 

LBA Logical Block Address 


HBA Host Bus Adapter 


Ist 


SCSI People 


Andy Hospodor 


HST Santa Clara Univ (408) 554-6853 


John Lohmeyer 


—— (316) 636-8703 
SCSI-2 Chairman 


SCSI Bulletin Board 


: (316) 636-8700 
Download the latest spec 


LIST 


SCSI Phases 


SCSI Control Lines 
SEL BSY C/D 1/0 MSG 


Bus Free oO oO CO oO © 
Select @®@e@ooo 
Command 088 @oo 


| Data In (Read from Disk) O @ oO @ OO 
Out (Write to Disk) O @ O dO O 


Status Oee@ee@eoo 


Message O08 @ 6 @ 


SCSI Phases 


8 _ Mf Petrr 


a 

m 

Oo 
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SCSI Req-Ack Handshake 


Target Requests a Transfer 


Initiator Acknowledges Transfer 


Target Removes /REQ 


Initiator Removes /ACK 


ari am 


LIST 


What Made SCSI 
take off 


Industry moved away from ST-506. Customers 
wanted data separators on the disk drives, leading 
to integration of the Read/Write Channel. 


Microprocessor technology allowed a logical 
Interface, error correction, and defect mgt. 


Custom VLSI technology allowed electronics to 
fit on the disk drive. Controller, Servo, and 
Motor speed electronics integrated in 3-5 ICs. 


Embedded controllers allow a generic interface 
to the host, while provided an optimal disk 
interface. 


Hey, What About 
IPI ? 


rIST. 


A tLlaenadcar 


IPI-2 is a device level interface 


IPI-3 is an I/O channel or system level interface 


Currently two domestic manufacturers shipping 
IPI-2 products. Only one IPI-3 product shipping 


Many manufacturers announcing IPI versions of 
new products 


Chip sets are sole sourced and priced prohibitively 


"If tt doesn't look like a washing machine... 


and doesn't run on 220 volts... 


then it isn't a high performance solution" 


SGS|-2 - 


Standard Expected:1990 


New Mandatory Commands 


Separate command sets for: Disk, Tape, WORM, CD-ROM 


New Messages 


. Command Queuing 


Faster Data Transfer Rates 


SCSI Command Descriptor Block 


6 Byte Command 


~ | OPCODE 


LUN | LBA 


XFER len | 


10 Byte Command 


resrv 
L 
L 
L 
L 


BA 
BA 
BA 
BA 
reserved 
XFER len 
XFER len 


Control 


OPCODE - operation to be performed 


LUN - Logical Unit Number 


-0 for embedded SCSI drives 


LBA - Logical Block Address 


© 


oO aA nN OO CO fF WD ND — 


XFER len - Number of blocks to be transferred 


IIST 


A LUlaenadar 49/0 


SCSI Command Descriptor Block 


Examples 
6 Byte Command 10 Byte Command i 
| 2A jo 8 
2 
aE 
48 fa 
: 
Read one block 800 vA 
Starting at block 14C363 
oo fs 
: 


Write 8000h blocks 
starting at block 0123456F 


Ten Byte commands address larger devices and offer 
longer transfer counts, making them more suitable for 
backup operations. 


ae eae eee eee en ea RE ae nae ocr es Sn 


Mandatory Messages 


SCSI (1986) 


Command Complete 


SCSI-2 (1990) 


Abort 

Bus Device Reset 
Command Complete 
Identify 


Initiator Detected Error 


Message Parity Error 
Message Reject 


No Operation 


Mandatory Commands 


St ee) ~ SCSI-2 (1990) 
Test Unit Ready Inquiry 


Test Unit Ready 


Request Sense | 
send Diagnostic 


Format _ Request Sense 
Format 
Read > , 
Read 
Write 
Write 
Reserve 


Release 


1 
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Command Queuing 


Untagged queuing allows single commands from 
multiple initiators to be queued 


Tagged queuing allows multiple commands from 
one or more initiators to be queued. Each command has 
an associated Queue Tag 


a. Commands with Ordered Queue Tags must be 
executed in the order received. (FIFO) 


b. Commands with Head of Queue Tags are placed 
_ first in queue, to be executed next. (LIFO) 


c. Commands with Simple Tags are executed in 


an order determined by the target. 


Advantage: Seek Optimization 


Disadvantage: Indefinite Postponement 


Data Transfer Rates 


1. Asynchronous SCSI 2.0 MBytes/s 


2. Synchronous’ SCSI 5.0 MBytes/s 


3. Fast SCSI-2 * 10.0 MBytes’: 
4. Wide SCSI-2 * 40.0 MBytes/: | 


“ - Requires Differential Drivers on SCSI cable 
Maximum data transfer rate also limited by: 
Noise 
Cable Impendance 


Termination 


LIST 


New SCSI-2 
Mode Select Pages 


Caching Page 


Notch and Partition Page 
Read/Write Error Recovery Page 


Verify Error Recovery Page 


Pils f 2, 


SCSI-2 | 
Reconfigurability 


LIST 


@ Medical Imaging System 


@ Database system search/sort/query 
© Workstation used for I.C. layout 


© Verify Error Recovery Page 


Application software can use Mode Select to 


Dptimize caching strategy for specific job. 
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OQ Overview 
© New Features 
storage Architecture 
© Evolution 
@ Optimizing Data Transfers 


© Disk Arrays 


© Disk Buffers and Caches 
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Data Transfer Rates 
Comparison of SCSI and SCSI-2 


Asynchronous SCSI 


LIST 
A WHnenndor 12/9 


/ACK 
—> le— 500 ns — 
Synchronous SCSI 
/REQ 
/AQK 
+ oe —> — 200 ns 


Fast Synchronous SCSI-2 


= yh ta — l¢— 100 ns 


Thread Management and 
Data Path Management 


-_ 


210 MBytes/second 


Mee UT A 


Data Transfer Idle 


5/0 MBytes/second / 


/ 
lll ll” ll Ill 


MIs T 


DBO-7 


Shorter periods of high activity with longer idle periods 


Why not relinquish the bus during these idle periods ? 


1. Target must send disconnect message to initiator 200 


2. Iniatiator must send ID message to Target 20C 
3. Target must reconnect — 101 
4. Initiator must restore pointers & resume 30C 


eed 800 microseconds between blocks to use disconnect/reconnect 


| Efficient 


Thread Management and L 
Data Path Management , 


& 
£ 
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Data Transfer Idle 
5.¢ MBytes/second / | 


Disconnect Reconnect 


Large data transfers may be packed together, allowing efficient disconnect 
reconnect. In this example 4 blocks are transferred during each burst, then 
the SCSI bus is relinquished. 


A winchester disk with 1 block = 512 bytes, and 32 blocks/track, 3600 rpm 
Each block would take 16.67/32 = 520 microseconds to read from disk 


A group of 4 blocks would take 2080 microseconds to read from disk 
and would take 409 microseconds to burst over SCSI 


Leaving about 1500 microseconds to disconnect/reconnect 


Notched Recording 


Inner Track 


28 Sectors per Track 
3600 RPM, 7.5 Mb/s 


Outer Track 


36 Sectors per Track 
a 3600 RPM, 9.6 Mb/s 
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Notched Recording 


Data Rates 
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Zone 0 
Zone 1 


Zone 16 


Outer Track Inner Track 
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striped Disks 


Embedded 
Controller 


Host 


Embedded 
Controller 


Embedded 
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Embedded 
Controller 


Increases Throughput performance 


Files are broken up and spread out 


Least reliable way to store data 
may be implemented with mirrorin 
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@ Disk Buffers and Caches 


Where Should the Cache go ? 


On the Disk Read/Write channel 
On the Disk Controller 


On the Host Bus Adapter 


In Main Memory 


Semiconductor Memory 
On Disk Drives 


Minimum RAM data buffer - commodity drives 
Caching RAM data buffer - high end drives 


All RAM data buffer. no disk - Solid State Disk 


How much buffer does your system require ? 


Buffer Size versus Performance 
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Disk Drive Buffer Sizes 
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1981 1983 1985 1987 1989 1991 


Low end: 2k-8K Static Ram Buffers 


High end: 1Mx9 DRAM SIP array caching buffer 
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Cache Performance 


The access time equation: bs : 
VY). 
= 

| — 

access _ overhead + seek vs Detees transfer 
in the ideal case: 
transfer 


access overhead cache_overhead 


Let p(h) = probability of a cache hit 
p(m) = 1 - p(h) 


| t = + + 
access overhead seek rot_latency transfer 


+ p(h) [' ~t « ] 
cache_overhead seek rot_latency 


Design question: 
What is the minimum p(h) required to break even ? 


“| or How much Cache Overhead can be tolerated ? 
: : J 


Tuning a Cache Buffer 


Cache Mode Select Page 


Enabling Fast Write where appropriate 
Minimum prefetch length 


Maximum prefetch length 


Disconnect-Reconnect Page 


Buffer Full/Empty Ratios 


~Disconnect/Connect time limits 


Maximum length of data bursts 


Application software has the ability to control: 


What data goes into the cache 


How and when the data comes out of the cache 


| Predictions 


990 - 2000 


1} 5.25" read/write channels will expand to envelope ~ 
SCSI synchronous data rate (5.0 MByte/sec) 


2| Hardwired decoding of SCSI commands on chip sets 
31 Cache buffers using 1Mbx9 and 4Mbx9 DRAM arrays 
4| Integration of SSD & Mag Disk = Mixed media devices 
ot Multiple head/slider - Paired Read/Write channels 


6| Disk arrays - for those who can afford them 


“| Key performance issue: Rotational Latency 
Higher Rotation Rates (4800, 5400 RPM) 
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TRENDS IN DATA STORAGE 
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- high Kerr-rotation media 

- multi-layer films 
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COMPUTER MEMORY EVOLUTION 
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Offline Storage. 


Constraint: Storage Technology 
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CAPACITY SHIPMENT SUMMARY 
Worldwide Shipments in Terabytes 
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Secondary Storage Market Growth Estimates (in OEM $) 


$M = =$M 

Optical 1988 1993 CAGI 

Autochangers | 
Optical Autochangers 43 B51 82% 
Rewritable 9 1200 166% 
WORM 145 437 25% 
CD-ROM 42 404 57% 
Minl-cartridge 153 383 8% 
1/4° Cartridge 333 348 1% 
DAT 2 351 181% 
8 MM 18 120 46% 
1/2" CART* >1.5 MBS 407 786 14% 
1/2". CART* <1.5 MBS 156 — §45 28% 
7 1600 Tape 219 37 -30% 
Trad. 6250 Tape 358 31 -39% 


Low-Cost 6250 Tape 222 143 -~8% 


. 
| ow-Coat 6250 Low-Ooet 6250 TOTAL A171 5942 22% 
1988 1993 SOUrces snowing — Reports 


1987 BIS Macintosh 


et 
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Growth % Dominated by 3480 cartridges 


Industry Structure: 


LIMITED VERTICAL INTEGRATION... Still 
COMPONENT VENDORS OWN CRITICAL TECHNOLOGIES 


DRIVE MANUFACTURERS 


Today: Consolidation / Partnerships 
Non-Captive -- boxes 


Captive -- storage systems 


TRENDS: R&D ... USA. 


MANUFACTURING ... OFF-SHORE 
HEADS & MEDIA ... JAPAN 
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OEM Price per Megabyte 


PRICE ($/MB) 


1994 


A @ fin nti Mam 


State of the Art: 


Disk OXIDE to ‘THIN FILM (startups) 
5.25" to. 3.5" (Apple) 
OPENLOOP _ to SERVO 
DEDICATED to | IMBEDDED 
ST506 to SCSI 
MTBF 5,000-10,000hrs. to 50,000 hrs. 

Tape LONGITUDINAL to ROTARY 
Optical TALK to DELIVERIES 
Example of change: 5.25" disk 5 MB to 1.2GB 
| 85 ms to 11msec. 


5 MBits/sec. to 2+ MBytes/sec. 
$20-30/MB to $2-3/MB 


FIST 
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Technology Dynamics: 


MAINSTREAM 
Thin Film Disk (startups) 
MIG heads (Japan) 


UNDERWAY 
Glass Disk 
ZBR 
Rotary Drive 
2.5" Hard Disk 
2" Drive 
Optical Servoing 
MR Head 
Horizontal Head 
MET 


Disk Arrays 


. Areal (NSG) 

... Imprimis 

. Exabyte 
... Prairetek 

. Sony 
... Insite 
... PCI 

. Let 

. Japan, Inc 


; Everyone | 
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AREA DENSITY 
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ony Corporate Research Labs 
Institute for Information Storage Technology 
Workshop V October 11-14, 1989 
Presented by Senri Miyaoka 


YEARLY ADVANCE of RECORDING DENSITY 
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RAMAC350 


Drum Diameter 
Wrap Angle 
Track Width 
Gap Depth 


Head Tape Speed 
(ips) 


FC] 
Azimuth Angle 


Tape Speed 
(ips) 


*Hi8 
*Exabyte 


HELICAL SCAN RECORDING 


8mm 
40mm 


221° 


25 micron 


-15 micron 


148 
77,000* 


+10° 


96 


103,800 
54,000 


20 micron 


15 micron 


124 
76,200 


+20° 


2 


VHS (1/2") 


62mm 
180° 
32 micron 


30 micron 


230 


41,700 


£O 


1.31 
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CAPACITY INCREASE 
(DC-600 CARTRIDGE) 


CAPACITY (MB) 


400 
QIc-320 
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REMOVABLE HARD DISK DRIVES 


Evolution of Technology and Capacity Examples 


Drive Capacity/Pack 


IBM 1131 1.024 MB 
Burroughs 9480 4.68 MB 
Ampex DM-442 3.125 MB 


IBM 3340 34.9 MB 
CI] D120 10 MB 
CDC Lark 8.35 MB 
DMA 6.75 MB 
SQ306 6.38 MB 
507 12.76 MB 


9Q555 54.78 MB 


# of disks 


| re] a — — 


ee , ! rs oe es er eo a ae | 


TPI BPI 
100 1100 
100 4400 
200 2200 
300 5636 
500 4750 
237 10161 
454 8617 
435 12000 
741 12608 
1086 23642 
P3920 
iis T—— 
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Heads & Media -- Long Term 
As Track Density Limits Get More Intensive Focus 


a Will require "pre-formatted" substrates (a la optical disk) 
m Disks in turn will depend upon multi-element transducers 
To servo during read and write - 
To separately optimize read and write elements 


To permit erase option to improve overwrite 


a Will accelerate move to thin film head devices 
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Technology for the Mid and Late 1990's 


Media: Sputtered Thin Co-Alloy Films 
Glass Substrates 
Possibly With Discrete Tracks 


Heads: Integrated Magnetoresistive 
(MR) Read, Inductive Write Heads 


Source: Applied Magnetics Corp. 
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e The User 
: — Lowest cost/bit 


— Read/Write and nonvolatile 


_— Competitive marketplace-large choice 
of product offerings 


je The Manufacturer 


— Multi-billion dollar industry 
— Established production processes 


— Demand for capacity increases faster 
than storage density from technological 
progress—growing market 


: — Areal density still far from ultimate limit: 


— Fairly well understood phenomenon 
— Reversible process—inherently stable 
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Other Directions 


High Velocity | 


"'Contact’’ Recording 
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Vertical Recording | 
Is History Passing It By? 


Need high readback resolution to leverage very narrow transition 
Thus potential depends on very small separation 


Best opportunity is then flexible disk 
(ultimately could involve pole type write head) 


Major alternative is magneto-optical 
i.e., high density and removability with no contact 


Vertical recording now seen as only evolutionary extension 
-- are still far from limits on conventional magnetic recording 
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Magnetic Disk Storage Density 


Areal Density (Bits/in*) 


(factor of 5 every 5 years) Optical recording ae 7 7 
densities —* 
_} © 
o.60O UP 
Currently greatest ® 
rate of progress is i. J 
with 5 1/4 drives C) 
© 
7 
o © Winchester technology 
(lI 5 1/4 announcements 
a O IBM Gigabit 
demonstration 
© IBM disk 
® RAMAC products (high end) 


1960 1965 1970 1975 1980 1985 1990 1995 
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Gigabit Pemionsiranon 


BH Cell Size 
3380K 

@® Scaled MR Head 

@ Low Noise Film Disk 

@ Laboratory Demonstration 

Error Rate 
Gigabit Deme 
+o 1 2 
| Off-track Posttion (um) 
0.1 


36 Mbi/in2? =: 1180 Mbit/In2 


17.2 158 em 
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T1ST WORKSHOP V PRODUCT PROJECTIONS 
First customer shipments: 4th Quarter 1991 or sooner 


ena etal tinatiaiaien tla tant tai aatieeadieadaeine deat tend SEE Ae SSeS SSS See eee eee ee ee ee cae ate mayne sa ewe sew Sep ere a eb Sane Ae we SEN Sew, ew, ee se (See em, SS: IS Se SE, SIN Oe le a A SSS) SS SNA SP, Se OY SN Gen OD SS SR a ew Sem 


PAUMETERS oT END Wt Ft SINGLE SINGLE OPTICAL OPTICAL HWE CAPAC. =o PERF CAUT TAPE =CART TAPE) «=OCART TAPE OHIELICAL HELICAL SEMI - 


5. eo DISK 3.5" DISK 3.5" DISK StH- 4 5" 5. ea 1. 3.5" EAS: {LOPPY FLOPPY 1/4" LOW 1/4" MID 1/2" SCAN -OMM SCAN-DAT CONDUC TOR 
Areal 2 
Density Mb/in 160 160 | 160 | >=]5 550 580 40-45 | 70-100 6.2 12.0 1.4 15 114 30-800 
wanwncnnennnc enenecences fececcnconee|eanecencee eaeeaecnae [enc necccnes |e neeecnce |onccecenee |eoeencnoes Doeeeterr] Doeeereren eae Eoresecren! Eau seeancees Dees 
DPI x 000 65 | 65 | 65 50 31.8 31.8 N/S 40 52 80 | 38 57 61 5.5 - 28 
Stu tecun lees wessekue pehiceue tee aoe ee tees Nee eee enon eueuaee non [aes te eee Pewee boa eeuete ceo lowes oeewane |----------- -~----------|----------- Canwemedwen Peakumenouue 
PI 2500 2500 2500 >= 1500 18.2 18.2 N/S 1.6-2.4 12 16 036 1640 1870 5.5 - 28 
area eee [anna nnn nnnnfenen nen nc ne fence seen feecce seen cece eee [eee snc nce fence tae eeencencene Penceetncnan | eoncece cane ene ennecen | aceecsenane [peecennnean 
On-t ine 
Capac ity-M):F 2400 4200 160 40: 750 200 40 100-160 1350 2000 480 5000 1300 4-12 
( teeceerecenee [----------- |~---------- [areeteeenes [sees eeee eats [asenaeesees |----------- | ----------- |----------- |----------- |----------- |----------- |----------- |----------- 
Average 
Seck Time-ms 10 12 15-17 20 25 22.5 50 15-20 37500 43000 N/S 40000 20000 -l- .2 
wana nnnnnnnn ede ae Ect eee Eee] oc En cc] Ee oe 
Average : . 
Latency-ms 5.6 5.6 8.3 9 +/- 2 12.5 12.5 30 8.3-10 ee ae a ye ae ae 
wanonsnennons |-n----nnnen | enenesencnn | eeec cece efeceenc cee fee cecceccn ence ste cen [ese cnecnnc cence tccnnee [enecencnnce [ecencencnne Penceeescnee Pecnecencnne eccesccsene fenceectecss 
Transfer | 
Rate MB/s 6.0 3.0 2.5-3.0 2.0 1.28 0.83 0.4 1.25-2.0 | 1.5 3.6 6.0 0.5 0.183 4 - 32 
ens ean |- nanan enn en fone n een ce na fence ccc ne cn feces ec cee ea ceecatcce fecnectencen [ecetennatan Peccnsccenan fescensenean| ecoenetecce |onnecetente | oneeetecnen [eeceenencns 
Capacity/ 3 
Volume-*1/ in 15 30 y 4 52 98 24 40 100 260 46 1400 540 6.7 - 20 
a [on -n nanan fencnnennen fence tcc ncee ese cec ences Penns nce cnccenccce eccencnnc feececeeccen fecencceccne [enecescnace | actecseenan|eoeeeecenen [eneennenene [eeeteeanens 
MIBF - Khrs | 
2 250 100 100 70 20 N/S 60 50 5() N/S 30 10 Very high 
esacesseceaes |sox-22 2225 )eqeesaecass|esesccsesse [uses nesses | esse aessees |-eessereces |eeeccneceas |eersecsemee | eeoeesemes [cement eens | aeeesseeees | eemenesnnen | exeeeensees | eoeeseeenne 
OFM Unit 
Drive Cost-$ 3000 1500 300 200 1600 600 125 300-500 550 650 20000 2300 800 256 - B00** 
wn nanncnen nn oeeeteerees peresbesorre] Eersabeecn! Enussucteen) Eoeaticete Exot ecaneny Eaeesusinec De SUGSnaiece occ saiicce] Heese EuspSSOcec| Ene unSiStcane] Euan saa peSiile 
OM Drive 
Cost /M)--$/0 1.35 1.25 2.00 4-5 2.13 3.00 6.25 3-4 0.41 0.2/. 42 0.46 0.62 64 - 200 
eee |= -a-nn nan enne en tecee fee e cence Penecc neces fen cec cece ecctcnc ence | een ee cence fecctc cece ferececennne ene ter cnn eecene cence |eceeteencne [eenetenncen |oeecenen ens 
Media Unit 
Cost-$ (E.U.) eu | =a | a | a 90 } 25 , 10 | 40 - 60 , 25 34 7.0 8.0 6.0 ee 
Sieuia Pena wee LNbw twee eee eee be ee eae ee Cameo ee |-----+----- | -----------] -----------]----------- ----------- |------.----]----------- TEN METERS REN: emer Tene l ten ee 
Media Cost $ * 0.062 | | | 
per MB (E.U.) es mi , ge | nee , 0.06 | 0.125 | 0.25 , 0.35-0.4 0.018 0.0°8 0.015 0.0016 0.0061 -- 
— Jann snanean fee cecnnenn | enen cnn cneeecececceeee fenee tees en fee nce ce ece fencncenccne feat ccetcne Penne tecnan eoeatencees [enetecennne|aceteee cece fenetennnene |eeneennnnes 
Power (Walts) <1 1.5 <: 1.0 | | 
RRR Re eee iteszesaeeee ss Fee eee ee oe: foe et we ow oe oe ee 8: = | it =: Euadenekd teas 1B es2|eccesecence |commeeneene etemenccens | eocnseseese |qceeceueses | eecemnenas | ceneenerees | encenecnses | eeecenacece 
Cunments FS fone 1.625" form 1.625" form 2.5" single Capacity ts for 1 side. 3.5" 3.5" QIC 1350 3480 DR conf ig- 
factor. factor. factor. disk. Bernoullf, ured as a 
2 disks. disk drive, 
8 disks. 8 disks. * For 2 sided media Non-removble. 
TR media & Wemedia & WF media & TF media, | 
heads. heads. leads. MIG heads. **Sub-ass'y 


Bue pees Ja re ya a ca a ce es eRe PORN RN oR TRC Tan! PaAer rE S Tt ROPE rom ane: pe smaeaseme 
Notes: Areal density In Mbits/square Inch. | | | | | 


On-line capacity: U> unformatted, Fe formatted 
N/S = Not specified, -- = Not appropriate 
Media pricing at end user. levels 


ULTIMATE LIMITS 


a» Magnetic Recording 
-- Minimum Domain Size 
-- Medium Signal to Noise Ratio 


Spacing Dominant Parameter 


a Optics 
_-- Limits due to Wavelength of Light 
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MANUFACTURING PROCESSES KEY 


-m Medium quality 
=» Head medium interface 


PRODUCT SIGNIFICANCE 


s Capacity demand 
40% CGR 
Technological progress 
26% CGR 
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SCENARIOS... 


UISIONS OF THE FUTURE | 


THE DISASTER 


OPTICAL 
MAG TAPE 


THE REALITY 


MAG TAPE | OPTICAL 


TOORY 


ARG DISK \ 


THE HOPE 


SS 
MAG DISK 


THE FERR 


$s 
AG TAPE 


THE INCONCEIURBLE 


THE DRERM 
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Mass _ storage 


Increase in storage density <--- 


----% Decreasing cost/bit 


New applications 


Growth of market 


L--—-- [Economies of scale 
Increase in R&D ------------- 


Alternative Storage Technologies 


GEOFFREY BATE : 


INSTITUTE for INFORMATION STORAGE TECHNOLOGY 
Santa Clara University 


CEI, September 1991 


High-Performance Data Storage Systems 


m 50% of the total system cost is in the storage system 
m Bank of America, San Francisco, (4 large computers, IBM 3090) 


e 1 GB of main storage 

e 1 GB of extended storage 

e 416 MB of solid-state disk 

e 96 MB of disk cache (within controller) 
e 750 GB of magnetic disk 

e 2000 x 200 MB tape cartridges/day 


24 hour operation: 12 hour on-line; 12 hour batch 
e 20 - 30% growth in capacity 
e cost/MB and MB/cu.ft. are key metrics 


; 
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R. Katz 


Problems with Optical Recording 


Present long-term corrosion resistance, archival storage 
overwrite data rate 

access time 

standards 


writing by thermal effects rather than photon effects 


Future m high-power, 50 - 100 mw long-life, cheap laser-diodes 
m laser diodes emitting in blue (only gain factor of 2) 
m = multi-layer films 
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Semiconductor Lasers 


ULTRAVIOLET VISIBLE 


Pbx Sm-x Se 
Pb, Sny,Te 
PDS 1-xSey ; 2 | 
Cd, Hgqi-xTe 
Ga AS, SDq-x 
In As, Pt-x 
(ALyGat-x yIny-yAs , | | | | - 
Aly Gay, ASr— 3 | ee 
Ga AS4_y Py jp 7 300K Fe 
In, Gd4_,As 7 : 
. (Al, Gdt-x) yIni-yPe 
Cd Sx Se {-. =" 
Cd, Zn4.,S —————4 
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LATTICE CONSTANT a (A) 


Energy bandgap and lattice constant for two IlI-V solid solutions. 


TMO Drives : Trends 


Access Time: split hds. with galvo. mirror <25ms 
head / surface 


caching controller (Some applications) < 10 ms 


Data Rate: | > 10 Mb/s 


higherrpm > 3,600 rpm 

multichannel; laser array heads 
one head / disk surface 

pulse width modulation 


Direct Overwrite 

Higher Capacity 

New Form Factors : 1/2 - height, 5.25" 
New Substrate Materials for > 3,600 rpm 
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-MAGNETIC RECORDING MAGNETO - OPTIC READING 


Dielectric 
——_$_ aa 
| —— ——s Chromium Dioxide Particles 


~ Substrate 


MAGNETIC RECORDING RESURGENT > 
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1987, |lwasaki et al, 680,000 bp! on CoCr 
1987, Lambert et al, 0.5 pm wide, discrete tracks _ 
1987, R-DAT; digital recording of audio signals, 58,000 bpi x 2000 ” | 
helically on 8 mm tape (1.3GB) | 
1987, Exabyte; digital recording of data 
2.3 GBytes on standard 8 mm video tape cartridge 
1987, StorageTek; 4400 Automated Storage System 
using IBM 3480, 1/2" data tape cartridges 
19.2 Terabytes on 96,000 cartridges 
1988, Insite Peripherals; flexible disk with optical tracking 
1,250 toi, 24,145 bpi, 21.6 MB on 3.5" diameter disk 
1989, Opticard, 50 MB, arced tracks 
1989, Domain Tech., PrairieTek, 2.5" plated disks 
1989, Superconducting Super Collider, 1 GB/sec. 
1990, IBM laboratory demonstration of 1:2 x 109 bits/in? 
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IBM Gigabit Demonstration 


158,000 bpi / Pe 
sr 4.180 x 10 bits /in 


Flying height 0.04 um (1.6 microinch) 


Head gap 0.20pm (8 microinches) 


Medium -- sputtered thin film 
-- 5 layers 
Head -- magneto-resistive read 
-- "giant M/R" material ? 
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Alternative Storage Technologies 


Semiconductor RAM 


Ferroelectric RAM 


Magnetoresistive RAM 


Vertical Bloch Line Memory 


Holographic Memory 
Spectral Hole-Burning Memory 
Ultra-violet Write -- Infra-red Read. 


Two-photon Recording In 3-dimensions. 
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D-RAM: Bits Per Square Inch versus Time 
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D-RAMS versus Disks 


D-RAM cost reduction factors: 


e higher resolution processes 


lower defect densities 


e larger wafers 


fewer components/bit 


multilayer processing 


Dollars per Megabyte 
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14” Hard Disk (-21%) _ 
Write-once optical — 
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SOLID STATE MEMORY — D-RAM 


Year of 
Introduction 


1976 
1981 
1983 
1985 
1987 
1990 
1993 
1996 


Capacity ~Minimum 


(bits) Line-width 
16K 5-7 pm 
64 K 3.0 


2956 K rf 


1M 1.3 


4AM. 07-08 


16 M 0.4- 0.5 
64 M 02-03 
256 M —o, 
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D-RAMS versus Disks 


D-RAM cost reduction factors: 


higher resolution processes 
lower defect densities 
larger wafers 

fewer components/bit 


multilayer processing 


Dollars per Megabyte 
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Write-once optical 
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D-RAMS versus Disks 


D-RAM cost reduction factors: 


higher resolution processes 
lower defect densities 
larger wafers 

fewer components/bit 


multilayer processing 
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100 
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SOLID STATE MEMORY — D-RAM 


Year of Capacity 
Introduction (bits) 
LoeG: 2 °° 16 K 
1981 64 K 
1963 200K 
1985 - 1M 
1987 4M 
1990 16M 
1993 64 M 
1996 256 M 
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FERROELECTRIC MEMORY TECHNOLOGY 


Ferroelectric Polarization Hysteresis of a Ferroelectric Capacltor 


Polarity 
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Ferroelectric Memory Schematic Cross-Sectlon of Ferroelectric 
Drive Memory Cell 
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U.S. Patent No. 4,759,823 
Patents Pending and other Patents Pending 
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Spectral Hole—Burning 


10° bits (holes) in 
frequency domain 

at each 10u diameter 
spot 


10°® spots In each 
1 cm? chip 


ne ws 
cor 16 to 64 chips in Each arm has one 


— py] each Storage wafer wafer with 2-8 GB 
eee Ih : | | 
e| | 50-1000 wafers per PHB 
memory system - 100-80C0 GB 
total. 8'x8'x4' 
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‘1 hree-Dimensional Optical S.vorage Memory 


Dimitri A. PARTHENOPOULOS AND PETER M. RENTZEPIS 
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addressing. BS, beam splitter. 
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| A New Erasable Optical Memory 


Joseph Lindmayer 
Optex Corp., Rockville, Maryinad 
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HIGH-POWER LASER EXCITATION LASER 
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PROCESS 
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PHOTOREFRACTIVE RECORDING 


PHOTOCARRIERS Standing wave (fringes) 
; created by intersecting 
. coherent beams 


rd aie leas Readout beam generates 


Detector NOEX CHANGE reconstructed image beam 


Array 7 WITH PHASE SHIT 
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Readout Beam 
Reconstructed Image 


MCC HOLOSTORE PROJECT 


MCC PATENTED INNOVATIONS 
ARRAY OF CRYSTALLITES (FIBERS) 


- Scalable storage capacity (larger array) 


- Larger capacity (more pages per stack) 
+ Low crosstalk | 


- Lends itself to low cost production 


-Image Beam 
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Vrs -9 Reterence Beam 


NON-DESTRUCTIVE READOUT TECHNIQUE 
- Allows prolonged readout in photorefractive material 
- Billions of reads without signal / noise degradation 


- May result in archival storage 
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OPTICAL CONFIGURATION 
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_ Magnetoresistive RAM | 


Anisotropic magnetoresistance: 
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e high density 
1Mb@ 1.531mx 5 41m cell demonstrated 


e high speed 
-e radiation hard. 


e IC compatible 
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messenger RNA 
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PROTEIN SYNTHESIS 


bond formation comolete new protein 9 tansler 


new protein chain forming 


OUTER PART OF 
CELL (CYTOPLASM) 


messenger RNA 


TWO PAIRS OF BASES IN DNA 
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A = adenine 

G = guanine 
S = sugar 

T = thymine 

C = cytosine 
P = phosphate 


------ ~- hydrogen bonds 
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STM for Storage? | |= 
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SE 
Starting Assisumptions Ste Orearic Molecuie 
St-3S (e.g. DEP”) 
Ss 
«1 STM can reid on atomic level Fe: Sueiin 
; (6.9. Graphit 
» Potential to write on atomic scale : ss 
— 25A - 2504 bit size O Possibility of 10'? - 10" bits/sq. in. 
4),; c 1 | 
10'7-10'4bits/sq.in. O Demonstrated Read/Write/Erase Feasibility 
- Voltage Pulse to Write and Erase 
- Tunneling Current to Read 
Storage Requirements - Induce topographical features 
| on metallic glass 
STM Current DASD (Rh,.Zr_,) 
e real Density-Mb:sq.in. 10°-10° | 62 
| 350A features Induced (melting?) 
e ata Rate | a~is 
Writing speed 100ns | JOns or i 
Reading speed ? 30ns 
e ,\ccess Rate ? 12ms 
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XENON ATOMS spell out the name of the researchers’ corporate benefactor. 


Disk Roughne:s 2 100A on Ni ak AK Stith) wey ee 


Head Dimensions ? 10% 
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STM for Storage? ad 


Starting Asssumoptions 


e STM can read on atomic level 


® Potential to write on atomic scale 
— 25A - 250A bit size 


O Possibility of 10"? - 10" bits/sq. In. _ 
— 10'2.10'bits/sq.in. 


O Demonstrated Read/Write/Erase Feasibility \ 
- Voltage Puise to Write and Erase | 


- Tunneling Current to Read 


Storage Requirements 


' Induce topographical features 
on metallic glass | 


_ STM Current DASD >» (Rh,.Zr7,) 
¢ Areal Density-Mb/sq.in. 108-108 62 Wee 
Ke 350A features Induced (melting?) 
¢ Oata Rate . tad t~ 1s 
Writing speed 100ns Wns . 


Reading speed Z 30ns ; 
¢ Access Rate ? i2ms j 
e Signal-to-Noise ? 30dB i 
© Reliability | 
lifetime ? 7 years 
error rate 9 . 19°'2 
2 eee XENON ATOMS spell owt the name of the researchers’ corporate benefactor. 
Disk Roughness ? 100A 
ac gal Head Dimensions ? 10% 


Types of Limits 


Limit a Concerned With Determined By 


FUNDAMENTAL _ Whatis possible Laws of Nature 
TECHNOLOGICAL What is practical Man’s Ingenuity 
ECONOMIC What is profitable Costs and Markets 


SUBJECTIVE Whatis acceptable -— Personal Preferences 
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